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FOREWORD 

This document presents results of an investigation by 
Lockheed Missiles & Space Company, Inc., Huntsville Re- 
search & Engineering Center under Contract NAS8-26668 for 
the National Aeronautics and Space Administration, George 
C. Marshall Space Flight Center. The work was undertaken 
to perform feasibility studies of wake vortex monitoring 
systems for air terminals. 

The work was performed during the period January 
1971 through July 1972. The NASA Contracting Officer's 
Representatives (COR) for this contract are Mr. Edwin A. 
Weaver and Mr. R.M. Huffaker of S&E-AERO-AF. 
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ABSTRACT 

Wake vortex monitoring systems, especially those using laser Doppler 
sensors, were investigated under this contract. The initial phases of the 
effort involved talking with potential users (air traffic controllers, pilots, 
etc.) of a wake vortex monitoring system to determine system requirements 
from the user’s viewpoint. These discussions involved the volumes of airspace 
to be monitored for vortices, and potential methods of using the monitored 
vortex data once the data are available. A subsequent task led to determining 
a suitable mathematical model of the vortex phenomena and developing a mathe- 
matical model of the laser Doppler sensor for monitoring the vortex flow field. 
The mathematical models were used in combination to help evaluate the capa- 
bility of laser Doppler instrumentation in monitoring vortex flow fields both 
in the near vicinity of the sensor (within 1 kilometer and at long ranges (10 
kilometers). A number of computer simulations were made which scanned 
the modeled laser Doppler sensor through the modeled flow field and plotted 
the sensor output along with the actual (simulated) flowfield velocity as a 
function of position with respect to the vortex core. These plots indicated 
that a one -dimensional laser Doppler sensor will be satisfactory for viewing 
a vortex at large angles to the vortex axis. The plots also indicated that a 
two-dimensional sensor having relatively coarse resolution along each sensor 
optic axis will provide satisfactory vortex data to a range of 10 kilometers. 
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Section 1 

INTRODUCTION AND SUMMARY 

1.1 THE WAKE VORTEX PROBLEM 

Accidents associated with wake turbulence have become a significant 

reported proportion of the total aviation accident profile within the last 

fifteen years. In fact, within the i960 to 1970 period, 128 general aviation 

* 

and one air carrier accidents were attributed to wake turbulence. Approxi- 
mately one-half of these accidents resulted in fatalities. 

On 1 March 1970 the Federal Aviation Administration (FAA) implemented 
separation criteria for the express purpose of separating aircraft from wake 
turbulence. These criteria categorized as "heavy jets" those aircraft whose 
takeoff weights are greater than 300,000 pounds. Aircraft not in this category 
and which are radar vectored are required to follow the "heavies 11 by a mini- 
mum separation of five miles instead of the previous minimum of three miles. 
This introduction of another set of separation standards has a direct effect on 
terminal traffic in the arrival, departure and ground phases of operations 
(Fig. 1-1). 


The above paragraphs mention two aspects of the wake turbulence 
problem: 

• Safety — loss of life and equipment, and 

• Terminal capacity — delays in handling aircraft due to 
imposed separation requirements. 


The safety aspect of the wake turbulence problem has been improved 
significantly with the new (1970) separation standards. However, potentially 


Thomas, David D., President, Flight Safety Foundation, Inc., Presented at 
FAA Turbulence Symposium, Washington, D.C., March 1971. 

Krupinsky, E., FAA Air Traffic Service, Presented at FAA Turbulence 
Symposium, Washington, D.C., March 1971. 
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Fig. 1-1 - Aircraft Minimum Separation Requirements Influenced by Wake Turbulence 
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damaging vortices can be generated from moderately sized aircraft (Ref. 1-1) 
not explicitly covered by the new separation standards. Also, the imposed 
increased separation definitely decreases the traffic capacities of busy 
terminals during peak traffic periods, 

A goal of the authorities during the decade of the 1970s is to: 

• Increase air terminal capacities, 

• Maintain or improve terminal safety. 

An effective wake turbulence monitoring system would provide a means to 
achieve this goal without purchasing new real estate and building new runways. 
A prior knowledge of wake vortex conditions in the terminal area would allow 
aircraft to be spaced according to the prevailing vortex conditions. Thus, 
separation minimums would be increased during atmospheric conditions when 
vortices remain in the corridor and would be decreased when the vortices 
rapidly dissipate within or depart from the corridor. Air terminal capacities 
would be increased during all periods except when vortices remain in the cor- 
ridors as a hazard to arriving or departing aircraft. Safety could be improved 
by regulating spacing according to the vortex transport and decay conditions 
instead of relying on the standard separations of today which are independent 
of meteorological conditions affecting the vortex hazard. 

1.2 WAKE MONITORING SYSTEM ECONOMICS 

The fact that major air terminals are built near major population centers 
dictates that major air terminals be constructed on very expensive real estate. 
The locations of John F. Kennedy International, O’Hare, Los Angeles Inter- 
national, Atlanta, etc., exemplify this. Major air terminals will continue to 
be located on very expensive real estate until ground transportation improves 
by several factors, since air terminals are located for the convenience of the 
populace . 
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As air traffic increases in the coming years, improving the utilization 
of these major air terminals becomes increasingly important. One means of 
increasing the capacities of major air terminals without purchasing additional 
real estate for more runways is to decrease the separation minimums for land- 
ing and departing aircraft. These decreased separation standards must be 
accomplished without simultaneously degrading safety. A logical first separa- 
tion standard to attempt to decrease is the last significant separation standard 
imposed (1970) — the minimum IFR separation for n non-heavies M following 
"heavie s , M 

Figure 1-2 depicts the projected increases in landing fees which could 
be collected at major air terminals through the use of effective wake turbulence 
monitoring systems. The effective wake turbulence monitoring system would 
allow the separations behind the n heavies H to be safely decreased in specific 
air corridors when the wake turbulence hazard was minimal. Backup data 
which were used to derive the increased -landing -fee data of Fig. 1-2 are de- 
picted in Fig. 1-3. The number of hourly operations (landings or departures) 
during peak traffic hours and the approximate percentages of "heavies" using 
these terminals were obtained from discussions with air traffic control (ATC) 
tower personnel at these facilities. From these data the increase in airport 
capacities resulting from reducing separations behind "heavies” from five 
miles to three miles was projected. The capacity increase was multiplied by 
the typical landing fee to determine the increase in landing fees per hour at 
each terminal. The hourly increase in landing fees was multiplied by four 
hours per day and 250 days per year to project the annual landing fee collection 
increase through utilization of an effective wake turbulence monitoring system. 

Another economic boost from an effective wake turbulence monitoring 
system is the savings to airlines from the reduction in holding delays at air 
terminals (Fig. 1-2). Backup data for these statistics arepresented inFig. 1-4 
which were projected from surveys made from the operations of three major 
air carriers. 
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Fig. 1-2 - Economic Benefits from Wake Turbulence Monitoring Systems 
at Major Air Terminals 
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Fig. 1-3 - Projected Increased Air Terminal Landing Fees with Effective Wake Turbulence 
Monitoring Systems 
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1.3 WHAT IS A WAKE TURBULENCE MONITORING SYSTEM? 


A wake turbulence monitoring system performs the following three 
functions: 

• Detects and/or predicts current locus of aircraft trailing 
vortices in the air corridors of an airport, and 

• Forecasts the future hazard to aircraft trailing (arriving or 
departing) the generator aircraft, and 

• Interfaces with traffic controllers and/or pilots to prevent 
a hazardous vortex encounter. 


The first of the above functions deals with locating vortices in the air 
corridors leading to and from an air terminal. This might be done with an 
active device that transmits radiation which is reflected or altered by the 
vortex presence, by a passive device which detects radiation emitted by the 
vortex itself, or by a predictive system which predicts the vortex locus from 
current meteorological data. Active and passive techniques which might be 
considered for a wake turbulence detector along with some of their character- 
istics are listed below: 

TECHNIQUES FOR PERFORMING FIRST FUNCTION: 

"Detect and/or Predict Current Vortex Locus" 


Radiation 

Frequency 

Passive Technique 

Active Technique 

0 to 100 Hz 

Detection of pressure wave 
generated by vortex. 
Characteristic: Low ampli- 
tude signal. Detects vortices 
in ground effect. 


100 Hz to 20kHz 

Detection of audible sound 
generated by vortex. 
Characteristic: Low ampli- 
tude signal compared to 
background noise. 

Detection of sound back- 
scattered by vortex. 
Characteristic: Affected 
by noisy acoustical en- 
vironment. Is being tested 
at air terminal. 
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Radiation 
F requency 


Passive Technique 


Active Technique 


Microwave Measure incoherent radiation 

associated with elevated 
vortex temperatures. 
Characteristic : Low energy- 
emission at RF wavelengths. 


Infrared, Visible Measure radiation associated 
and Ultraviolet with elevated vortex tempera- 
ture. Characteristic : Back- 
ground radiation and at- 
mospheric effects upon 
transmission. 


Detection of RF energy 
scattered by vortex tur- 
bulence of similar wave 
length. Characteristic : 
Vortex is poor RF re- 
flector. Much transmitted 
energy required. 

Detection of backscattered 
energy by aerosol particles. 
Characteristic : Fundamental 
questions are answered for 
IR frequencies. Method 
should be evaluated in air 
terminal environment. 


While the above chart depicts active and passive vortex detection devices, the 
totally predictive system would rely on none of the above devices but would 
utilize local meteorological data and aircraft coordinates to predict the pres- 
ence of the vortices. The major portion of this report is devoted to wake 
turbulence monitoring systems which use active vortex detectors operating 
at infrared frequencies. 


The second of the above functions — "Forecasts the future hazard to 
aircraft. . . " involves predicting vortex movements and their decay and fore- 
casting the hazard to trailing aircraft. This forecasting function would be 
performed either by a minicomputer or a module of the existing computer 
system at the air terminal and would depend upon inputs from the detector 
sensors and/or meteorological inputs of the above paragraph. Air traffic 
controllers might also supply inputs to the forecasting subsystem of the 
turbulence monitoring system. The forecasting portion of the wake turbu- 
lence monitoring system would initially rely upon a very simple vortex model, 
but would likely become more sophisticated as more knowledge is gained about 
the vortex decay and transport phenomenon. 


The third function of the monitoring system — "Interface with traffic 
controllers, pilots, etc. . . " is equally as important as the first two. For a 


3 
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wake monitoring system to be effective, there must be a means of effectively 
communicating with the personnel involved. Communication with what per- 
sonnel, and regarding what information should be relayed is yet unresolved. 

The air traffic controllers (ATC) prefer that the system communicate directly 
with the pilots and leave ATC out of the loop. However, if the system is to 
increase terminal capacities, aircraft separations must be decreased, and 
the controller becomes involved since he directs minimum spacing. The con- 
troller also becomes involved if a pilot decides to abort a landing, since not 
only the aborting aircraft is affected, but also the whole landing sequence. 

1.4 USER REQUIREMENTS FOR A WAKE TURBULENCE MONITORING 

SYSTEM 

Lockheed interviewed potential users of a wake turbulence monitoring 
system in the early phases of the subject contract. These potential users 
included local and departure/arrival air traffic controllers at Atlanta Municipal 
and John F. Kennedy International airports, Delta Airline pilots in Atlanta, and 
district FAA officials in Atlanta. Lockheed felt that inputs from potential users 
in the early phases of system concept were essential for the development of an 
effective wake turbulence monitoring system. All of these persons recognized 
the need for a wake turbulence monitoring system and were most cooperative 
in helping to develop criteria for such a system. Figure 1-5 is a summary of 
user requirements developed from interviews with the above personnel. 

The requirements developed from interviewing the potential system users 
included the "volume of interest" and the "communications/command process" 
of Fig. 1-5. The volume of interest includes the air corridor to the altitude 
above which an aircraft encountering a vortex would have a high probability 
of safe recovery. Other factors affecting this recovery altitude include the 
different aircraft configurations for arrival and departure and the configuration 
effect upon the likelihood of the aircraft stalling or becoming uncontrollable 
upon encountering a vortex. Since the aircraft is in an accelerating mode and 
is normally climbing rapidly, less range is probably required for a system 
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Volume of Interest 

Range (beyond end of runway) 

Two miles — Departure 
Five miles — Arrival 
Horizontal (Azimuth) Angle — 

2.5 degrees to either side of runway centerline plus area to 
accommodate vortex drift between scans 

Vertical (Elevation) Angle — 

2 degrees to 15 degrees —Departure 

2 degrees to 3 degrees —Arrival 

Communications/Command Process 

Traffic Controller’s Viewpoint — Controllers are very busy; 
system should communicate directly with pilots. 

Pilots Viewpoint — Pilots are busy in terminal area; let 
controllers direct spacing interval. 

Compromise — Controllers would be involved anyway if pilots 
perform evasive maneuver; therefore, controllers should 
probably control spacing interval according to prevailing vortex 
conditions and pilots take only evasive action when required. 

* 

6 to 1 1 degrees would cover approximately 90% of cases 
ILS and VASI glide paths nominally 2.0 jf0.5 degrees 


Fig. 1-5 - User (Pilots, ATC, etc.) Requirements 
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monitoring departure corridors than is required for monitoring arrival cor- 
ridors where aircraft are in near stall configurations coming in at a very- 
shallow (^2^ deg) angle. The horizontal width of these volumes was projected 
from the instrument landing system (ILS) beam width. Some additional width 
should be included in the volume to prevent vortices from drifting into the flight 
corridor between successive sensor scans of the volume. 

The communications/command process as discussed in Section 1.2 in- 
volves relaying the information obtained by the system to the user in an effective 
manner. The differences in the two depicted viewpoints must be resolved for 
the development of an effective system. 

1.5 SYSTEM CATEGORIES 

Wake turbulence monitoring systems are divided into three basic cate- 
gories to facilitate discussion in this report. They are: 

• Predictive 

• Hybrid 

• Detective 

Each of the categories of wake turbulence monitoring systems performs the 
three functions described in Section 1.4: (1) detects and/or predicts current 
vortex locus; (2) forecasts future hazard; and (3) interfaces with ATC and/or 
pilots to prevent hazardous encounter. 


The basic similarities and differences in the three types of systems 
involve the first two functions. They are depicted in the table on the following 
page. As seen in the table, all three types of systems perform measurements 
and forecast predictions. The different nomenclature for the three categories 
is derived from the relative degrees of prediction and measurement (detection) 
performed by the three basic systems: 
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Method of Performing System Function 

System Function 

Predictive 

System 

Hybrid 

System 

Detective 

System 

1. Determine Current 
Vortex Locus 

Predict 

Measure and 
Predict 

Measure 

2. Forecast Future 
Hazard 
a. Determine 

meteorological 

conditions 

Measure and 
Predict 

Measure 

Measure 

b. Determine 
future vortex 
locus 

Predict 

Predict 

Predict 


• Predictive : Measures meteorological data with presently avail- 
able instrumentation; forecasts additional meteorological data; 
forecast future vortex hazard. 

• Hybrid : Measures vortex movement in several vertical planes; 
forecasts vortex movement between planes ; forecasts vortex 
future hazard. 

• Detective : Measures vortex movement throughout the air cor- 
ridor of interest; forecasts future hazard. 

All three systems both measure and forecast. The detective systems 
have more sophisticated instrumentation and therefore provide higher quality 
data to the forecasting computer. Therefore, the detective systems are likely 
to be more expensive to develop; yet they provide a more accurate forecast of 
the future vortex hazard. The predictive system relies upon presently avail- 
able, limited meteorological data and therefore would provide a less accurate 
forecast. The hybrid system relies upon viewing the vortex at a limited num- 
ber of loci and therefore would provide data to the computer of a higher quality 
than that of the predictive system but of lower quality than that of the detective 
system, 
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Since the detective system was said to provide the most accurate hazard 
forecast of the three system configurations discussed, the installation of a de- 
tective system would allow air terminal capacity to be increased to its maxi- 
mum. The predictive system, which would provide a less accurate forecast, 
would require a greater safety margin than would the more accurate forecast. 
The greater safety margin requirements of the predictive system would de- 
crease the effectiveness of the system in increasing air terminal capacities. 
The effectiveness of the hybrid system will likely lie between that of the pre- 
dictive and detective systems. These facts are illustrated in the following 
table. 



Predictive 

Hybrid 

Detective 

Sensor Sophistication 
(Therefore Expense) 

Least 

Median 

Most 

Hazard Forecast 
Reliability 

Least 

Median 

Most 

Required Safety 
Margin 

Most 

Median 

Least 

Increase in Terminal 
Capacity 

Least 

Median 

Most 

Complexity of Integration 
into Existing System 

Least 

Median 

Most 

Application to Other Air- 
Terminal Problems 


Wind Shear 
Monitor 

Wind Shear 
Monitor 

Development Time 

Least 

Median 

Most 


A suggested path of development of the wake turbulence monitoring system 
is to field the predictive system first, while simultaneously developing tech- 
niques for the hybrid and/or detective systems. Next, additional sensors are 
incorporated into the predictive system, which evolves into a hybrid or de- 
tective system. This results in immediate improvement in air terminal capa- 
bilities with further improvements forthcoming as the systems become more 
sophisticated. 
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Figure 1-6 depicts an illustrator r s concept of each of the three systems 
to be discussed. The left-hand system is a detective system, the right-hand 
system is a predictive system and the middle system is a hybrid system which 
uses some of the components of both the detective and predictive systems. 

• Predictive System 

The predictive system conceptually illustrated in Fig . 1-6 utilizes available 
meteorological data as input to a minicomputer which predict s lengths of time vor- 
tex elements will remain in the pertinent air corridor s . The available meteor- 
ological data are used to calculate wind conditions aloft which are in turn used 
to calculate displacements of vortex elements in the air corridors. The computer 
output — a prediction of how long vortex elements will remain as hazards in the 
pertinent air corridors — is displayed to local controllers who in turn either in- 
crease or decrease minimum separation requirements according to the prevail- 
ing vortex conditions. The controller action is taken before aircraft enter the 
final approach corridors or, in the case of departing aircraft, before the air- 
craft enter the runway. This is the least expensive of the three systems. It re- 
quires only a minicomputer, display, and existing meteorological measurement 
facilities. The predictive system is discussed in Section 2 of this document. 

• Detective System 

The detective system briefly illustrated in Fig. 1-6 consists of, for ex- 
ample, a pair of laser Doppler velocimeters (LDV)for the runway with a range 
of 8 to 10 kilometers. These LDVs are tied together through a computer which 
directs their scanning and processes their outputs. Formatted data from the 
computer are transmitted directly to the cockpit and displayed to the pilot as well 
as to the controllers in real time. The velocimeters provide a real-time two- 
dimensional picture of the vortex condition in the air corridor directly to the 
aircraft pilot. The pilot can then take appropriate action to avoid vortex hazards. 
Similar displays would be provided to air traffic controllers to coordinate the 
actions of the pilots. This is the most expensive of the three illustrated systems, 
but it provides the highest resolution, longest -range data. The detective system 
is discussed in Section 3. 
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• Hybrid System 

The hybrid system depicted conceptually in Fig. 1-6 is composed of 
certain elements of both the predictive and detective systems. Development 
time and implementation costs of the hybrid system will be less than for the 
detective system, yet longer and more costly than for the predictive system. 
The hybrid system utilizes for example, a single component LDV which scans 
the surface of a cone overhead to provide wind profiles for the prediction of 
vortex element movements. A single -component LDV located several kilo- 
meters beyond the end of the runway is also used to scan a vertical plane 
normal to the runway axis to locate vortices above the velocimeter and to 
determine the sizes and transport velocities of the vortices as they pass 
through this plane. Acoustic sensors may be feasible for locating vortices 
at lower altitudes. Section 4 of this document presents a detailed discussion 
of the hybrid system. A typical hybrid system will interface Doppler sensor 
outputs into a minicomputer which will predict the movements of the vortices 
in the air corridors from the available data. The computer output will be dis- 
played to the air traffic controllers who in turn control the traffic spacing. 

The prediction of the vortex movement by the hybrid system is of finer reso- 
lution than the prediction by the predictive system since more pertinent input 
data are supplied to the computer. Thus air traffic spacing can be controlled 
to a finer tolerance (i.e., smaller minimum separations) than with the purely 
predictive system. Section 4 presents a detailed discussion of the hybrid ap- 
proach to the wake turbulence monitoring problem. 

Section 5 of this document, the conclusion, presents an attempt to answer 
four basic questions about the wake turbulence monitoring system. The four 
questions are: 

• Is a wake turbulence monitoring system feasible or within the 
near state of the art? 

• Would a wake turbulence system be cost-effective to develop, 
operate, and maintain ? 

• What type system should be developed? 

• What are the logical steps to develop the system? 
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The primary sections of the report are followed by three appendixes. 
Appendix A describes the vortex models chosen by Lockheed to be used in 
system simulation studies. Development of the programs used for these 
simulation studies is outlined in Appendixes B and C. 
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Section 2 

PREDICTIVE WAKE VORTEX MONITORING SYSTEM 
2.1 INTRODUCTION AND SUMMARY 

t 

The five -mile minimum spacing regulation for an aircraft lighter than 
300,000 pounds following a "heavy" (> 300,000 pounds) is presently enforced 
for all flying conditions, regardless of weather phenomena. However, our 
present knowledge of the vortex phenomenon indicates that for certain meteor- 
ological conditions the aircraft wake is either swept out of the flight path or 
breaks up internally in less time than that imposed by five -mile minimum 
spacing. If these conditions could be accurately predicted from measurements 
of local meteorological phenomena, spacing could be decreased during appro- 
priate meteorological conditions. A predictive system based on considerations 
of local meteorological conditions, air terminal layout, aircraft involved, etc., 
could be beneficial to major air terminals. The initial wake turbulence pre- 
diction system which could be operational within a year or two could reduce 
spacing requirements from five miles to four miles or three miles for, per- 
haps, 50% of the operational conditions. As the state -of -knowledge of the 
vortex phenomenon improves, the spacing reduction might be extended to 60 
or 70% of the operational conditions. 

As previously discussed in Section 1, potentially damaging vortices can 
be generated from moderately sized as well as "heavy" aircraft. With a turbu- 
lence prediction system, more adequate warning of wake turbulence conditions 
could be given to small aircraft operators following a moderately sized aircraft 
in both approach and departure configurations. 

A wake turbulence predictive system using existing instrumentation at 
the air terminal coupled with computer capability for predicting the vortex 
hazard (or lack of hazard) and a simple communications system could thus 
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provide an increased capability for handling more traffic without degradation 
of safety. 

A typical wake turbulence predictive system is depicted in block dia- 
gram form in Fig. 2-1. The wake turbulence predictive system consists of 
several subsystems which are discussed below. Section 2.2 discusses the 
vortex phenomenon and its behavior. This will include definition of all param- 
eters which enter into forecasting and displaying aircraft wake vortex strength, 
persistence and transport in terminal area operations. The acquisition of some 
of these parameters such as wind speed and direction at various altitudes and 
temperatures is illustrated in Fig. 2-1, Design of a wake turbulence predictive 
system is outlined in Section 2.3. This system will be the software utilized by 
the computer illustrated in Fig. 2-1. Section 2,4 discusses the information 
acquisition, storage and processing subsystem that will execute the decision- 
making process. The subsystems (as illustrated in Fig. 2-1) will likely con- 
sist of; 

1. The sensors together with a means for entering their data 
into the computer; 

2. Data from the aviation (weather) forecaster together with a 
means for entering this data into the computer, and 

3. The computer itself. 

The communications and action command process is discussed in 
Section 2.5. This process is illustrated typically as an air traffic con- 
troller display which forecasts local vortex conditions. The controller uses 
this display and knowledge of aircraft in the vicinity to project spacing re- 
quirements for departing and arriving aircraft. He can also enter data into 
the computer such as runways which are closed to heavy aircraft, aircraft 
departure angles, etc. 

2.2 THE TRAILING VORTEX PHENOMENON AND BEHAVIOR 

This subsection discusses pertinent properties of the trailing vortex 
system, including its generation, structure dissipation and transport. 
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Fig. 2-1 - Typical Wake Turbulence Predictive System Block Diagram 
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• Mechanisms of Formation of Aircraft Trailing Vortices 


The general principle by which an aircraft generates lift is directly re- 
lated to Newton's Second Law (F = d(mV)/dt), where the lifting force is a result 
of a downward acceleration of air. The process is usually described in theo- 
retical analyses as producing circulation about the wing. This produces a con- 
tinuous vortex sheet which is shed from the trailing edge of the wing. Because 
the wing is of finite span length, the high pressure air of the lower surface 
flows around the wing tip into the region of low pressure on the top of the wing; 
thus, a vortex core is created with its center near the trailing edge of the wing 
tip. The strength of these tip vortices is directly related to the lift produced 
by the aircraft wing. 

The existence of these trailing vortices has been known for many years. 
Two decades ago Spreiter and Sacks (Ref. 1 ) presented what is still probably 
the most significant analysis of trailing vortices. Their approach, which was 
based upon a solution of the linearized potential equation, revealed that the 
distance required for the trailing vortex sheet to become established is 


e^ 

c 


0.28 


A b 

C L C 


( 2 . 1 ) 


where 

e/c = roll up distance in chord lengths 
c = wing chord 
A = wing aspect ratio 
b = wing span 

= aircraft lift coefficient 


For a typical takeoff of a C-5A, this means that the trailing vortices will be 
completely formed approximately 300 feet behind the wing. 
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• Vortex Structure 

A vortex consists of a rotating region of irrotational flow and a smaller 
central core of rotational flow. The core may be thought of as rotating as if 
it were a solid cylinder. Spreiter and Sacks (Ref. 1 ) estimated the diameter 
of this core by equating the total kinetic energy of the fluid inside the vortex 
core to the induced drag of the wing. This resulted in a core diameter of 
approximately 16% of the wing span for an elliptic wing load distribution. 
Although the assumptions involved in this computation are rather dubious 
and believed to be invalid during the decay phase of the vortex, it does pro- 
vide an indication of initial core size. Core sizes are discussed further in 
Appendix A, where recent vortex models are considered. 

Just as lift of a wing is generated by creating downwash, the drag of a 
wing results in a change of momentum of the air in the axial direction. Ex- 
perimental data (Refs. 2 and 3) indicate that relatively high axial velocities 
exist in the center of the vortex. Quite possibly, most of the axial momentum 
of the vortex is restricted to this vortex core. However, the nature and origin 
of this axial component of velocity in the core of the vortex is apparently an 
area requiring additional investigation and comprehension. One explanation 
(Ref. 4) is that it is due to the pressure buildup along the vortex core due to 
vortex decay because of the coupling between core pressure and tangential 
momentum. 

Outside the core the tangential (transverse) velocity, at a radial distance 
(r) from the center of the vortex, is estimated from simple theory as 

r Q 

v n - ~ , where J f = magnitude of circulation (2.2) 

6 27 rr’ o j .. . o 

around the wing (Eq. 2.5) 


2-5 


c 


LOCKHEED - HUNTSVILLE RESEARCH & ENGINEERING CENTER 



LMSC-HREC D 225936 


or in terms of aircraft weight (W), velocity (V), and wing span (b), and density 
(P) 

2 W 

v e = — pTbr • (2 ' 3) 

7T 

By assuming the maximum velocity occurs at the radius of the core which is 
taken as 8% of the wing span, tangential velocities of up to 65 ft/sec might be 
expected to exist in a fully rolled -up vortex generated by a C-5A aircraft 
during takeoff. Recent flight measurements (Ref. 5)however, have indicated 
smaller cores and peak velocities considerably in excess of this figure, sug- 
gesting the simple model described above is inadequate for practical purposes. 
Recent models are discussed further in Appendix A. 

Because of the rotational nature of the flow, the tangential velocity in 
the vortex core is more difficult to describe. Lamb (Ref. 6) presented the 
following approximate relation for computing the tangential velocity through- 
out the vortex: 

0 2ir r 

where V - kinematic viscosity, t = time. 

This expression not only includes the viscosity effect in the core, but allows 
this region to increase with time. For large values of r, the viscous con- 
tribution to this expression becomes negligible; hence the irrotational velocity 
is retained. For practical purposes it has been found necessary to define an 
eddy viscosity rather than the kinematic viscosity since the vortices are prob- 
ably invariably turbulent. Such an approach is due to Squire (Ref. 7). 

• Vortex Strength 

Spreiter and Sacks explain that n the strength of one of the trailing vortices 
is equal to the sum of the strengths of all the vortices shed from one -half of the 
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wing; hence it is equal to the magnitude of the circulation P around the wing in 
the plane of symmetry." By assuming an elliptic lift distribution over the wing 
(an idealized distribution of high lift devices), this circulation is 

T 0 = 2bY^ . (2.5) 


The circulation (i.e., vortex strength) may also be written in a slightly differ- 
ent form by assuming that the lift generated by the wing is approximately the 
same as the aircraft weight, W and thus 


r 


4W 

7Tp V b 


(2.6) 


where p and V are the ambient density and aircraft velocity, respectively. 


• Vortex Drift 


The vortex cores are established slightly inboard of the wing tips at a 
position which is stabilized by their mutual interaction. The stable position 
of the two vortex cores is a function of the lift distribution over the wing, but 
a typical stabilized position, for the case of elliptical lift distribution, is at 

t 

a distance b, as given in Ref. 1 as 

b = ( 7T / 4)b (2.7) 

which is approximately 80% of the wing span. 


The vortex system also drifts downward relative to the aircraft flight 
path at a vertical velocity, Z, determined as follows: 



( 2 . 8 ) 


Assuming an elliptic lift distribution and that the wing lift equals the aircraft 
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weight, the formula for the settle rate of the vortex system is restated in 
aircraft terms as follows: 


_8_ W 

-i 2 ' 

7T pvb* 


(2.9) 


For a typical C-5A takeoff, this vertical velocity will be approximately 7 ft/sec. 

By neglecting any buoyancy effects due to density gradients in the atmos- 
phere, these vortices are expected to continue this downward drift at essen- 
tially a constant rate until they either decay or approach the ground. As the 
vortices approach the ground the vertical velocity is reduced and the system 
starts to spread horizontally. This ground effect may be computed by con- 
sidering the "mirror image" concept which is commonly used in potential 
solutions. In this approach the ground plane may be thought of as a plane of 
symmetry with two image vortices of equal strength drifting upward at identical 
velocities toward the ground plane. The real vortices will continue to drift 
downward until the influence of the image vortices slows and eventually stops 
their descent. The downward momentum of these vortices is conserved and 
serves to spread the vortices horizontally at approximately the speed of their 
original descent. The stabilized height of these vortices is established at b/2 
above the ground. Vortices generated closer to the ground will rise to this 
height while spreading horizontally. 


Ambient winds also significantly influence vortex drift. A uniform cross- 
wind forces the vortices in the direction of the wind. For a crosswind from 
right to left, the system of vortices could be thought of as stationary relative 
to the crosswind and moving to the left relative to a fixed point on the ground. 

If the magnitude of the drift velocity of the vortices is equal to the wind speed, 
then the right vortex will become stationary relative to the ground. This vortex 
would remain stationary until it has lost its effect through normal decay. Should 
this stationary point be over a runway, the hazard to other aircraft attempting 
to take off or land is obvious. 
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A more general treatment of vortex motions is now given, considering 
the behavior of one vortex pair (Ref. 8). 

In the analysis that follows, the vortex-induced velocities are assumed 
not to vary with time; i.e., the effects of eddy viscosity on the circumferential 
velocity distribution in the vortex have been neglected, as these effects are 
considerable only in the region adjacent to the vortex core and the mutually 
induced vortex motions are dictated by the velocities prevailing several core 
diameters away from the vortex center. This can of course be easily modified 
to the decaying vortex case. For typical aircraft under approach conditions 
Fig. 2-5 suggests that it will be several minutes before the relevant velocities 
decay to a harmless state. 



Fig. 2-2 - Vortex Drift Near the Ground 


The analysis follows the methods of classical hydrodynamics. Let the 
trailing vortex pair be located a distance 2y apart, and at height z above the 
ground, then the vortex system with its image may be represented as stated 
in Fig. 2-2. 


Now the circumferential velocity, v, due to a single vortex of strength, 

K, is 


v = K/2 j r r 


( 2 . 10 ) 
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where r is the distance from the vortex axis. 

If we consider the velocities induced at Q by the remaining vortex and 
the two image vortices, and resolve velocities (i) horizontally and (ii) vertically 
we have 

(i) 


and 




therefore, 


and 


Y = 


— I—*—) 

lirz L 2 + z 2 ' 


( 2 . 11 ) 


z = 


I sl J \ 

ly 2 +z 2 ' ' 


( 2 . 12 ) 


Combining Eqs. (2.11) and (2.12) 


i 

• dz 


3 

_ X_ 

3 

z 


or, 


dz 

3 


dx 

3 

y 


and hence, 


J_ + J_ - A 
2 2 ’ 
V z 


(2.13) 


The term A may be found by substituting y Q and z q for y and z, respectively. 
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Equations (2.11) and (2.12) may now be rewritten as 


K 

4 ttAz 3 
-K 

4jrAy 3 


By eliminating z from Eq. (2-11) and rearranging: 


or, 


Y 


K 

4?rA 



4n 

K 



dy 


dt 

A • 


Hence, it can be shown that, 


4ir 

AK 



= t + B 


B may be found by substituting y Q for y at t = 0. 


Rearranging, 

2 

y = 


2[64tt 2 /AK 2 + A(t 


64jt 2 /K 2 


j, + 


6 45r 2 /A 2 K 2 (t + B) 2 + l' 


- 1/2 


and, by symmetry, z is given by: 


2 

z 


2[64ff 2 /AK 2 + A(t +B) 2 ] ( / 64?r 2 

64tt 2 /K 2 ( 'A 2 K 2 (t+B) 2 
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In each of these expressions, the upper alternative sign applies when t < B f 
and the lower alternative sign applies when t > B. 

Hence, the vertical and horizontal displacements of the vortex cores 
can be calculated as a function of time to give the theoretical vortex positions. 

The effe ct o f wi nd on these results can be obtained by simply super- 
imposing the wind velocity components on the above velocity components. 


As an illustration of the application of the foregoing, it is interesting to 
recall test results obtained on the behavior of trailing vortices near the ground 
(Ref. 8). The test configuration is as shown in Fig. 2-3. Behavior of the two 
smoke -filled vortices generated by a Hunter aircraft were monitored over a 
measurement plane by both ground and airborne stations. The motions of these 
vortices were mapped and compared with the results of the previous theory after 
a measured wind velocity was incorporated. Typical examples of data obtained 
are shown in Fig. 2-4. From the good agreement between theory and practice, 
one must deduce that knowledge of the wind velocity is sufficient to determine 
the path of the vortex to good accuracy. The basis of an approach would 
therefore be to map the wind field in the pertinent regions either by multiple wind 
sensors or by a limited number of wind sensors plus inference (e.g., a 1/7 power 
law for the variation of wind strength with height was found to be satisfactory in 
the aforementioned study). A model for the spatial turbulence pattern might also 
be incorporated to yield more dependable predictions. 


• Vortex Decay and Breakdown 

The process by which vortices decay and eventually break down is the 
least understood phenomenon associated with vortices. Under certain con- 
ditions a vortex will completely disintegrate a few seconds behind the aircraft; 
whereas under different conditions a vortex might persist for more than five 
minutes after generation by the same aircraft. 
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Fig. 2-3 - Diagrammatic View of Test Area 

2-13 


LOCKHEED - HUNTSVILLE RESEARCH & ENGINEERING CENTER 


LMSC-HREC D225936 




Crosswind = -3.5 ft/s (Ref. 8) 

Fig. 2-4 - Comparison of Measured and Theoretical Vortex Positions 
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There are several theories which attempt to isolate the mechanisms 
which cause decay. The previously presented time -dependent expression v^ 
represents the original approach. This method assumes that the viscous core 
grows in size and the core velocities decrease with time. However, this ex- 
pression predicts a growth of the core which is much slower than that which 
has been observed for trailing vortices. The explanation of this core growth 
was that turbulence existed in the core in addition to viscous effects. This 
would result in much higher shear forces than those of the expression for v^. 
This turbulence could easily arise from atmospheric turbulence and from the 
turbulent wake of the airplane. A viscosity term known hs "eddy viscosity" 
was derived empirically in an attempt to account for the turbulence effect on 
core growth. Wetmore and Reeder (Ref. 9) present an alternate expression 
for Vq which is based on such an eddy viscosity in calm air; 


v 


0 


£o 

2ir r 


1 - exp 


i—2—i 

r 


(42b 2 + 1.2c r t) (10" 4 ) 



(2.16) 


where c is the mean aerodynamic chord of the aircraft wing. This expression 
represents a maximum condition for vortex persistence. More rapid vortex 
decay would be expected if additional turbulence levels were created by such 
factors as: 

• Cross winds 

• Convection currents 

• Friction resulting from vortex-ground interaction 

• Irregular span -load distributions 

• Jet exhaust entrainment. 


One obvious disadvantage of this method of predicting vortex decay 
stems from its empirical nature and the neglect of the effect of the items 
listed above. A case when none of these items would be expected to influ- 
ence the turbulence level would rarely occur. On the other hand, this method 
is easy to compute. If additional tests or more thorough analysis of existing 
test data can provide empirical constants to include the effect of the above 
listed items on the vortex decay, then this method can provide a practical 
engineering approach toward predicting vortex decay and breakdown. 
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Several theories have been proposed for providing a quantitative analyt- 
ical description of, and isolating the mechanisms leading to vortex decay and 
eventual breakdown. Squire (Ref. 10) has suggested that breakdown will occur 
when long axisymmetric standing waves are sustained in the vortex. Another 
approach, suggested by Ludwieg (Ref. 11) and Jones (Ref. 12), suggests that 
breakdown is a result of a hydrodynamic instability with respect to spiral dis- 
turbances. Still another idea proposed by Crow (Ref. 13) says that vortex decay 
and breakdown result from vortices undergoing a symmetric and nearly sinu- 
soidal instability. Vortex breakdown is assumed when the vortices eventually 
join at intervals to form a train of vortex rings which quickly disintegrates into 
a harmless turbulent state. 

Another entirely different approach to describing vortex breakdown was 
proposed by Benjamin (Ref. 14). This theory, which appears quite popular 
with British researchers (i.e., Hall in Ref. 15 and Harvey in Ref. 16), states 
that vortex breakdown does not result from an instability or any other infini- 
tesimal disturbance alone. Instead, he describes two basically different types 
of axisymmetric flow which are analogous to subcritical and supercritical 
flow in hydraulic open-channel flow. The transition from the supercritical 
form of the axisymmetric flow to the conjugate subcritical flow results in 
phenomena analogous to the hydraulic jump in an open channel flow. This 
transition requires an energy loss which means that in practice a region of 
vigorous turbulence is generated; hence, vortex breakdown occurs. 

Which, if any, of these theories correctly describes the conditions which 
determine the mechanism by which vortices decay and break down is not now 
known. Those derived analytically from such ideas as unstable disturbance 
propagation or the hydraulic jump analogy and sinusoidal instability certainly 
provide the most logical and scientific approach. However, the empirical 
technique using eddy viscosity probably provides the best engineering tool for 
predicting vortex decay. This approach, with perhaps a slight modification to 
include improved empirical constants obtained from recent flight tests, could 
easily be computerized and combined with a digital computer program for 
tracking aircraft trailing vortices. 
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The vortex lifetime data which should probably be utilized for immediate 
development of a predictive system is that presented by McGowan (Ref, 17). The 
curve of that reference is presented in Fig, 2-5 on which is also superimposed 
data on vortex lifetimes obtained from recent FAA/NAFEC tests on the vortices 
generated by a Convair 880. McGowan's curve conveniently yields an envelope 
to the experimental data and is probably a satisfactory state of the art curve. 
These vortex lifetimes can be considered as the time taken by the vortex to 
dissipate below some threshold peak velocity condition, the decay being due 
either to viscous effects or to vortex bursting, etc. The final model would 
probably allow for an eddy viscous decay up to the time limit given by the above 
curve at which time it will be assumed that the vortex no longer presents any 
danger . 

• Previous Efforts to Define the Vortex Hazard in the Airport Environment 

Previous efforts (Refs. 17 and 18) have mapped vortex motion in various 
wind fields to determine hazardous areas in airports for normal terminal opera- 
tions. These efforts have resulted in plots like those shown in Figs. 2-6 and 
2-7. 


In this way caution areas surrounding a particular aircraft rotation point 
have been defined as being the envelope of such a system of curves. A similar 
system of curves exists for each touchdown and rotation point and each point 
on the flight path. The hazard area then becomes the envelope of the system 
of curves in Fig. 2-6 for each point of the flight path for a given wind velocity 
history. Such hazard areas may well include a parallel runway or taxi strips 
as indicated in the example of Fig, 2-7. 

Similar considerations have been applied to the landing approach and 
climbout phases of flight. Figure 2-8 shows how a crosswind component of 
various magnitudes can modify and shift the vortex velocity pattern. From 
this, the importance of considering the possibility of the vortex drifting into 
the glide or takeoff path of a parallel runway is demonstrated. 


2-17 


LOCKHEED - HUNTSVILLE RESEARCH & ENGINEERING CENTER 




2-18 

EARCH & ENGINEERING CENTER 


LOCKHEED - HUNTSVILLE RES 


LMSC-HREC D225936 


LATERAL DISTANCE 



Wind Speeds 
0-40 Knots 


Wind directions 




♦►90 

60° 


0° 30° 


Fig. 2-6 - Vortex Movement Near Ground with Wind (Ref. 17) 


Federal Aviation Administration separation regulations and Air Traffic 
Control operations have been established using curves and plots such as these 
for guidelines. The purpose of a wake turbulence predictive system develop- 
ment would be to develop a near real-time capability — a software accented 
predictive system — for establishing accurate separation requirements to allow 
increased utilization of air terminals without degradation of safety. 

• Recent Vortex Models 

Prediction of vortex transport requires only a simplified model of vortex 
structure since one vortex moves in the field of another at relatively large 
separations, the inverse radius approximation being satisfactory under such 
circumstances. However, for aircraft-vortex interaction estimation, a more 
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Fig. 2-8 - Envelope of Vortex System Displacement with Wind 
(No Ground Effect) (Ref. 17) 


realistic vortex model is required — one that yields predictions of peak velo- 
cities and core diameters more in line with recent FAA data (Ref, 5), Lockheed, 
during this investigation, has surveyed current vortex models both empirical 
and theoretical. These current models are discussed in Appendix A. 

2.3 DESIGN OF A WAKE TURBULENCE PREDICTIVE SYSTEM 

This section discusses the design of a wake turbulence predictive system. 
The aim of this design was the simplest system capable of performing a useful 
role in vortex prediction. Also, flexibility was sought, so that sophistication 
could be added after the initial systems are operational. For the proposed de- 
sign, the vortex pairs are divided into vortex elements, each of which consists 
of a cross sectional slice of the vortex. A segment of the vortex pair is treated 
as a single element when out of "ground effect." Once in "ground effect," the 
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pair is considered as two elements — one port element and one starboard ele- 
ment. A number of such elements are used to represent vortex pairs corre- 
sponding to vortices over each runway and in each air corridor leading to the 
airport. 

The discussed system design assumes that at time t = 0 the vortex sys- 
tem is frozen along the path of the generating aircraft (see Fig. 2-9). The 
vortex elements are then allowed to settle, move in the windfield and undergo 
ground effect. Equations in the computer provide vortex element positions as 
a function of time (t = 10, 20, 30 sec, etc.) for the input wind field. 

Volumes of the air space where presence of an active vortex would con- 
stitute a hazard; i.e., the active air corridors leading to and from the airport 
and the near vicinity of the runways themselves are specified in the same co- 
ordinate system as the vortex elements. These areas are entered into the 
computer as simple line equation coordinates. For the example illustrated, 
and probably the first system design, these danger areas are not extended be- 
yond the middle marker; however, for operational type systems these would 
probably be extended at least to the outer marker. Initially, the danger zones 
are also considered to be two-dimensional. For an operational system this 
space might be reduced by including a "z M coordinate in the danger zone speci- 
fications. 

After the vortex element initial conditions, meteorlogical conditions 
and danger areas are entered into the computer, the vortex elements are 
programmed to move according to the conditions for a specified time incre- 
ment. After this displacement, the individual vortex elements are tested to 
determine if any elements are still in danger areas. The time increment and 
danger area identification are stored and the vortex elements are displaced via 
computer to their loci corresponding to the second time increment. Again, 
the time and danger areas are stored. This process is continued until all 
vortex elements have moved for time increments corresponding to their cal- 
culated lifetimes (Section 2.2). 
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Fig. 2-9 - Danger Areas for Vortex Hazard Example: Atlanta Airport 
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Next, the stored time/danger-area data are examined to determine the 
length of time taken to clear all vortex elements from each danger area. This 
ti me to c lea r all e lements should correspond to the time required for separa - 
tion between aircraft using the corridor. 

To repeat , the basic assumption of this proposed system software is that 
vortex elements at time t = 0 are spatially at the points where they were gen- 
erated. They are then allowed to drift in space for given time increments. 

Their positions at the end of each time increment are noted and compared to 
danger zones in air space which have been entered into the computer. The 
computer notes when elements are in specific danger zones and stores these 
data. After the elements have been processed for time increments equaling 
their lifetimes, the computer scans the data to determine when it would be safe 
for an aircraft to follow the wake -generating aircraft. Thus, a safe but not 
excessive separation distance is established. 

Departing and arriving aircraft rotate and touch down in various segments 
of the airport and depart and arrive at different angles with respect to the ver- 
tical, (Departing aircraft rotate over the upwind two-thirds of the runway and 
climb out at 2 to 15 degrees, depending upon aircraft, load, atmosphere density, 
etc. Arriving aircraft usually aim for the 1000-foot marker and a _+ \ degree 
glide slope but may land over the first third of the runway.) The proposed sys- 
tem could have feedback from the controller concerning the following: 


• runways closed to heavy traffic (this would likely reduce 
vortex element drift to adjacent corridors), 

• whether aircraft are arriving or departing (this would vary 
the calculated initial altitude of the vortex elements), 

• the weight of the generating aircraft; etc. 


The refinement of this controller input into the computer would come with use 
of a prototype system. 
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• Simplified Flow Chart 

Figure 2-10 is a simplified flow chart of the computer operation for pre- 
dicting the turbulence conditions. The initial inputs into the computer are the 
vortex element coordinates and danger area coordinates. The element pairs 
are spaced longitudinally approximately 1000 feet apart along the runways and 
with increasing separation along the approach and departure corridors. The 
total number of elements is approximately 80 in the example or 20 per runway. 
To extend the elements to the outer marker, either more elements or wider 
spacing could be used. The danger areas are entered as coefficients for the 
line equations which form boundaries of the danger zones. For runway 9L in 
the example, the danger zone would be defined in two dimensions as the following 

x > 13.5, and 

10.8 -0.033x < y < 10.40 +0.033x 


As the model is extended, this two-dimensional representation of danger zones 
would be extended into three dimensions. 

Meteorological data, next input into the computer model, will include 
whatever real time data are available as well as upper air windfield data when 
they are available. The upper air windfield data can be supplied by the Aviation 
Forecaster who supplies local forecasts to the field. The local forecaster in 
turn will receive his data from available instrumentation (radiosonde, etc) plus 
data supplied by pilots using the terminal. (The air data systems of the wide- 
bodied aircraft (i.e., inertial platform, airspeed and temperature sensors) 
could supply valuable, near real-time data to the forecaster.) 

The wind field is next calculated. For the purpose of predicting the move- 
ment of the trailing vortices, the assumption will be made that the wind field in 
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Fig. 2-10 - Flow Chart of Computer Software for Wake Turbulence 

Predictive System 
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the vicinity of an airport can be represented by a deterministic steady- state 
mean wind. The mean wind will be assumed to be a function of height only (i.e., 
variations in the horizontal plane and over time intervals less than one hour are 
assumed to be negligible). First, the wind near the ground will be averaged over 
a 10- to 20-minute period. Then the winds at altitude will be provided by the 
forecaster as indicated in the previous paragraph. 

Wind speed and direction from near the surface to several hundred feet 
altitude will be fed into the computer on a periodic (likely hourly) basis. The 
computer will then generate least- square s polynomials (probably second- or 
third -degree) to fit the wind component data. The polynomials can then be 
used in the vortex transport prediction model to estimate the steady-state 
wind vector at any desired altitude. 

There are, of course, many approaches which may be taken in estimat- 
ing the wind profile in the lowest several hundred feet of the atmosphere. The 
approach taken here, however, eliminates the need for any elaborate and ex- 
pensive wind measuring systems such as instrumented towers, frequent re- 
leases of pilot balloons, and so on. An experienced forecaster, using the 
information available to him, should be able to produce reasonably accurate 
estimates of the winds up to 2000 feet — even when conditions are changing 
rapidly as with a frontal passage. The prediction model will be set up to handle 
more precise wind inputs, if they become available through improved instru- 
mentation or extrapolation techniques. 

After the wind field is calculated, each vortex element is individually 
operated upon by the computer for a sequence of time increments equal to its 
calculated lifetime. These operations are represented in abbreviated form 
in the remainder of Fig. 2-10. The computer determines if the element is in 
ground effect; if it is not in ground effect, the computer treats the vortex pair 
segment as a single element. If it is in ground effect, the computer handles 
the individual (port and starboard) elements independently. The computer 
next determines the element positions at each time increment and checks to 
see if the element is in the danger zone. The process is repeated for each 
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time increment and each vortex element until all vortex elements are observed 
throughout their active lifetimes. 

In particular, two areas deserve further study: 

1. The relationship between observed decay rates and the 
spectrum of the winds near the level of the vortex; and 

2. The relationship between decay and the Richardson 
number which is a nondimensional number containing 
the vertical shear of the mean wind and thermal sta- 
bility, both of which are important in determining 
atmospheric turbulence. 

• Confidence (Safety) Time Interval 

Since many of the variables used in the computations are subject to vary- 
ing degrees of accuracy, a safety factor will be built into the model to avoid 
prediction of safe conditions when hazardous conditions actually exist. This 
factor, in the form of a confidence interval in time, will be determined from 
estimated inaccuracies in aircraft trajectories, meteorological parameters, 
etc. (Note: The confidence time interval will be less for the hybrid system 
of Section 3.3 than for the predictive system due to the added input data of the 
hybrid system.) 

• Use of the System 

From the data of the persistence of vortex elements in the individual 
corridors, the Air Traffic Controller can mentally or from tables determine 
the spacing that is required between various classes of aircraft for a safe 
landing or departure. The controller would see how long is required for all 
vortex elements to clear each of the aircraft corridors for a specified gen- 
erating aircraft type, departure, arrival, departure angle, etc. He would 
then know the conditions of wake turbulence persistence and judge the sep- 
aration requirements following an aircraft of the indicated type. Use of the 
system is further discussed in Section 2.5. Communications and Action 
Command Processes. 
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2.4 INFORMATION ACQUISITION, STORAGE AND PROCESSING 

The information acquisition, storage and processing subsystem 
should be compatible with the field meteorological instrumentation, system 
software design and display/command process, and yet will be as in- 
expensive and simple to interface with existing air terminal facilities 
as is possible. 

• Information Acquisition 

The information acquisition subsystem should gather the data required 
to forecast the wake turbulence and enter these data into the wake turbulence 
predictive computer (Section 2.3). The information to be entered into the 
computer may be of the following forms: 

• Meteorological data from sensors presently installed at the 
airport 

• Meteorological data from sensors to be installed at the airport 

• Meteorological data from observers 

• Aircraft operations data (aircraft type, departure or landing, 
departure angle, etc.)* 

Meteorological data from sensors presently installed at the facility must 
be formatted to be compatible with the wake turbulence predictive system 
computer. This formatting may require analog -to -digital converters, digital- 
to-digital converters, multiplexers, etc. However, for an interim system, 
input of these data into the computer should probably be via computer 
console and operator to read the indicated values and enter the data* Sensors 
presently available at air terminals are located near or on the control tower 
and at remote locations around the terminal. Consideration must be given 
to cabling lengths, cable routes, remote analog-to-digital converters, reli- 
ability, maintenance, etc., if present instrumentation is used with an opera- 
tional prediction system. 


2-29 


LOCKHEED -HUNTSVILLE RESEARCH & ENGINEERING CENTER 



LMSC-HREC DZ259#* / 

. . . • ; -*.-*• * 

Economic factors should be considered in deciding whether to' use existr- 
ing meteorological sensors and converting their outputs to be compatible with 
the data processing subsystem or purchasing new sensors to be compatible'. 

In addition, new sensors might be required to provide additional meterological 
data to these data already available at air terminals. 

In addition to sensor -provided data, the system should be capable of re- 
ceiving data through the forecaster for the local meteorological conditions. 

An experienced forecaster will filter information from pilots, radiosondes, 
etc., to provide the required input into the model. Wind shear data from the 
air data systems of the 747s, LlOlls, and DC-lOs would constitute a valuable 
input into the system which could come through the local forecaster. The 
trained forecaster should be a valuable system link for providing forecasts 
of wind profiles to several hundred or even several thousand feet. 

It is envisioned that in an operational system, the air traffic controller 
will have a limited capability for communicating with the computer. He might 
supply the computer (Section 2.3) with the following data from a single console 
with a half-a-dozen push buttons: 

• Runways closed to heavy aircraft 

• Type aircraft generating vortices 

• Landing configuration 

• Departing configuration 

• Departing angle with respect to horizon 

• Information Storage and Processing 

Information storage and processing will be performed via digital com- 
puter for the wake turbulence predictive system. After the computer memory 
size and speed are defined, existing and planned computer facilities at the air 
terminals should be surveyed to determine if it is practical to add this burden 
to existing facilities or to determine if a "minicomputer’ 1 should be used for 
the predictive system which is independent of existing facilities. 
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The computer operation will be divided into three operations: 

• Collect meteorological and other pertinent data 

• Calculate wake turbulence conditions 

• Provide data for display in proper format. 

The computer size can be determined from summing the following storage 
requirements: 


(A) Pertinent Input Data Storage 

(B) Program Format 

(C) Arithmetic Functions 

(D) Display Data 

(E) Housekeeping 


Initial estimates indicate 8 to 16K of memory should be sufficient for storage 
and performing the desired calculations. 

Once computer storage, word length and computation time are speci- 
fied from knowledge of the wake turbulence prediction system software de- 
sign, a survey of existing and planned airport computer facilities would 
indicate the feasibility of integrating these data processing operations into 
existing or proposed (ARTS II) air terminal computational facilities. If 
existing or planned facilities are not compatible with requirements, a "mini- 
computer" will most probably perform the computational requirement at a 
cost of $25,000 to $50,000 per unit. 

2.5 COMMUNICATIONS AND ACTION COMMAND PROCESSES 

The communications and action command process for the wake pre- 
diction system should be compatible with the present National Aerospace 
System (NAS). The primary difference between the communication/ command - 
action systems for detective and predictive systems is that detective systems 
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may provide short time-period data for possible aircraft wave-offs while the 
outputs of a predictive system should provide sufficient warnings of hazardous 
conditions to allow spacing procedures to provide a mechanism of hazard 
avoidance. The following paragraphs will discuss several methods of com- 
munications and action command processes which have been or should be 
considered for integration into the NAS. These methods include the following: 

• Cockpit Display: Pilot takes precautionary action 

• Runway Display: Pilot takes precautionary action 

• Controller Display: ATC relays turbulence data; pilot takes 
precautionary action 

• Controller Display: ATC regulates traffic spacing according 
to turbulence forecast 


• Cockpit Display with Pilot Reaction 

M A cockpit display would remove the overburdened controllers from 
the loop” is a view advanced initially by some air traffic controllers when 
discussing the subject of wake turbulence advoidance. However, ATC is not 
really removed from the loop in landing operations because any turbulence 
avoidance action taken by a pilot in the terminal area affects other aircraft 
in the area; therefore, ATC is automatically involved. In addition the pilot, 
who is very busy when landing an aircraft, must decide what kind of action to 
take. For a small aircraft he could possibly land over the turbulence; but 
for larger aircraft, unless the runway is very long, he can only abort the 
landing thus affecting the whole traffic pattern. In addition, pilot prejudice 
would enter into the picture with some pilots aborting landings and others 
continuing their landings under the same turbulence conditions. From the 
equipment standpoint, a turbulence indicator in a larger aircraft equipped 
with sophisticated landing equipment might be a minor additional expense; 
however, for a smaller aircraft which is affected more by turbulence than 
the larger aircraft, the additional expense of the warning device would likely 
pose an economic burden upon the owner of the craft. 
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The runway lighting display eliminates the purchase of equipment burden 
upon the individual aircraft owner. A system similar to VASI (Visual Approach 
Slope Indicator) has been recommended as a possible turbulence warning mech- 
anism. Different colors could represent various degrees of turbulences. For 
departure this might work quite well. However, for approaches this system 
suffers from the same deficiencies as the cockpit display in that the pilot has 
to make a sudden decision at a time when he is already overburdened, and his 
decision likely affects the remainder of the airport traffic. ATC becomes 
involved if hazard avoidance maneuvers are necessary. 

• Display for Air Traffic Controller 

Today air traffic controllers are required to maintain a five -mile 
separation behind a "heavy" (>300,000 pounds) jet when the heavy aircraft is 
followed by a lighter aircraft. The controllers are also required to provide 
a turbulence warning if, in their judgment, wake turbulence constitutes a 
hazard. A typical tower message for landing behind a heavy jet on the same 
runway is "Number two to land, following Lockheed C-5A on final. Caution 
wake turbulence." For landing behind a heavy jet on a parallel runway closer 
than 2500 feet the message is similar, "Cleared to land runway 9-R — Caution 
wake turbulence —747 on final 9-L." The courts have declared the tower re- 
sponsible for providing judgment decisions of when it is safe to depart or 
land under turbulence conditions. Lockheed -Huntsville personnel, after 

consulting with numerous pilots and controllers have concluded that a system 
which would provide meaningful turbulence prediction data to the controllers 
would be most beneficial in helping the controllers to make more accurate 
decisions about hazardous conditions. These data would allow the controllers 


FAA Wake Turbulence Advisory Circular, AL No. 90-23B, 19 February 1971. 

"How the Courts Look at Wake Turbulence," R. H. Jones, FAA Turbulence 
Symposium, Washington, D.C., 22-24 March 1971. 
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to provide adequate spacing between the aircraft — both arriving and departing — 
to minimize turbulence hazard risks and maximum utilization of the air term- 
inal. A number of pilots and controllers alike have agreed that adequate but 
not excessive spacing, even without a turbulence message to the pilot, would 
be the best method of controlling wake turbulence hazard avoidance for the 
immediate future. 

For a wake turbulence predictive system using state-of-the-art hardware 
and vortex theory, predictions of conditions conducive to vortex dissipation 
should be valid for several minutes at least. The data into the model will, for 
most cases (i.e., except in frontal conditions), not vary excessively over a period 
of several minutes. Therefore, a forecast of the conditions should be valid for 
the same period of time. Control of aircraft final landing spacing begins 
several minutes prior to aircraft touchdown. The controllers, given a vortex 
condition forecast, should be able to vary spacing intervals according to the 
hazard condition. No sudden wave -off should occur unless some unforeseen 
phenomenon occurs. In this case, the wave-off would be handled in the con- 
ventional manner. 

A more advanced prediction system might provide the controller with 
additional data. Instead of forecasting meteorological conditions conducive 
to turbulence dissipation from meteorological data, the more advanced system 
might also take inputs of aircraft type, velocity and location and predict the 
hazard from and to other aircraft in the system. This type system would 
semi-automate the controllers judgment decision. 
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Section 3 

DETECTIVE WAKE VORTEX MONITORING SYSTEM 
3.1 INTRODUCTION AND SUMMARY 

The detective wake vortex monitoring system is very likely the ultimate 
wake vortex warning system. The detective system will aim at providing a long 
range (10 km from runway threshold) warning of actual measured vortex condi- 
tions in a real time manner. To be most effective, the status of the detected 
vortex is to be displayed simultaneously in an easily interpreted format to both 
the pilot and to the departure -arrival and local air traffic controllers. The air 
traffic controllers are to maintain aircraft spacing according to vortex move- 
ment and decay trends and the pilot, in coordination with the controllers, would 
make real time decisions from his cockpit monitor if action were required to 
avoid hazardous contact with vortices which were determined to be in his pro- 
jected flight path. While the predictive or hybrid system might increase maxi- 
mum traffic handling rates in say 50 to 70% of the operational conditions, the 
detective system should increase maximum traffic handling capacities in a 
higher percentage (say 80 to 90%) of the operational conditions because of the 
increased knowledge of vortex conditions and less margins for error in pre- 
dicting the hazard. 

Although the detective wake vortex monitoring system is the most 
effective type of vortex monitoring system from the standpoint of increasing 
air terminal capacities, it is the most expensive system to both develop and 
install in the field. In order to build a system which will provide effective 
real-time, long-range vortex data, a research program must be planned to 
develop more basic knowledge about the vortex characteristics and presently 
available hardware (lasers, detectors, etc.) must be "Mil-spec'ed 11 for long 
life operation with minimum maintenance. The state-of-the-art of data trans- 
mission from the controller to the pilot should evolve from the present verbal 
communication system to a digital-data link system or equivalent with a simply 
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formatted (possibly CRT) cockpit display of data communicated from the ground. 
Considerable experimentation should be performed in probing vortices at long 
ranges with laser Doppler instrumentation. In addition to development costs, 
the detective system — although probably much more effective than the predic- 
tive or hybrid systems — will cost more than the other systems to install be- 
cause of the additional hardware components. 

A considerable amount of this contractual effort (NAS8-26668) has been 
expended in the investigation of detective vortex monitoring systems. This 
chapter of the report will first briefly describe a typical detective system (3.2) 
and then discuss studies which led to the selection of this typical system. Sec- 
tion 3.3 discusses the advantages and disadvantages of one -component (one- 
dimensional) and two-component (two-dimensional) detective systems. The 
program used to develop these data is discussed in Section 3.4. Conclusions 
and recommendations for additional effort toward predevelopment of the detec- 
tive system are covered in Section 3.5. 

3.2 TYPICAL WAKE YORTEX DETECTIVE SYSTEM 

Figure 3-1 illustrates a typical system for monitoring wake vortices with 
a detective laser Doppler sensor. Figure 3- 1 , a two-component laser Doppler 
sensor — actually two each of single -component Doppler sensors scanning 
synchronously — provide data which are computer-processed to output two- 
dimensional visualizations of the vortex flow fields. The outputs of these sen- 
sors are processed either by a mini-computer or by a module of the Automated 
Radar Terminal System (ARTS III) which also performs a decision process in 
determining which vortices may be hazardous to oncoming aircraft. These 
vortices are displayed to the departure/arrival and local air traffic controllers 
via the plan position indicator (PPI) displays. Data are also relayed to the air- 
craft via the proposed digital data link and displayed in the cockpit by the multi- 
function display (MFD). These various functions are discussed in the following 
paragraphs. The subsystem subsections of Section 3 present the reasoning 
which led to the choice of the depicted system. 
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Proposed Multifunction 



Fig. 3-1 - Typical Detective Wake Vortex Monitoring System 
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3.2.1 Sensor Subsystem 

The sensor for the long range wake turbulence monitoring system is 
the laser Doppler velocimeter (LDV). This subsection will first describe 
the basic principle of the LDV; next a typical system layout at an air terminal 
will be described; and finally a number of alternate optical configurations for 
the LDV will be described. 

The basic principle of the LDV is described below: 


"An electromagnetic wave having an angular frequency, uk, 

propagating in the direction k^( |k.| = Uh/C = Zir /X^, C and X . 

being the velocity of light and the incident wavelength, re- 
spectively) and scattered^ by a particle moving with a velocity, 

V, into a direction k (|k 1 = Zn /X.) experiences an angular 

sc SC* 1 

frequency shift Ad which for nonrelativistic velocities is 
given by 

Ad = (It -k.)-V. 
sc i 


The detected frequency shift is due, therefore, to the velocity 
component of V in the direction (k — k.). Confining considera- 

S C X 

tions to a backscatter system, the Doppler shift Af ( = Au)/2 7 r ) 
is given by 


Af =2|V|cos0A, 


where V is the target velocity, X the wavelength, and 0 the 
angle subtended by the wind direction and the optical system 
line of sight" {Ref. 19). 


Figure 3-2a depicts atypical (Atlanta) airport layout and the arrange- 
ment of laser Doppler velocimeter (LDV) sensors to provide measured 
vortex status data for the air corridors entering and departing the air termi- 
nal. In this illustration, one of the parallel runways is the primary approach 
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Scan Region 



Fig. 3 -2a - Typical Airport Layout of Detective, Long-Range, Laser Doppler 
Wake Monitoring System 



Fig. 3 -2b - Typical Layout and Range of Two Component Laser Doppler 
Velocimeter for Specific Runway 
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runway and the other is the primary departure runway. LDV sensors are 
oriented to scan the appropriate corridors. Vortices are to be tracked within 
these corridors; therefore, flight within the corridors would not be considered 
safe in the presence of the vortex hazard. Aircraft entering or departing the 
terminal outside these corridors would do so at their own risk. The LDV 
sensors would be capable of scanning either end of the runway or both ends 
alternately depending upon the requirement. 

A typical range and dimension of an LDV sensor system is illustrated 
in Fig. 3-2b. The sensors are typically separated by approximately one kilo- 
meter . Separation requirements and range requirements are discussed in 
Sections 1.4 and 3.3 respectively. The requirement for two sensors to provide 
a two-dimensional picture of the vortex and the minimum separation between 
the sensors are discussed in Section 3.2. Basically, the two sensors are re- 
quired to obtain a clear picture of the vortex flow field under all ambient wind 
conditions and maximum practical angular separation is required between the 
sensors and the data point being examined. Sensor separation is determined by 
airport real estate limitations, in the figure to approximately 1 km separation 
between the sensors and 6 deg angular separation between the sensors at the 
most distant data point (10 km). 

• Typical Detective System Optical Configuration 

Selection of a specific optical configuration for a detective wake turbu- 
lence monitoring system sensor is probably premature at this point; however, 
several options are presented which might fulfill system sensor requirements. 
Basic optical configuration options are the following: 

• Coaxial (monostatic), focused, continuous wave 

• Coaxial (monostatic), pulsed 

• Coaxial (monostatic), frequency modulated, continuous wave 

• Bistatic focused, continuous wave 

• Bi static pulsed 
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Of the five configurations mentioned on the previous page, the coaxial 
focused laser Doppler velocimeter depicted schematically below is the optical 
configuration which has been most used to date. 


Plane Distant L from Telescope Primary 



1 


Recording 

Electronic 


"Radiation from the laser is focused at the required range via a New- 
tonian or Cassegrain telescope. Target range variation is provided by 
varying the position of the focusing lens which also serves to expand the 
small diamter laser beam to fill the primary telescope objective. The 
laser beam is divided by a beam splitter before it enters the telescope. 

A small portion of the laser energy is split off toward the detector. The 
Doppler shifted radiation scattered from the focal region of the telescope 
is collected by the telescope and allowed to reenter the laser cavity where 
it is amplified. Upon reemerging from the laser, a fraction of the fre- 
quency shifted component takes the path to the detector where is is photo - 
mixed with the original laser oscillation. The output of the detector is 
thus modulated at the frequency of the Doppler shift introduced by the 
motion of the particulate matter in the atmosphere at the telescope focus. 
For a 10.6 jl wavelength (carbon dioxide laser wavelength) the Doppler 
shift Af is given by 

Af = 188.5 v kHz 

where v is the velocity component along the optic axis in m/sec (Ref. 19). 

"Spatial resolution from such a system is obtained by the requirement 
that wave fronts must match at the detector surface to be efficiently 
photomixed. Since the original laser radiation emerges from the laser 
as a plane wave front, only the scattered radiation collected by the tele- 
scope that emerges from it as a collimated beam is heterodyned effi- 
ciently. Essentially this means that only radiation scattered from the 
region near the focus of the telescope contributes to the Doppler signal" 
(Ref. 19). 
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Figure 3-3 indicates the spatial resolution of this optical configuration 
as a function of range. This figure indicates that the spatial resolution of the 
coaxial configuration is probably too poor at long ranges for this system to be 
considered as a candidate for the long range wake monitor. 


The signal-to-noise ratio (SNR) for the coaxial system in the near field 
appears (Ref. 20) as: 


SNR 


1 naAN 
~ 4 V B 


where 


rj = efficiency of the optical system 
n a = scattering cross section of the medium per unit volume 
X = laser wavelength 

N = number of photons per second emitted by laser 
B = electronic bandwidth 

This equation appears independent of range and optic diameter. However, as 
the range to the focal volume increases with a given coaxial optical configura- 
tion, the focal volume size increases and the velocity distribution within the 
focal volume spreads, requiring the bandwidth of the electronics to be increased. 
Thus, from the above equation the resultant SNR is decreased proportionately 
to the increased bandwidth requirement. 

Signal-to-noise is affected by scan rate since scan rate and range deter- 
mine time -on-target which in turn contributes to the signal bandwidth. This is 
depicted in Fig. 3-4. 

The coaxial pulsed system depends upon pulse length for spatial resolu- 
tion. However, as pulse lengths are shortened to improve spatial resolution, 
velocity resolution decreases proportionately (see Fig. 3-5). Interpulse 
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Scan Rate, rad/sec 

Fig. 3-4 - Signal-to-Noise Ratio for Coherent, Continuous -Wave 
Coaxial Configuration as a Function of Scan Rate 
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modulation is used in microwave radar systems (Ref. 21) to improve spatial 
resolution. However, interpulse modulation with a CO^ laser system would 
add sophistication to the overall system design and likely increase costs and 
reduce reliability. 

Signal -to -noise per pulse for the coaxial pulsed configuration can be 
given by the same equation as for the standard coaxial system when the data 
point of interest lies in near field (R < D /A) with N redefined as the number 
of photons per pulse. However, when the data point lies in the far field 
[R > D 2 /A), the following equation holds for SNR (Ref. 22) 

SNR = r)[J /hv)no c TD 2 /r 2 

where 


A, rj and no are, as previously defined 
J = net output energy in observation time 
c = speed of light 
T = pulse duration 
hv = laser photon energy 
D = optic operative diameter 
R = range to data point 

In both near and far field the SNR is proportional to the average energy illumi- 
nating the volume of interest. To obtain sufficient average energy to provide 
adequate SNR at long ranges with present technology requires sophisticated 
laser amplifiers or a high power laser operating in a Q-switched or 
mode-locked configuration. The power efficiencies of these types of laser 
configurations are nominally much lower than for continuous wave configura- 
tions. Since S/N depends upon the energy illuminating the target, S/N is also 
dependent upon scan rate as depicted in Fig. 3-6. 

The frequency modulated coaxial system depends upon linearly varying 
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Fig. 3-6 - Signal-to-Noise Ratio for Coherent, Pulsed, Coaxial Configuration as a Function of Scan Rate 
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the frequency of the iaser transmitter with time so as to oscillate above and 
below the mean frequency. The distance to the target is determined by com- 
paring the phases of the received signal and local oscillator. The Doppler 
shift also translates the phase of the returned signal with respect to the trans- 
mitted signal; so that the difference in frequency in the two halves of the 
modulation period is a measure of the Doppler shift, which can be incorporated 
into the range calculation to determine Doppler shift (velocity) as well as 
range. Range error for this type system can be shown to be c-f Af, where c 
is the speed of light, and Af is the frequency excursion of the modulation. For 
Af = 5 MHz, the range error is 50 feet. This type system is again more so- 
phisticated than the coaxial system and is in its infancy with regard to appli- 
cation at infrared frequencies. 

The bistatic focused laser Doppler velocimeter, being developed for NASA 
by Lockheed (Ref. 23) offers hope: of providing good spatial resolution at long 
ranges (Fig. 3-3) without sophisticated pulsing or modulation techniques. The 
bistatic laser Doppler velocimeter consists of a separate transmitting and 
receiving set of optics designed to intersect the transmitted beam’s focal 
volume and receiver telescope's field of view (Fig. 3-7). The primary so- 
phistication of the bistatic system involves the me chani cal/ optical difficulties 
of maintaining this intersection. If the bistatic system is constructed with a 
fixed mirror surface as indicated in Fig. 3-7, this mechanical alignment 
problem is lessened. The signal -to-noise ratio for the scanning bistatic 
system is depicted in Fig. 3-8. The scanning of the rather large mechanism 
presents interesting mechanical problems which, along with other aspects 
of the bistatic system, are discussed in detail in the recent Lockheed report 
of Ref. 24. 

Bistatic pulsed systems offer the resolution of the bistatic configuration 
at short ranges and the resolution of the pulsed configuration at the longer 
ranges. This type system would require the combined complexities of the 
bistatic system and the pulsed system, which does not appear practical at 
the present time. 
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Fig. 3-7 - Fixed -Surface Bistatic Configuration 
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3.2.3 Computer Subsystem 

The computer subsystem for the typical detective wake monitoring 
system will consist of a mini-computer interfaced with the sensor pre- 
processor outputs and the automated air traffic control system. The mini- 
computer performs two functions: 

• Control of the scanning process 

• Processing of sensor output data 

In performing the control function, the mini-computer forces the indi- 
vidual sensors to scan a preprogrammed pattern. A coarse scan is performed 
first at predetermined ranges to locate any irregularities in the normal wind 
patterns. Once an anomaly is discovered, the computer forces the scanners into 
a detailed scan routine to further determine if the anomaly is a vortex, and if 
it is a vortex, determine the vortex characteristics. Once a vortex is located, 
the sensors will likely be controlled to scan rapidly along the vortex to deter- 
mine its position within the air corridor. 

The data processing function of the mini-computer receives the pre- 
processed outputs of the individual (one -dimensional) sensors and records 
them in storage. After sufficient data are in storage, the computer combines 
the one -dimensional data to form a two-dimensional picture of the coarse scan 
area. The computer looks for anomalies (mentioned in the above paragraph) 
to determine if a vortex exists in the corridor. The outputs from the detail 
scan are processed to provide a two-dimensional picture of the vortex and to 
calculate circulation, position, movement and other pertinent parameters. 

The mini -computer of the detective system, in conjunction with the local 
automated air traffic control system, integrates vortex data with traffic data 
to determine if a hazard exists to aircraft in the area. The computer per- 
forms this operation with stored vortex data and inputs of projected aircraft 
flight profiles. If the computer determines that a vortex of hazardous strength 
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lies in the projected path of an oncoming aircraft, a warning is displayed 
via the display subsystem to the controllers and pilot involved. 

3.2.4 Data Display/Communi cation Subsystem 

The method of performing a real time display of wake vortex hazard 
data to both the pilots and controllers is most important in the development 
of a detective wake monitoring system. The data display must provide data 
to the persons involved in a simple yet effective format. If the pilot is to 
make the decision to abort a landing or maneuver to avoid a vortex encounter, 
he must have sufficient data to know the most appropriate move to make. 
Equally important, the pilot should not be required to make an evasive ma- 
neuver if the vortex to be encountered is not of sufficient strength to constitute 
a hazard. To be considered along with these factors, aircraft crew members 
are busiest during approach and departure portions of flight and must not be 
over burdened with complicated display formats of vortex conditions. 

An effective real time display can perhaps be incorporated into the pro- 
posed Multifunction Display (MFD) (Ref. 25), a cathode-ray tube monitor proposed 
to display data provided from the ground via a digital data link. The digital 
data link will provide computer -to -computer communication between the 
ground-based traffic control computers and airbourne flight computers. The 
display could perhaps provide an indication of vortex hazards along the pro- 
jected flight path together with other pertinent navigation data already planned 
for the MFD. 

3.3 CONSIDERATIONS OF MULTI-COMPONENT (MULTI-DIMENSIONAL) 
AND SINGLE -COMPONENT (ONE -DIMENSIONAL) SYSTEMS 

Should a detective wake vortex monitoring system be a one-, two-, or 
three -component detection system? This question has arisen since a laser 
Doppler detective wake monitoring system was first conceived. An attempt 
is made here to clarify the advantages and disadvantages of the single and 
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multicomponent systems and to indicate the preferred choice of system with 
today's knowledge of the vortex monitoring problem. 

3.3.1 Single -Component (One -Dimensional) Vortex Sensor 

A number of computer simulations have been made of a laser Doppler 
type sensor scanning through a vortex. The simulation consists of looking at 
the vortex flowfield with a course resolution instrument and attempting to 
recognize profiles from the resulting output. The simulated output of the 
sensor is plotted by an SC 4020 plotter as velocity versus locations within the 
vortex. Figure 3-9 defines the terms which were used in plots of the simu- 
lated sensor output. This simulation program is discussed in Section 3.4 and 
in Appendixes B and C. 

A single -component sensor of vortex velocities can provide a realistic 
picture of the vortex tangential flow field of a vortex if it views the vortex 
at an angle perpendicular or nearly perpendicular to the axis of the vortex. 
Figure 3-10 illustrates this fact. The W (sensed velocity) and T, L (actual 
velocities) curves agree very well in this figure to indicate a good measure- 
ment of tangetial velocities with a one -component sensor viewing the vortex 
at angles approaching 90 deg with respect to the vortex axis. However, as 
the viewing angle with respect to the vortex axis (angle between "site line" 
and "vortex axis") approaches zero degrees, the ability of a single component 
sensor to sense accurately the vortex tangetial velocity degrades rapidly. 

This results from the fact that a coaxial LDV measures the along -line - 
of-site velocity component. Figure 3-11 shows a comparison between a sim- 
ulated one-dimensional sensor* s output and the actual tangential flow within 
the vortex as the angle between the line -of- sight of the sensor and the vortex 
approach zero degrees. In the upper set of curves of Fig. 3-11, this angle is 
approximately 2 deg, in the middle set, the angle is approximately 6 deg and 
in the lower set the angle is approximately 10 deg. Figures 3-10 and 3-11 
show that as expected, the velocity resolution capability of a one -dimensional 
sensor viewing a vortex tangential flow field is greatly dependent upon the 
viewing angle, which improves as the angle (between sensor line-of-sight and 
vortex axis) approaches 90 deg. 
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Definitions for Simulations 


Vortex Axis; 
Line of Sight; 


Scan Plane; 
Site Line: 


runs along center of vortex core Reference Sweep Line: line in scan plane which 

passes through vortex 

line of sight of the sensor as it scans axis and is perpendicular 

across the vortex to site line 


formed by line of sights T 

line of sight which intersects 

vortex axis L 

W 

X 


actual tangential velocity of the vortex at the 
focal point of the simulated system 

actual line of sight velocity of the vortex at 
the focal point or range rate point of the 
simulated system 

Simulated velocity output of sensor resulting 
from logic described in Section 3.4 

Output from data processor which provides 
maximum absolute value of velocity within 
sensor focal volume or pulse length (X is 
not plotted for all figures). 


Fig. 3-9 - Laser Doppler Velocimeter Making Single Line Scans Through a Vortex 
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Fig. 3-10 - Viewing Vortex Tangential Flow Perpendicular to and at +60 Degrees 
to Vortex Axis (Range 300 m) *” 
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Viewing Vortex Tangential Flow at Angles of 1.7, 5.7 and Approxi- 
mately 10 Degrees (top to bottom curves, respectively) to vortex 
axis (Range 300 m) 


3-22 


LOCKHEED - HUNTSVILLE RESEARCH & ENGINEERING CENTER 







LMSC-HREC D225936 


For a short range system, positioning the sensor to view the vortex 
with a line -of -sight nearly perpendicular to the vortex axis is a relatively 
simple matter. The sensor could be placed under the flight path while look- 
ing vertically (Fig. 3-12a)or to the side of the flight path to look-over to the 
vortex (Fig.. 3- 12b). In either case the sensor would view the vortex at an 
angle approaching perpendicular to the vortex axis. 

However, for a long range sensor (the topic of Section 3), the location 
of the sensor such that the vortex could be viewed at a large angle (approxi- 
mately 90 deg between the sensor line -of-sight and the vortex axis) becomes 
a most impractical task. In order to obtain this approximate 90 deg viewing 
angle, the sensor would necessarily be located 5 to 10 kilometers from the 
runway edge and scan back toward the flight path. Because of the shortage 
of air terminal real estate and the normal physical obstacles — terrain, 
buildings, etc., —in the sensor line -of-sight, locating the sensors to view 
the vortex perpendicularly at long ranges will probably be impractical. The 
vortex will normally be viewed as depicted in Fig. 3 -12c with a small angle 
between the sensor line -of -sight and the vortex axis. The resultant measured 
velocities will be much smaller than actual velocities as indicated in Fig. 3-11, 

# Axial Flow Considerations for Single -Component Systems 

The above paragraphs indicate difficulties in characterizing a vortex 
with a long-range, single -dimension sensor if only the tangential velocities 
are considered. Figure 3-13 illustrates a one -dimensional sensor's output as 
it scans through a vortex with simulated axial (Newman's model, Appendix A), 
as well as tangetial flows. It is immediately evident that an axial flow com- 
ponent similar to the one modeled in Fig. 3-13 will greatly enhance the ability 
of a single component sensor to detect a vortex at small angles with respect 
to the vortex axis and at long ranges. However, a minimum of experimental 
evidence exists to verify the existence of an axial flow such as the one de- 
picted in Fig. 3-13. There is very little theory or data to correlate either 
the strength or dimension of the axial flow 'with the generating conditions; 
i.e., aircraft circulation, etc. Also, tangentail — not axial velocities — upset 
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Fig. 3- 12a - Sensor Located Under Flight Path Looking Up (Viewing 
vortex perpendicular to vortex axis) 



Fig. 3-12b - Sensor Located to Side of Flight Path Looking Toward 

Flight Path (Viewing vortex perpendicular to vortex axis) 



Fig. 3- 12c - Long Range System — Angle Formed by Sensor Line-of- 
Sight and Vortex Axis (a) Necessarily Small Due to Lack 
of Real Estate for Locating LDV far away from Flight Path 
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Fig. 3-13 - 


One-Component System with Combined Axial and Tangential Flows 
Viewed at Angles of 1.7, 5.7, and ~ 10 Degrees to Vortex Axis {jNTote: 
"T" represents total (axial plus. tangential) simulated velocities in- 
stead of tangential velocity as in; previous plotsj Maximum velocity 
was 38 m/sec; maximum tangential flow was 34/ sec 
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aircraft. Therefore, it is hardly worthwhile to begin immediately to design 
a wake turbulence detection system based upon the existence of vortex axial 
flow. 


• Conclusions about Utility of One -Component Long-Range Vortex Detection 

If the orientation of the vortex were known, it might be possible to derive 
a true, long-range (2 to 10 km) picture of a vortex from one -dimensional data 
taken at a small angle (such as illustrated in Fig. 3-11). The vortex flow could 
be derived by dividing the single component velocity data by the sine of the 
viewing angle (with respect to the vortex axis)F. The drawbacks to this type 
system are the following: 

• Difficulty in knowing exact orientation of the vortex. 

• Difficulty in distinguishing between sensed vortex velocity 
component and wind component. 

• Difficulty in resolving axial and tangetial velocity contributions 
to measured data. 

With a high spatial resolution sensor, it would be possible to determine 
the orientation of the vortex with some degree of accuracy. However, at small 
angles, a small error in angle determination would mean a large error in the 
resolved tangential velocity; (i.e., since sine 4 deg = 0.07 and sine 5 deg = 0.09; 
a 1 deg error in vortex angular orientation would result in an approximately 
25% error in calculated vortex tangential velocity). Since the spatial resolu- 
tion of a high resolution bistatic velocimeter (l/3 meter diameter optics, 2 
meter separation) is approximately 500 meters at 5 km range and 2000 meters 
at 10 km range (Fig. 3-3), it would be difficult to determine the vortex orienta- 
tion to within +1 degree at the larger ranges. 

Viewing a vortex at long ranges at a small angle with the axis results in 
sensing a small velocity component of the tangential flow. Figure 3-11 illustrates 
these sensed velocities being typically below 6 meters/sec for 747 vortices 
1000 meters behind the aircraft at 10 deg viewing angle and below 4 meters/ 
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sec at 6 deg viewing angles. A 5 to 10 meter/sec wind gusting directly down 
the line -of-sight of the LDV sensor would not be an uncommon atmospheric 
phenomenon at many airports. With a one -dimensional sensor, the difficulties 
of unscrambling a vortex velocity component less than 5 meters/sec on top of 
a turbulent wind component of 10 meters/sec are obvious. 

Since vortex axial velocities are presently ill defined, it is difficult 
to perform meaningful analysis of viewing vortices with a one -dimensional 
sensor from a direction such that the axial component is likely to be the pre- 
dominant velocity component. Figure 3-13 presents an attempt to simulate 
such a case with axial flow according to Newman's model (Ref. 26 and Appendix 
A). The presence of the axial flow might greatly facilitate the monitoring of the 
vortex at long ranges with a single component velocity sensor. At present the 
pressence of an ill defined axial flow only confuses the concept of a one- 
dimensional, long-range vortex sensor. 

3.3.2 Two Component (Two-Dimensional) Vortex Sensor 

A two-dimensional LDV concept would use two each of one -dimensional 
LDV sensors viewing the same data volume and vectorially computing a resultant 
vector velocity from the one -dimensional sensor outputs to provide two- 
dimensional data in the plane formed by the sensor line -of- sights. Since the 
two sensors do not necessarily give orthogonal coordinate components, the 
summing of their outputs is somewhat more involved than a normal vector 
summation. This can be illustrated graphically (Fig. 3 -14a) by constructing 
perpendiculars from the vector tips of the velocities sensed by the individual 
LDVs. The intersection of the two perpendiculars defines the vortex velocity 
in the plane defined by the two LDV sensors and the data point. 

Figures 3- 15a and b compare vortex tangential flow visualizations by 
two-component (Fig. 3-15a) and single -component (Fig. 3-1 5b) sensor systems. 

In both figures the range is 5 km and the angle between the system sensor(s) 
and the vortex line -of - sight is six degrees. In Fig. 3-15b the vortex is viewed 
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by a single sensor with a 6 deg angle. In Fig. 3- 15a each sensor views the 
vortex with a 6 deg angle with respect to the vortex line-of-sight; in this 
figure the sensors are located on opposite sides of the vortex line-of-sight. 

Note the improved agreement between measured (W) and actual (T) vortex 
velocities in Fig. 3-15a over Fig. 3-15b. 

Figures 3-l6a and 3-l6b illustrate the effect upon the two-dimensional 
system by increasing and decreasing the angular separation (angle j3 of Fig. 3- 
14a) of the sensors. The abscissa of these curves can be expressed as the 
distance from the vortex axis (Fig. 3 -14b) in meters. The ordinates are the 
resolved velocities from the mathematically combined outputs (Fig. 3 -14a) 
of the two sensors. The curves marked T are actual tangential velocities at 
the data points. The curves marked L are the actual line-of-sight velocities 
at the data points which have been resolved by combining as depicted in Fig. 3- 
14a. The W curves are the sensor output velocities combined mathematically 
as again depicted in 3 -14a. The X curves are the absolute values of the sensor 
outputs combined mathematically (positive and negative velocities are not dis- 
tinguished from one another). Figure 3-l6a indicates that the agreement 
between the W and L-T curves degrades as the range to the data point ( n Vortex 
Range") increases from 1 km to 5 km and again to 10 km. This deviation be- 
tween actual velocity and sensed velocity increases for the 300 meter separation 
between sensors (Fig. 3-l6b). Both of these figures illustrate that measurement 
errors are magnified as the angular separation subtended by the two sensors as 
the data point decreases. 

Figures 3-l6a and 3-l6b are the resolved outputs (theoretical) of the 
LDV sensors with capabilities of resolving the sensed velocities within plus 
or minus one meter per second. With a better velocity resolution for example, 
(+1/3 m/sec) for each of the LDVs, the agreement of the W and T-L curves 
would improve. However, it was felt that since natural atmospheric (Kolmo- 
gorov) turbulence would be added to the LDV outputs in a real system (and has 
not yet been added in this simulation) that + 1 m/sec is a good conservative 
velocity resolution to expect from a LDV sensor operating at long range. 
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Fig. 3- 15a - Two-Dimensional Visualization of Vortex Derived from Two 

One -Dimensional Sensors. Sensors are Located on Either Side 
of the Vortex Axis and View the Vortex at Angles of Six Degrees 
with Respect to the Vortex Axis (range is 5 km) 



Fig. 3-15b - One -Dimensional Visualization of Vortex Derived from a One- 
Dimensional Sensor (range is also 5 km) Viewing Vortex at 
Six Degree Angle with Respect to the Vortex Axis 
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Fig. 3-l6a - Two-Dimensional Data Resolved from Single Azimuth Scans Through 
a Vortex by Two Laser Doppler Velocimeters Separated by 1000 
Meters (no axial velocity) 
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Fig. 3-l6b - Two-Dimensional Data with Velocimeter Separation of 300 Meters 
(no axial velocity) 
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• Misalignment Effects on Two -Component Sensors 

Requirements to maintain independent optical pointing accuracies be- 
yond certain limits for scanning type systems greatly increase the cost and 
reduce the reliability of system hardware. Thus, investigation of point accu- 
racy requirements early in the conceptual phases of system development was 
an important aspect of the study. Beams from sensors located on separate 
platforms must intersect spatially to provide data from the same points in 
space (Fig. 3 -14a) for a true two-dimensional picture of the vortex. The 
beams need not pass through the same spatial point at the same instant of 
time. However, if they do pass through the points at different times, sensor 
outputs must be stored until both sensors have taken data from the spatial 
points of interest. The accuracy required for steering these beams is the 
subject of this subsection. 

Figures 3-17, 3-18, and 3-19 illustrate the effects of misaligning the 
sensor beams in the vertical plane by 1, 3, and 5 meters respectively at the 
target points. With a vertical misalignment, the sensors' beams no longer 
intersect at the target points; but one beam is skewed above the other. 

Figure 3-17 depicts a one meter misalignment. The "W" curves of this 
figure are practically indistinguishable from those of Fig. 3-15 (no misalign- 
ment). Thus, it can be concluded that a one-meter misalignment has little 
effect upon the system measurement accuracy. The vortex peak tangential 
velocity and core diameter are identifiable as in Fig, 3-15 where there is no 
alignment error. 

Figure 3-18 depicts a 3 -meter alignment error. In Fig. 3-18, peak 
tangential velocities have dropped approximately 10%, yet core diameters have 
remained approximately constant. The center of the core has shifted approxi- 
mately one meter. A three-meter misalignment is probably permissible for 
a long-range, two -component, wake turbulence measurement system. 
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DATA POINT POSITION WITH RESPECT TO SITE LINE IN METERS 
VORTEX RANGE IS 1 .000X10* 0? 40 



Fig. 3-17 - One-Meter Vertical Misalignment of LDVs for Ranges of 1, 5 and 
10 Kilometers (1 km separation between sensors) 
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Fig. 3-18 - Three-Meter Vertical Misalignment of LDVs 
for Ranges of 1, 5 and 10 Kilometers 
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Fig. 3-19 - Five-Meter Vertical Misalignment of LDVs 
for Ranges of 1, 5 and 10 Kilometers 
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A five-meter alignment error is depicted in Fig. 3-19. At the 10 km 
range, the peak tangential velocity sensed by the system has dropped approx- 
imately 40%, and the core diameter has increased from 6 meters to approx- 
imately 11 meters. At shorter ranges, the peak tangential velocity error is 
decreased somewhat. 

From these figures, it can be seen that a 3 meter pointing accuracy is 
adequate and possibly too stringent. This amounts to a one milliradian pointing 
error at a range of 3 km and a 0.3 milliradian error at a 10 km range. 

Runs were also made with a horizontal misalignment rather than the 
vertical misalignment discussed above. The horizontal misalignment merely 
intersects the two beams at a point slightly in front of or behind the nominal 
intersection rather than causing the beams to skew as discussed above. The 
horizontal misalignment did not have a noticeable effect for the 1, 3, and 5 
meter alignment error cases. The vertical misalignment represents a "worst” 
case misalignment. 

• Axial Velocity Effects on Two-Component Systems 

Figure 3 -20a depicts a two-component scan through a vortex (with sensor 
outputs combined as illustrated in Fig. 3 -14a with simultaneous axial and tan- 
gential flows. Figure 3 -20b depicts the same sensor scanning through a vortex 
with the same tangential flow but no axial flow. The similarities and differences 
in system output can be seen. These outputs are listed below for several data 
point positions: 


Distance from 
Site Line (m) 

Combined Axial and Tangential 

Tangential Only 

Actual 

Measured 

Actual 

Measured 

+ 20 

5 m/ s 

5 m/ s 

6 m/ s 

6 m/ s 

+ 10 

12 

15 

12 

15 

+ 5 

25 

25 

22 

22 

+ 2 

37 

37 

32 

32 

0 

37 

37 

0 

0 
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Fig. 3-20a - Two-Component Scan Through Vortex with Combined Axial and 
Tangential Flows 



Fig. 3-20b - Two-Component Scan Through Vortex with Tangential Flow Only 
Simulated 
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The outputs of the sensors for the vortices both with and without axial flows 
agree quite well until the data points move into the cores. Since axial flows 
are confined primarily to the core in Newman’s model (Ref. 26) and Appendix 
A) this would be expected. Once the sensors probe the core, the vortex with 
the axial flow displays velocities which continue to increase, while the vortex 
with tangential flow only displays decreasing velocities which approach zero 
at the core center. 

From the above data the two -component system can be seen to be capable 
of resolving tangential vortex velocities even in the presence of axial flow. 

This important property of the two-component system would allow vortex 
circulation to be calculated and thus vortex hazard somewhat independently 
of an axial flow presence within the vortex core. 

• Ambient wind effects upon component systems 

Figures 3-21 and 3-22 depict two-component scans through vortices with am- 
bient winds superimposed on their flows. Figure 3-21 simulates vortex tangen- 
tial flow without axial flow while Fig. 3-22 simulates combined axial and tangential 
vortex flows. These two figures are discussed in the following paragraphs. They 
represent a first cut at adding ambient wind effects to the simulations. 

Figure 3-21 depicts two -component scans through vortices at several 
ranges with an ambient wind of 5 meters/second superimposed upon the tangen- 
tial flow field of the vortex (no axial flow is included). The ambient wind in this 
case was directed along the vortex axis. This would correspond to a headwind 
in a real situation. A comparison of Fig. 3-21 with Fig. 3 -20b — the same con- 
figuration without ambient wind — indicates that the introduction of ambient wind 
increases the actual flowfield velocity (T) in the wings of the vortex (at radii 
greater than 20m with respect to site line) but has little effect upon core veloc- 
ities. Since the ambient wind simulated in this figure is normal to the vortex 
flow (tangential flow only simulated), the ambient wind vector (5 m/sec) — small 
when compared to the core max velocity (33 m/sec) — is oriented at 90 deg to the 
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Fig, 3-21 - Two-Component Scans Through Vortex Tangential Flows at 
Several Ranges with Ambient Headwind of 5 m/sec (no axial 
flow) 
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Fig. 3-22 - Two-Component Scans Through Vortex Tangential and Axial 
Flows at Several Ranges with Ambient Headwind of 5 m/sec 
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core velocity vector and has very little influence upon the combined resultant 
vector magnitude (34 m/sec). However, in the vortex wings, the ambient 
wind vector (5 m/sec) is comparable to the vortex flow vector (approximately 
4 m/sec) and the two combine to form a resultant of approximately 8 m/sec. 

The measured velocity (W) of Figure 3-21 closely agrees with the actual 
velocity (T) at close range (1 km) but falls 30% below the actual at the longer 
range (5 km). This can be attributed to the means of computing the two-component 
flow field from outputs of two one component sensors (Fig. 3 -14a). The outputs 
of each single component sensor are extrapolated by the computer to be the 
highest velocity magnitude within the focal volume. If the velocity is positive 
a positive sign is given to the magnitude and if it is negative, the velocity mag- 
nitude is termed negative. The positive and/or negative sensor outputs are then 
added according to Figure 3 -14a to develop a velocity magnitude for the resultant 
vector velocity in the plane formed by the line -of- sights of the two sensors. This 
mode of computing the resultant vector was not exact but provided a vector mag- 
nitude which compared very favorably with the simulated flow velocities for all 
vortices simulated until this particular case which included effects of ambient 
winds. 

As can be seen from Fig. 3-21, another means should be used to improve 
the agreement at long ranges (narrow view angles) between actual simulated 
velocities (T) and measured velocities (W) where winds are involved. One means 
of improving this agreement would be the subtraction of ambient wind velocity 
components from the outputs of each of the two sensors before combining the 
two sensor outputs to form the two component velocity magnitude. The ambient 
wind velocities could be measured beyond the influence region of the vortex 
phenomenon. 

Figure 3-22 depicts two -component scans through vortices at several 
ranges with an ambient wind of 5 meters/ second superimposed upon the combined 
tangential and axial flow fields of the vortex. The ambient wind is directed along 
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the vortex axis as a headwind as in Fig. 3-21. The figure corresponding to 
3-22 but without the ambient winds is 3 -20a. In Fig. 3-22 the ambient wind 
translates the flowfield velocity upward in the wings (20 meters outside the 
vortex axis) but does not have much effect within the vortex core. The meas- 
ured (W) velocities and actual (T) velocities again agree very well at the short 
range (approximately 1 km) but drop below the actual on the edges of the core 
(3 m to 10 m radii) at the longer ranges. However, since the ambient wind is 
in the same direction as the axial component, the measured (W) data agrees 
very well with the actual (T) at the vortex axis. If the ambient wind were a 
cross wind, this agreement may not be quite as good. Cross wind case should 
be simulated in the future to evalute the effectiveness of this type of system. 

3.3.3 Three (Three-Dimensional) or More Sensed Components 

The objective of the long-range detective system is to detect the vortex 
and determine the hazard created by the vortex flow field. Since the hazard 
from a vortex is considered to be due to its. tangential flow, one objective of a 
detective system might be to determine this tangential flow field; i.e., the vortex 
tangential flow versus distance from the vortex axis. The two-component sensor 
configuration defines the vortex flow in the plane formed by the lines of sight of 
the sensors as they probe the vortex. The output of the two-component configura- 
tion has been shown to be adequate in defining the vortex tangential flow field as 
long as scans are made through the vortex (i.e., the scan plane is moved through 
the vortex) in such a manner that the vortex tangential flow vectors lie in the 
scan plane at some time during the scan through the vortex. Since the two- 
component scan has been shown to be adequate to provide the desired data for 
the long-range detective system, the requirement for a more complex three 
(or more) component system is unjustified. 

3.4 SIMULATION PROGRAM: VORTEX FLOW FIELD AND FLOWFIELD 

MONITORING WITH LASER DOPPLER SENSOR 

Early in the contract the decision was made to emphasize mathematical 
modeling of laser Doppler sensors to detect simulated vortex flow fields to 
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develop further knowledge of sensor requirements and capabilities for a wake 
turbulence monitoring system. This modeling effort involved first choosing 
a simple but realistic vortex flow field model and, second modeling a simu- 
lated laser Doppler sensor scanning through the flow field. The vortex flow 
model would provide tangential flows since these were well defined and option- 
ally axial flows (not well defined) according to Newman 1 s model (Ref, 26). The 
flowfield simulation will have the capability of superimposing ambient winds. 
The sensor simulation is to consider optionally coaxial and bistatic optical 
configurations operating in both the pulsed and continuous wave laser con- 
figurations. 

3.4.1 Vortex Flow Field/Sensor-Detection Model 


The vortex tangential flow field (V^) for the sensor simulation is repre- 
sented by the equation (Ref. 26): 


where 


and 


£a_ 

27TD 


1 - exp 


y 4pz ) 


r = circulation 
o 

= freestream velocity 
V = eddy viscosity 
z = vortex axial coordinate 
D = distance from vortex centerline 


r n = 


4W 


o npVb 


W = V 
co 


V = k [l/p 
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W is total aircraft weight 
p is density of atmosphere 
V is aircraft velocity 
b is wing span 
t± is coefficient of viscosity 

k is the coefficient that relates laminar kinematic 
viscosity to eddy viscosity for a particular 7^ 

The vortex axial flow field (V^) is provided optionally as the following 
(Ref. 26): 

D / W D 2 \ 

v a = i^Tz-p \-frr) 

where D q is profile drag of the aerofoil and the other variables are those 
defined above. 

In order to relate the tangential (V^,) and axial (V^) flow fields of the 
vortex to some velocity sensor for a simulation of the sensor/flowfield inter- 
action, a coordinate system was set up to describe the sensor to the flowfield. 
The coordinate system chosen was spherical polar with the sensor at the origin 
(O) and the vortex axis direction defined by the vector PQ. 


Z 
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The velocity at point S is determined by D of the above equations which 
is the distance of S from the vortex centerline (PQ). The distance from S to 
PQ is: 

lsp X pq! 

|PQ| 

where 

SP = OP - OS 


The tangential velocity vector at S is 


VT 


V, 


(SP x PQ) 
SP x pqI 


and the magnitude of the component of VT in the direction of SO (the sensor) 
is 


V 

V ST 


VT * OS 
OS 


The magnitude of the axial component Vg ^ sensed by the sensor can be com- 
puted in a similar manner. Ambient winds can also be added vectorially to 
the vortex flow field. 

3.4.2 Laser Doppler Sensor Model 

The laser Doppler sensor model allows the simulated sensor to view the 
flow field as a typical laser Doppler sensor would, and it provides an output of 
velocity in terms of range and angular position of the sensor. The type of out- 
put depends upon the type of sensor (continuous wave or pulsed) and the type of 
processor (spectrum analyzer, filter bank, etc.) modeled. 
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• Coaxial Focused, Continuous Wave Laser Sensor 


The focused sensor does not make a point measurement, but samples a 
finite volume of space along a line (the focal volume, Fig. 3-23a) and con- 
sequently observes a variety of velocities (Fig. 3-23b). The signal is the 
strongest at the focal point (Fig. 3-23c) and tapers off as the distance from 
the focal point increases. The functional dependence of the signal on range 
L is given by 


I 



K 

(L - f) 2 ~ 
AL 2 J 


(Ref. 20) 


where 

K = constant 
f = distance to focus point 
L = distance to point under consideration 
AL = AL f/7T R 2 
X = wave length 
R = radius of telescope aperture 


In this program the interval (f - Af, f + Af ) is observed along the line 

2 2 

of sight of the laser where Af = 2 X f / 7T R . Within this interval a finite 
number of equally spaced points are -sampled. The velocity at each point 
(Wj) (Fig. 3-23d) is multiplied times the weighting function (1^) (Fig. 3-23c) 
at that point; these weighted velocities (V^I^) are then summed for each 
velocity present to determine the returned signal spectrum (Fig. 3-23e). 


• Unfocused, Pulsed Laser Sensor 


For a pulsed sensor there is a weighting analogous to the previous dis- 
cussion for the focused sensor, except that the weighting function is triangular 
in range with the base of the triangle equal to pulsed length (see Appendix B). 
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Fig. 3-23 - Example of Physical Derivation of Signal Properties for l/3 
Meter Diameter, 1 Meter Separation Bistatic Configuration 
Focused at 3 km Range 
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The sensor simulation was designed to handle either pulsed or CW-focused 
configurations. The simulations of both types are listed in Table 3-2. 

• Preprocessing of Returned Signal 

Once the returned signal is received by the detector of the laser Doppler 
sensor, the frequency spectrum output of the detector must be processed to 
provide a value which is representative of the desired data from the spatial 
data point. (The spatial data point is the focal volume for a focused configura- 
tion or the pulsed width volume for a pulsed configuration.) Two basic processes 
have been simulated to decipher a single velocity value for each spatial data point. 
These processes involve determining: 

• Centroidal velocity, and 

• Maximum velocity existing in the spatial cell 

The centroidal velocity refers to the centroid of the velocity spectrum 
as it appears in the example of Fig. 3-23e. For short ranges where the 
vortex velocity signal intensity is as great as the ambient wind signal intensity, 
this type of processing is adequate. However, for longer ranges such that the 
vortex velocity signal intensity is very low in comparison to the ambient wind 
signal (as depicted in Fig. 3 -23e), the centroidal velocity is most inadequate 
for representation of the vortex flow field. 

Determination of maximum velocity above the noise level requires an 
adequate signal-to-noise ratio in order to distinguish between the portion of 
the velocity spectrum attributed to returned signal from the vortex and the 
signal attributed to noise. This type of preprocessing would provide a good 
representation of the vortex flow from the spectrum as illustrated in Fig. 3-23d. 

Figures 3-24 and 3-25 illustrate centroidal and maximum velocity above 
the noise level preprocessing, respectively. Both preprocessing configurations 
provide adequate data at the 100-meter range. However, as the range increases 
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to 300 meters, note that the measured centroidal (W) velocity drops approxi- 
mately 50% below the actual (T) velocity in Fig. 3-24 while the maximum velo- 
city above the noise level (W) of Fig. 3-25 remains at the level of the actual 
velocity (T). 

3.4,2 Computer Simulation Runs 

Many computer simulation runs were made during the contract. A num- 
ber of the runs, made during the early portion of the contract, involved multiple 
scans through the vortex and displayed sensor output versus sensor position 
and velocity contours from multiple scans. These runs are listed in Table 3.1 
and discussed in Appendix B. Another larger group of runs was conducted 
involving single scans through vortices with results plotted as velocity versus 
a coordinate referenced to the vortex itself. The second group of runs is listed 
in Table 3-2 and discussed further in Appendix C. 

3.5 ADDITIONAL EFFORT REQUIRED BEFORE THE LONG RANGE 
WAKE TURBULENCE MONITORING SYSTEM IS IMPLEMENTED 

Several investigations should be performed before the long range wake 
turbulence monitoring system is implemented. These investigations can be 
grouped according to system simulations, measurement programs, and hard- 
ware development and are listed below: 

• System Simulations 

• Simulate single -component sensor systems with ambient winds 
and turbulence, 

• Simulate two -component sensor systems with ambient (cross) 
winds and turbulence. 

• Simulate sensor systems with true axial velocities, after experi- 
mental axial velocity data are available. 

• Develop decision-making philosophy and simulate entire system 
including sensors, computer, and displays to determine system 
effectiveness . 

• Determine the sizes of corridors into air terminals from flight 
path statistics of aircraft using these corridors. 
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Fig. 3-24 - Data Preprocessed Using Output Centroid of Return Spectrum 
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OATA POINT POSITION WITH RESPECT TO SITE LINE IN NEtERS 


Fig. 3-25 - Data Preprocessed Using Output Maximum Velocity (above noise level) 
within Spectrum of Returned Signal (same scan as Fig. 3-24 but dif- 
ferent data preprocessing) 
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Measurement Programs 


• Measure axial flow velocities of full-scale aircraft vortices 
with remote sensors. 

• Experiment with a single -component laser Doppler velocimeter 
detecting vortices in an air terminal environment. 

• Experiment with a two -component laser Doppler velocimeter 
detecting vortices in an air terminal environment. 


Development Programs 


• Develop MIL-Spec laser velocimeter hardware to meet high 
reliability standards. 


3-55 


LOCKHEED - HUNTSVILLE RESEARCH & ENGINEERING CENTER 



LMSC-HREC D225436 


Section 4 

HYBRID WAKE VORTEX MONITORING SYSTEM 

4.1 INTRODUCTION 

The hybrid wake vortex monitoring system relies less on sensors than 
on the detective system and less on prediction than on the predictive system. 

It incorporates some features of both as a system which should provide fore- 
casts with more accuracy than the purely predictive system without quite as 
much sensor sophistication as the detective system. Since the sensors of the 
hybrid system are not as sophisticated as those of the detective system, with 
maximum effort the system could become field -operational within a very few 
years. Although the predictive system does not provide quite the air terminal 
capacity increase that the detective system does, it is expected to enhance the 
terminal capacity over the purely predictive system because of the improved 
accuracy of the vortex hazard forecasting. 

4.2 TYPICAL HYBRID WAKE TURBULENCE MONITORING SYSTEM 

A typical predictive wake turbulence monitoring system is depicted in 
Fig. 4-1. This artist's concept of the system depicts both laser Doppler and 
acoustic sensors performing the measurement function. The laser velocim- 
eter measures vortex size and movement in vertical planes normal to the 
runway axis and located several hundred or thousand feet beyond the end of 
the runway. The acoustic sensor locates the vortex position in a plane nearer 
the runway end where the range requirement is less. System sensory mech- 
anisms are controlled by the system computer which also processes sensor 
data and performs the hazard forecasting function. The hazard status is dis- 
played to the traffic controller who maintains spacing between aircraft accord- 
ing to prevailing vortex conditions. These various functions and the equipment 
for performing them in the hybrid system are discussed in the following para- 
graphs. 
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4.2.1 Detection of Current Vortex Locus 

A single -component laser Doppler velocimeter should be capable of 
readily detecting and determining the tangential velocity distribution of a vortex 
if the system is configured to view the vortex perpendicularly to the vortex axis. 
This premise was discussed in Section 3.3.1, and is illustrated in Figs. 3-10 
and 3-11. The sensor could either be under or to the side of the vortex (i.e., 
under or to the side of the flight corridor) and would accurately measure the 
tangential flow field of the vortex. 

A typical air terminal layout with a hybrid laser Doppler wake turbulence 
sensor system is depicted in Fig. 4-2. The wind is shown in this figure blowing 
from left to right. Therefore, for both arrival and departure, aircraft would be 
traveling into the wind; i.e., from right to left. The approach corridors for this 
wind condition would be at the right (in the illustration) end of the field and the 
departure corridors at the left. The boundaries above and below ILS approach 
corridors are typically at 3 and 2 deg, respectively. The horizontal spread is 
2 2 ” deg to either side of the runway. Since planes descend onto the air terminal 
at widely varying vertical angles independent of the ILS corridor, the boundaries 
on the departure corridors include a larger vertical slice of air space (2 to 
16 deg) and the same horizontal spread (Zj deg to either side of the runway). 

The NASA report of Ref. 27 discusses feasibility loci for vortex sensors 
for the hybrid wake turbulence monitoring system. A first cut at locating 
the sensors has indicated that the vertical scan planes at 5000 feet beyond 
each end of the runway would provide much useful vortex position data to feed 
into the wake turbulence forecasting computer. Sensors at these ranges are 
depicted in Fig. 4-2. By using Fig. 4-2 and simple geometry, typical maxi- 
mum ranges for the sensors to scan the overhead planes adequately can be 
determined. The sensors located 5,000 feet from the runway end and scan- 
ning the arrival corridors should have a maximum range of approximately 
1000 feet for maximum glide slope of 6 deg. Sensors located 5,000 feet from 
the runway end and scanning the departure corridor should have a maximum 
range of approximately 2500 feet. From the requirements subsection of 
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Fig. 4-2 - Typical Physical Layout of Hybrid Wake Turbulence Monitoring System 
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Section 1 of this document, it was indicated that the altitude of an aircraft 
above the ground greatly influences the controllability of the aircraft after it 
encounters a wake turbulence. A 2500 -feet altitude above the ground could be 
used as a typical maximum altitude for a hybrid system monitoring capability. 

Figure 4- 3a illustrates a vortex scanner performing two or more near- 
vertical scans to determine the angular slope of the vortex axis between the 
two scan planes. High spatial resolution sensors would be required to 
perform this vortex axis slope determination unless the scan planes are 
widely separated. 

Figure 4-3b illustrates the use of a conical scan from a single -component 
laser Doppler velocimeter to provide a three -component average velocity pro- 
file. The laser beam is focused at a specific range and rotated through a to 
scan a cone of angle (3 . The average vertical wind component will appear 

on the returned signal as V ; the average (dc) value of the signal. As the 
sensor scans into and out of the horizontal wind component, a sinusoidal 
output will be detected with the peak of the sine wave at the apposition indi- 
cating the direction of the horizontal wind. The average horizontal wind velo- 
city can then be calculated by dividing the peak-to-peak amplitude of the sensor's 
velocity output by twice the sine of (3, the conic angle. 

At heights of up to a couple of hundred feet and at horizontal distances 
of 800 to 1000 feet, acoustic sensors have proven to be effective in tracking 
vortices (Refs. 28 and 29). The acoustic sensor vortex tracking configura- 
tion developed by the Department of Transportation's Transportation Systems 
Center is depicted in Fig. 4-4. The transmitter sends out an acoustical pulse 
of energy. The path of this pulse is altered by the flow fields of the vortices. 
One vortex bends the acoustic path upward; the other bends the acoustic path 
back toward the ground. The vortex which bends the acoustic path toward the 
ground is tracked by measuring the time delays of the acoustic pulses in reach- 
ing several receivers (1 and 2 of Fig. 4-4). Data of vortex tracks taken with 
this type system are shown in Fig. 4-5. This type sensor is designed to track 
vortices and not measure the velocities of their flow fields. 
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Plane B 



Velocimeter Van 


Fig. 4-3a - Determination of Vortex Loci Within Scan Planes (Multiple 
planes used to determine vortex trajectory) 
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(0 to 360 ) 


Laser Beam is steered 
through circular scan; 
beam moves to map out 
cylinder at Focal Volume 


» i 

\r 



Beam position on circular scan (a) 


= vertical wind component 
= V^, average wind velocity 
measurement 

V Q = horizontal wind component 
= V /z sin(3 

p ./ -K 

a, = direction of horizontal wind 
h 


Fig* 4-3b - Determination of Average Wind Profile from Conical Scan 
with Single-Component Laser Doppler Velocimeter 
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Fig. 4-4 - Pulsed Acoustic Bistatic Radar (Ref. 29) 



LOCAL TIME IOOI 

RADAR ALTITUDE 200 FT 
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4.2.2 Forecasts of Future Hazard 


The computer of Fig. 4-1 receives the inputs from the sensors discussed 
in Section 4.2.1 and processes these data to provide a forecast of the vortex 
persistence in the active air corridors of the air terminal. The forecast from 
the hybrid system should be more accurate than that from the purely predictive 
system because of the increased accuracy of the input data. The real-time wind 
input capability of the hybrid system (from the conical scan of the laser Doppler 
sensor) would alone greatly improve the system accuracy in forecasting wake 
turbulence persistence over the "predictive system" since the wind inputs into 
the predictive system are provided as a forecast by the local aviation weather 
forecaster. The additional sensing of the vortex locus and movement within 
the vertical scan planes of Fig. 4-2 should further enhance the system ! s fore- 
casting capability. The computer must then predict the loci of the vortices 
between the scan planes and forecast their movements for several minutes. 
The computer then would compare the predicted vortex movement with the 
boundaries of the active air corridors to determine the safe separation 
distance required between successive aircraft about to move into the corridor. 
The primary improvement to be gained by going from this hybrid system to the 
"detective system" is that the requirement to predict vortex loci between the 
scan planes (the detective system would measure these directly) will be elim- 
inated, 

4.2.3 Interface with Air Traffic Controllers and/or Pilots 

From the computer the vortex information must be fed to the ATC and/or 
pilots involved in the hazard avoidance. This operation could be handled for the 
hybrid system in the same manner as for the predictive system (discussed in 
Section 2.5). The communication must be simple and precise in order not to 
confuse the user of the data. The computer should be used to make as many 
of the decisions as practical to reduce, not add to, the burden of the user. 
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4.3 HYBRID SYSTEM SIMULATION 

A digital computer program was developed to simulate a single -component 
LDV scanning overhead searching for vortices. The simulation first performs 
a coarse scan overhead in 0 (Fig. 4- 3a) and range to locate anomalies in the 
flow field. If anomalies are located, the simulation, after it makes the coarse 
scan, goes back for a fine scan in the regions of the anomalies. If vortices are 
present, the simulation maps the vortices' velocities, calculates their circula- 
tions and locates the centers of their cores. 

Figure 4-6 represents the output of the simulation after the first coarse 
scan. The asterisks mark the loci of the high velocity points within the scan. 

Figure 4-7 is a detail scan mapping of the regions of the high velocity 
anomalies located by the coarse scan. The mapping represents velocity in 
terms of the angular position (0 ) of the sensor with scans at a number of dif- 
ferent ranges. The ranges begin at 575 meters (bottom curve A) and are 
incremented in one meter increments out to 625 meters (top curve A). The 
velocity scales for each curve are offset by 8 meters/sec from the previous 
curve (starting at the bottom) so that the viewer can distinguish between the 
curves. High velocity activity can be depicted near the center of this figure. 

Figure 4-8 represents a contour mapping of the velocities of Fig. 4-7. 

The contours are plotted in terms of sensor angle (0 ) and range. The F 
contour is the outermost and represents velocities of approximately 5 meters/ 
sec. The E contour represents 10 meters/sec; D represents 15 meters/sec, 
and so on until the innermost contour A represents 30 meters/sec sensed 
velocity. 

Once the computer has stored the data of Fig. 4-8, computations are 
then made to determine peak velocity, circulation, etc. The computer picks 
the peak velocity at each range of Fig. 4-7 by applying a cubic least squares 
curve fit to the velocities measured at the specified ranges and differentiates 
the curve and tests for the maxima. The peak velocities at each range are 
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next fitted by a quadratic least squares curve to determine the peak velocity 
of the vortex for all ranges. This curve is differentiated as. before to deter- 
mine the peak. 

After the peak velocities within the vortex are determined (one positive 
peak and one negative peak), the computer determines the diameter of the 
vortex core (distance between these peaks) and the locus of the core center. 

The computer next determines the vortex velocity at a point, radii from 
the vortex center (Z\ radii arbitrarily chosen to provide velocity data from 
the region of low velocity gradients), and multiplies this value by Zttt where 
r is the distance from the center of the vortex to the point. This provides 
the circulation. (In a real-world situation the velocities should be integrated 
along a line integral at some radius, r, to accurately determine the vortex 
circulation.) 

Figure 4-9 is a map of the loci of the centers of the vortices in the ex- 
ample case of one vortex). The peak velocities, core radii and circulations, 
both actual and as detected, are displayed along the bottom of the figure. 

4.4 IMPLEMENTATION OF A HYBRID SYSTEM 

Personnel at NASA (Ref. 27) have been statistically studying the effects 
of varying spatial resolution upon the probability of predicting the vortex posi- 
tion accurately enough for hazard avoidance. The preliminary results of these 
studies have indicated that a spatial resolution of approximately 40 meters is 
adequate. From Fig. 3-3 it can be seen that a one-half meter diameter coaxial 
focused optical configuration will provide a spatial resolution of 40 meters at 
a range of 600 meters (~2000 feet). If this spatial resolution proves to be 
adequate, the hybrid system could be constructed around state-of-the-art 
one -half meter diameter coaxial sensors with minimum of development 
costs. Scanning systems to perform scans similar to those required for the 
hybrid system are presently being fabricated by Lockheed and NASA. NASA 
is planning to have data processing equipment designed and fabricated in the 
near future which would be very similar to that required for the hybrid system. 
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Because of these present NASA efforts, the development costs for the hybrid 
system would be primarily devoted to system integration and hardware re- 
liability. 


4.5 WIND SHEAR AND TURBULENCE MONITORING AS A BY-PRODUCT 
OF HYBRID (AND PROBABLY DETECTIVE) SYSTEMS 


Wind shear in approach corridors has become a problem of increasing 
interest to pilots in the past few years. An example of this problem from 
Ref. 30 is illustrated below. 


"During a two -hour period on 4 January 1971 there were nine 
missed approaches to runway 4R at JFK. Included in the missed 
approaches were B-747s, B-727s, a B-707, a C-990, and a DH-6. 
During this same period an accident occurred at LaGuardia Air- 
port where one of the probable causes was stated as 'The failure 
of the pilot to recognize the wind shear condition and compensate 
for it."' (Ref. 30). 


The conical scan configuration of Fig. 4-3b would provide wind profiles 
to 2000 feet. Wind shears of interest and concern to approaching pilots lie in 
this region. Thus, a logical by-product output of the hybrid system would be 
a wind shear forecast from the wind shear observation. The instrument would 
be readily available. The effort to develop such a by-product would be devoted 
to processing the LDV sensor output to provide a suitable format for display to 
ATC and/or approaching pilots. 


4-17 


LOCKHEED - HUNTSVILLE RESEARCH & ENGINEERING CENTER 



LMSC-HREC D225936 


Section 5 
CONCLUSIONS 

When the need for a device (or system in this case) becomes recognized, 
a number of questions immediately arise. In the case of the wake turbulence 
monitoring system, some questions which arise are: 

• Is a wake turbulence monitoring system feasible or within the 
near state of the art ? 

• Would a wake turbulence system be cost-effective to develop, 
operate, and maintain? 

• What type system should be developed ? 

• What are the logical steps to develop the system? 

These questions are considered in tiiis conclusion. Attempts are made to 
answer those questions which are answerable at the conclusion of this limited 
study. 

5.1 IS A WAKE TURBULENCE MONITORING SYSTEM FEASIBLE? 


The studies performed by Lockheed -Huntsville have indicated that a wake 
turbulence monitoring system is feasible . In fact, two of the three types of 
systems conceived could be developed to field -operational status within two to 
four years, and the third, a more sophisticated system could probably be de- 
veloped in three to six years. The first two types of the system would require 
a minimum of new technology and the third would require that concepts presently 
on the drawing boards be developed. 


Safety is of concern in system feasibility. This is especially true when 
something as exotic in the public's eye as the laser is involved. Laser radia- 
tion in the wavelengths regime 0.4 to 1.4 microns is considered hazardous 

- 6 2 

(Ref. 31) for occular exposures of 10 W/cm and above. Laser radiation at 10.6 
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micron wavelength (far infrared) is considered hazardous (Ref. 21) for ex- 
posures of 10 W/cm and above, clearly five orders of magnitude more 
intense than the visible radiation hazard level. The wake turbulence monitor- 
ing system should be feasibly designed to maintain transmitted radiation levels 
well below the hazard level in order to assure safe operation in the air terminal 
environment. 

5.2 WOULD A WAKE TURBULENCE SYSTEM BE COST EFFECTIVE? 

For a system to be cost effective, the benefits derived from the system 
should outweigh the costs to develop, operate and maintain the system. The 
benefits to be derived from such a system were outlined in the introduction to 
this report as well as projected dollar savings to airlines from reducing spacing 
behind "heavies." Increased revenues to terminals from increased capacity and 
thus increased landing fees were projected. An additional note was made of the 
approximately $36M in litigations involving the Government in claims for wake 
turbulence -related accidents. The true benefits to be derived from the devel- 
opment and installation of wake turbulence monitoring systems is difficult to 
establish accurately because the capability of the system to increase terminal 
capacity, reduce delays and improve safety is not accurately known. Some of 
these values should be obtainable to a higher degree of accuracy if recommended 
computer simulation studies of the entire systems are performed. Some prob- 
ably will not be accurately known until a prototype system is developed and 
operationally tested in the terminal for a reasonable period of time. 

Table 5-1 presents a highly speculative picture of development prototype 
system costs. The dollars depicted could easily be misestimated by a factor 
of 50% in either direction. The costs were based upon developing one system 
to handle a set of parallel runways at a single air terminal. Costs for indi- 
vidual systems will be reduced for quantity production of the units. 

A comparison of Fig. 1-2 of the report introduction with Table 5-1 
indicates that over a period of several years the system should be quite cost 
effective . More accurate cost figures effectiveness will become available 
from additional development studies. 
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5.3 WHAT TYPE SYSTEM SHOULD BE DEVELOPED? 

This report discusses three types of systems: predictive, detective and 
hybrid. The relative effectivenesses of the system configurations is discussed 
in the introduction and the relative expense and development time are listed in 
the table on page 1-14. The predictive and hybrid systems contain many of the 
same elements. The predictive system requires the least development time, but 
it is also least effective. A logical development path is to begin development of 
the predictive and hybrid systems simultaneously with the idea of installing the 
sensors of the hybrid system into the predictive system once they are developed. 
The development of the detective system could proceed at a normal rate. This 
approach would allow a system to be introduced — the predictive system — into field 
operation within the shortest possible time. The system could .be updated to 
become a hybrid system with improved performance. The detective system 
could be introduced into field operation when it became available some time in 
the future. 

Additional studies might indicate that either the hybrid or detective sys- 
tem should be developed, but not both. This limited study does not go into 
enough depth to make this determination possible. 

5.4 WHAT ARE THE LOGICAL STEPS TO DEVELOP THE SYSTEM? 

Steps to develop the system can be grouped according to the system — 
predictive, detective, or hybrid — to be developed. Some of these are listed 
below: 

5.4.1 Predictive System 
• Simulation Studies 

• Develop mathematical model of predictive system. 

• Determine flight path statistics of arriving and departing 
aircraft at specific air terminals. 
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• Introduce wind rose to system model to determine percentage 
of cases where separations could be reduced. 

Measurement Program 

• Test prototype system at active air terminal. 

Equipment Development 


• Develop display concept. 

• Develop computer/sensors/display integration. 


.2 Detective System 


System Simulations 


• Simulate single -component sensor systems with ambient 
winds and turbulence. 

• Simulate two -component sensor systems with ambient (cross) 
winds and turbulence. 

• Simulate sensor systems with true axial velocities after experi- 
mental axial velocity data are available. 

• Develop decision-making philosophy and simulate entire system 
including sensors, computer, and displays to determine system 
effectiveness . 

• Determine the sizes of corridors into air terminals from flight 
path statistics of aircraft using these corridors. 


Measurement Programs 


• Measure axial flow velocities of full-scale aircraft vortices with 
remote sensors. 

• Experiment with a single -component laser Doppler velocimeter 
detecting vortices in an air terminal environment. 

• Experiment with a two-component laser Doppler velocimeter 
detecting vortices in an air terminal environment. 

• Determine LDV performance in inclement weather. 
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• Development Programs 

• Develop MIL-Spec laser velocimeter hardware to meet high 
reliability standards, 

5.4.3 Hybrid System 

• Simulation Studies 


• Continue development of total system model including display, 
logic, etc. 

• Continue development of model which performs overhead coarse 
and detail scans and calculates circulation, etc. 

• Develop flight path statistics 


• Measurement Program 


• Setup single -component LDV system in terminal environment 
to perform hybrid -type scans during aircraft flybys; record 
and analyze data. 

• Determine LDV performance in inclement weather. 

• Perform conical scan with LDV to determine wind profile. 


• Equipment Development 


• Further develop system hardware to meet high reliability 
standards . 
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Appendix A 

NUMERICAL VORTEX MODEL 
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Appendix A 

A model was sought that would approximate the vortex flow fields generated 
by large transport aircraft. Consequently considerable emphasis was placed on 
recent flight investigations of vortex parameters using such aircraft, in particular 
those initiated as a result of priority action on the part of the Director, Flight 
Standards Services, in February 1970. Two of these investigations (Refs. 1 and 
2) were concerned with aircraft response during vortex penetration while the 
third (Ref. 3) concerned itself with the measurement of vortex velocity profiles 
using tower-mounted hot wire anemometers. The latter study is the most appli- 
cable for the purposes of this investigation since the aircraft were operated in 
terminal-type configurations, and the vortices were monitored at altitudes of 
zero to several hundred feet above the ground. Such data, however, suffer from 
considerable scatter primarily because of the vortex not intersecting the loca- 
tion of the instrumentation, in most instances, through the center of the core. 
Consequently, considerable error exists in the estimation of core diameters 
and peak velocities. Figures 1 and 2 (obtained from Ref. 3) present values of 
peak tangential velocity versus vortex age and core diameter, respectively, 
for a variety of current airline transport aircraft. 
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Fig. 1 - Summary of Peak Vortex Tangential Velocities vs Age for All Aircraft. 
Data from Various Models Superimposed (heavy lines) 
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Fig. 2 - Summary of Peak Vortex Tangential Velocities vs Estimated Core 

Diameters for All Aircraft. Data from Proposed Model Superimposed 
(heavy lines) 
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The approach taken in the choice of a model was to choose from a com- 
parison of the test data of Ref. 3 with various theoretical vortex models that 
predict both tangential and axial velocity components, that model which yields 
the most satisfactory agreement. Of particular interest is that of Ref. 5 which 
presents a numerical method for solving the equations for a vortex core. The 
technique solves the system of equations of motion for a steady axially symmetric 
spiraling motion of an incompressible fluid at large Reynolds number with the 
additional assumption of the boundary layer approximation. The results of an 
application of the technique to a trailing vortex are presented in Figs. 3 and 4, 
for a Boeing 72 7 and 747, respectively. Peak velocity magnitudes and vortex 
core diameters are in reasonable agreement with the data of Ref. 3 (summarized 
in Figs. 1 and 2). 

Deficiencies in the model are the: 

• Requirement of an assumption of an eddy viscosity, and 

• Assumption of a set of starting conditions for the computation. 

A value for the eddy viscosity is estimated from the correlation that exists 
between the circulation developed by the aircraft and the measured eddy vis- 
cosity coefficient deduced in Ref. 6 from various wind tunnel and flight investi- 
gations on vortex decay. Figure 5 graphically summarizes this correlation. 

The initial conditions for the computation were taken from the approximate 
theory of Newman (Ref. 7), which solves the linearized form of the set of equa- 
tions used in the computational study of Ref. 5. 

A comparison of the proposed model with various popular vortex models 
is given in Fig. 6, for a given set of conditions (a 727 at 1000 feet downstream 
of the aircraft). 
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Profiles of Axial Velocity in a Trailing Vortex 



Profiles of Circumferential Velocity in a Trailing Vortex 


Fig. 4 - Boeing 747 Vortex Characteristics 
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Fig. 6 - Comparison of Vortex Models 
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Appendix B 

LASER VELOCIMETER VOLUME SCAN 
SIMULATION PROGRAM 
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Appendix B 


B.l DESCRIPTION 


Program simulates one vortex and scans it with one laser velocimeter 
Output is air velocity as a function of position. 


B .2 CAPABILITIES 

1. Vortex can be in any position. 

2. Velocimeter can scan any size volume in any direction. 

3. The tangential and axial velocity profile of the vortex 
can be input with different values. 

4. For focus system aperture size can be selected; for 
pulse system pulse length can be selected. 

B.3 RESTRICTIONS 

1. No wind exists. 

2. Aerosol is assumed to be constant throughout the 
vortex. 

3. The vortex is assumed to be uniform along its 
length. 

4. The vortex core is assumed to be a straight line. 

5. No change in vortex position or strength is allowed 
with change in time. 

6. Infinite signal-to-noise (S/N) ratio is assumed. 

B.4 SPHERICAL COORDINATE SYSTEM OF LASER 


The laser's focal point orientation is represented in a spherical coordi 
nate system centered at the laser. 
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The transformation from the spherical coordinate system to the laser's rec- 
tangular coordinate system is 


X = r sin0 cos 6 
s s s s 


Y = r sinG sin <4 
s s s s 


Z = r cosG 
s s s 


B.5 SCANNING OF LASER 


The bounds on the scan are input 


4> . lower bound of <f> 


4 > 7 upper bound of <£ 

S L S 

0 . lower bound of 0 

si s 

0 7 upper bound of 0 

S La S 

r _ lower bound of r 
si s 

r 9 upper bound of r 

S La S 
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Also input are the number of equal divisions for each of the above intervals. 

The scanning starts at the lower bound of each parameter. The simu- 
lated sensor first scans across <j) then moves along 9 while it continues to 

s s 

scan (f> and finally it moves along r while stiLl scanning 4> and 9 # For 
s s s s 

each point of the scan the simulation computes the velocity component of the 
vortex in the direction of the laser. 

B.6 COMPUTATION OF TANGENTIAL VELOCITY COMPONENT OF 
VORTEX IN DIRECTION OF LASER 

First the focal point of the laser is transformed from the spherical 
system to the rectangular system. 

X = r sinG cos6 
s s s Y s 

Y = r sinG sin<6 
s s s r s 

Z = r cosG 
s s s 

OS = (X , Y , Z ) focal point in laser's 
s s s 

rectangular system 


Z 



X 
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OP = (X^, Y^ f Z^) point on vortex core axis (input) 
PQ = (A, B, C) direction of vortex axis (input) 


where 


Distance from S to PQ is D = 


SP x PQ 



SP 



Total tangential velocity magnitude of the vortex at point S, according 
* 

to Newman 


V 


r 

o 

T ” 27 rD 


1 - exp 


4is z J 


where 


is normal of tangential velocity 

r is circulation 
o 

is freestream velocity 

V is eddy viscosity 

z is vortex axial coordinate (input 
value is presently being used 
fpr all positions) 


t-i . 4W 
O 7TpVb 

W = V 
00 

V =k n/p 


Newman, B.G., "Flow in a Viscous Vortex," Aeronaut. Quart., Vol.X, 
May 1959, p. 149. “ 
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where 


W is total aircraft weight 

p is density of atmosphere 

V is aircraft velocity 

b is wing span 

p is coefficient of viscosity 

k is coefficient that relates laminar 
kinematic viscosity to eddy viscosity 
for a particular r 


The tangential velocity vector at S is 


y T _ y (SP x PQ) 

T SP x pq| 


The magnitude of the component of VT in the direction of SO is 
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VT • OS 
OS 


COMPUTATION OF AXIAL VELOCITY COMPONENT OF VORTEX 
IN DIRECTION OF LASER 


The magnitude of the axial flow component, according to Newman, is 



D 

o 

4?r p v z 


exp 


W D 2 \ 

_°2 ) 

4i;z / 


where 

V A is norm of axial velocity 

D q is profile drag of the generating aerofoil 


The axial velocity vector at S is 


VA 



PQ 

IpqI 
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B.8 COMPUTATION OF VELOCITY INTERPRETED BY A COAXIAL 
FOCUSED LASER DOPPLER SYSTEM 


The laser does not make a point measurement. It effectively samples 
a finite volume of space along a line and consequently observes a variety of 
velocities. The signal is the strongest from the focal point and tapers off as 
the elemental volume considered increases its distance from the focus. 


The power contained in the photo-detector due to N particles per cubic 

centimenter is proportional to the current squared, as shown in the following 
* 

equation . 


.2 

1 


J_ 

4 


2 2 


.4 .4. 


T 7T J? OtoR A N 


00 

/ 


dL 


L 2 + 44 (£ _ L) 2 


X 2 f 2 


where 

i? is the quantum efficiency (electrons/photon) 

a is power level of local oscillator 

a is backscattering coefficient of particle 

R is radius of transmitter lens 
2 

A is total transmitted light flux in photons/sec 

3 

N is effective number of identical particles per (cm) 
f is nominal range of focusing 
X is optical wavelength 
L is range from transmitter lens 


The current squared due to an infinitesimal unit of length dL is 



(•£ 7T 9 2 a 2 o A 4 N^ 


R 

2 2 R 4 

L, 2 + 7T .. K - 
2 2 
X T 


(f-L)' 


jj; 

This is Eq. A. 8 on page A-6 of the Interim Report for Contract NAS8-25921, 
"Application of Laser Doppler Velocity Systems." 
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2 2 4 

assume 17, a, a, A, N are constants and let = 1/4^1? a a A N, then 



This is proportional to the equation used in the program. 



where (V Q jt/X) is a constant . 

Therefore, the equation in the program is proportional to the power in the 
detector due to a volume of space having an infinitesimal unit of length in the 
line of sight direction from the detector. 

In this program only the interval (f - Af , f + Af) along the line of sight 
of the laser is considered. 

Af = 2\f 2 A R 2 
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% 

This interval provides 50% of the signal. Within this interval a finite number 
of equally spaced points is sampled. The velocity at each point (V ) is multi- 
plied times the weighting function (1^) at that point; these points are then 
summed and divided by the sum of the weighting factors. 

t N 

V ! / V I 
n n' O n 

n=l n=l 


This V ^ function thus provides an LDV system output which is weighted accord 
ing to a calculated system spatial resolution. This type weighting function pro 
vides a data output similar to that which was recorded by using the spectrum 
analyzer where the one-dimensional LDV was field tested. 

B.9 COMPUTATION OF VELOCITY INTERPRETED BY A COAXIAL 
PULSED LASER DOPPLER SYSTEM 

The transmitted pulse is assumed to be a square wave. The detector is 
open for the same length of time as that of the transmitter; therefore, a pulse 
of length equal to the transmitted pulse enters the detector. 

The pulse length is also assumed to be short relative to the distance to 
the volume being observed. Therefore, the detected power of an elemental 
volume is directly proportional to the length of time its illumination passes 
into the detector. 

The power going into the detector relative to time is a square wave. 

But as a function of location of elemental volume, the curve looks like 

Power 



* 


Page 3-2 of the Interim Report. 
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"b" is the center of the pulse when the incoming pulse and the outgoing 
pulse exactly coincide. At this instant of time "a" is the back of the outgoing 
pulse and "c" is the front of the outgoing pulse. Figure B-l shows the outgoing 
and incoming pulses passing each other. The power curve shown above is then 
the weighting function (I n ) as in the focused case and is computed the same 
way. 

B.10 INPUT FOR LASER VELOCIMETER SIMULATION PROGRAM 

There are two sections of input of which the first is by way of namelist. 
The name of the namelist is (INPUT) and it has the following variables: 


Name 

Unit 

Description 

GAMMAC 

N-m-sec/kg 

circulation of vortex from gener- 
ating aircraft 

W8 

m/ sec 

freestream velocity 

UNU 

m / sec 

eddy viscosity 

Z 

m 

downstream distance of vortex 
from generating aircraft 

OP(3) 

m 

position vector of a point on 
the vortex axis 

PW (3) 

m 

direction vector for vortex 

PHI I 

deg 

starting PHI for sweep in PHI 

PHI2 

deg 

final PHI for sweep in PHI 

NPHI 

— 

number of evenly spaced PHI 

THETA1 

deg 

starting THETA for sweep in 
THETA 

THETA2 

deg 

final THETA for sweep in 
THETA 

NTHETA 


number of evenly spaced 
THETA in sweep 

DIST1 

m 

starting range for sweep in 
range 

DIST2 

m 

final range for sweep in range 

NDIST 


number of evenly spaced 
ranges in sweep 

PSI 


control magnitude of weighting 
factor. Set to one. 
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Name 


Units 


Description 


R 

m 

radius of telescope aperture 

FLAMDA 

m 

wave length of laser light 

PLENGT 

m 

pulse length 

DZERO 

N 

profile drag of the aerofoil 
(Used for axial velocity. If 
no axial velocity wanted, set 
DZERO to zero.) 

RHO 

kg/m 3 

density of atmosphere 

ICHOIC 

— 

flag for the following 


ICHOIC = 1 for focused con- 
tinuous wave laser 

ICHOIC = 2 for a non-focused 
pulsed laser 

ICHOIC = 3 for perfect velocity 
detection 


NWEIGH 

IOUT 


IPLOT 


IP LOTI 


half the number of evenly spaced 
points used in the illuminated 
volume for weighting 

flag for the printed output 

IOUT = 1 for no printed output 

IOUT = 2 for printing only 
velocity specified by IPLOT 1 

IOUT = 3 for printing all velocities 

flag 

IPLOT =0 no SC4020 plots 

IPLOT = 1 SC4020 plots 

flag indicates which velocities 
are to be plotted 

IPLOT 1 = 1 plots centroid of 
weighted velocity 

IPLOT 1 = 2 plots max-min or 
max (ABS(VL)) as specified by 
IVMAX 

IPLOT 1 = 3 plots total line of 
sight velocity at focus or at 
center of pulse volume 
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Name 


Units 


Description 


IPLOT 1 IPLOT 1 = 4 plots line of sight 

component of tangential velocity 
at focus or at center of pulse 
volume 

IPLOT 1 = 5 plots line of sight 
component of axial velocity at 
focus or at center of pulse 
volume 

IPLOT 1 = 6 plots normal of 
vortex velocity vector at focus 
or at pulse center 

IPLOT 1 = 7 plots normal of 
tangential velocity vector at 
focus or at pulse center 

IPLOT 1 = 8 plots normal of 
axial velocity vector at focus 
or at pulse center 

IVMAX flag 

IVMAX = 0 program computes 
maximum minus minimum of 
line of sight velocity within the 
focal volume or within the pulse 

IVMAX = 1 program computes 
the absolute maximum of the 
line of sight velocities within 
the focal volume or within the 
pulse 


The second input is for the plotting control. Each card for this input 
controls a group of plots. There are three variables input on each card, 
which are: 


1. Column 10 

2. Column 20 


0 for no lines between plotted points 

1 for lines between plotted points 

1 indicates PHI will be plotted on the horizontal axis 

2 indicates THETA will be plotted on the horizontal axis 

3 indicates RANGE will be plotted on the horizontal axis 
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3. Column 30 


1 indicates PHI will be plotted on the vertical axis 

2 indicates THETA will be plotted on the vertical axis 

3 indicates RANGE will be plotted on the vertical axis 


A blank card terminates the plotting and returns control to the beginning of the 
program for another case. 


B.l 1 OUTPUT FOR LASER VELOCIMETER VOLUME SCAN SIMULATION 

PROGRAM 


The input is printed on the printer and, in the case where plots are 
proudced, the input is also printed on the plotter. This is shown in Fig. B-2. 


NAME L I ST 


t A KMAC 

= 

C .*29000000 

03 

U8 


0.7730G000E 

02 

OP (1 > 

s 

0 . OOGOGCQOE • 

■3# 

OP 121 


0 . OOOOOOOOE 

-3d 

P« (1 1 

3 

0 . 1 OGCOOGQE 

01 

P« 12) 


0. OOOOOOOOE 

-3® 

PM11 

- 

Q.30Q0GQ0GE 

01 

PM J 2 


•0.30000000E 

Cl 

tmctai 


o.asgoqccoe 

02 

THE TA2 


0 .4 1999995E 

02 

CISTl 

* 

0.69C00DD0E 

02 

0IST2 


0-7 I C00C00E 

02 

P*l 

= 

0. 1DOOOOOOE 

01 

R 

= 

0 . 762 00000E 

-01 

c zz so 

= 

0 • 500DC00QE 

05 

RHO 

= 

0 .11920000E 

01 

IOUT 


3 


1 P LOT 

= 

1 



WNFUT) 


UNO 


0 . 9 1 4 000GGE -0 1 

2 

- 0. 1000000 0E 0* 

OP 13) 


0.50G000GGE 02 



F«13) 


0 •OOOOOOOOE -3d 



NFH I 


11 



NTHETA 


11 



NCI ST 

= 

11 



FLAMCA 


0.10600000E-04 

FLENCT 

= 0 .OOOOOOOOE -3d 

ICHOIC 

= 

1 

NWE1CH 

s 12 

I PLOT! 

r 

1 

1 VHAX 

s 1 


Fig. B-2 - Example for Volume Scanning Program 


The printer output of the velocities is controlled by the input variable 
IOUT. Figure B-3 is an example where IOUT =3. In this figure 

Distance = 70 represents the range of the focal point 
for the following measurements. 

The numbers following (PHI =) represent the scan angle PHI in degrees 
for the measurements in the column below the angle. 
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2.4 
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1.2 
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-0.0 

-0.6 
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-1.8 

-2*4 

-3,0 

nrr. 7 

THETA *t8.0 

-14.3 

- 1 4 . 6 

-14.7 

-i3.sr 

-riTc 

-7 .T~ 

-2.3 

2.6“ 

6.9 

9.9 

DVL 

26.8 

3 l • 4 

35.4 

37.0 

34.5 

27.1 

16.1 

4.8 

8 , 8 

14.1 

1 7 • 7 

VL 

- 1 3 t 3 

-12.9 

-11.8 

-9.9 

-7.0 

-3.3 

0.8 

4.8 

8.1 

10»S 

1 2 . 0 

VTL 

-12.6 

-11.7 

-10.0 

-7.3 

-3.9 

0.0 

3.9 

7.3 

10.0 

11.7 

1 2 • 6 

V AL 

-0.7 

-1*2 

-1.9 

-2.6 

“-3. 1 

-3.3 

3 • r 

-2.6 

-1.9 

-1.2 

• 0.7 

V 

27.1 

?8.9 

30.5 

31 .6 

32.3 

32.6 

32.3 

31 .6 

30.5 

28.9 

27.1 
~2 7 . 1 

VT 

27.1 

20*? 

30.4 

31.4 32 . r 

3 2T3 - 

”3ZV1 

31.4 

30.4 

28.9 

VA 

0.9 

1 .6 

2.5 

3.4 

4.1 

4.4 

4.1 

3.4 

2.5 

1.6 

0.9 

RADIUS 

4.2 

3.8 

3.6 

3.3 

3 .2 

3.2 

3.2 

3.3 

3.6 

3.8 

4.2 

THETA 

-16.5 

• 1 7 ♦ 3 

-17.5 

- 1 6 • 5 

- 1 3 . 8 

-9.3 

-3.7 

2.1 

7.2 

1 0 • ® 

1 2.9 

DVL 

27.5 

32.0 

35.7 

37.0 

34.2 

26.’ 

16.0 

4.7 

8.7 

14.0 

I 7.8 

VL 

-16. 4 

- 1 6 . 5 

-lb. 9 

- 14.0 

- 10.8 

-6.3 

-1.1 

4.2 

0.6 

n •& 

l3.7 

VTL 

<•15.0 

-14.2 

- 1 2. 2 

•9.1 

--.if. 9 — oio mr~ 

“ 9.1 

12.2 

14.2 

1 5 . 0 

VAL 

"I 

-2.4 

-3.6 

-4.9 

-5.9 

• 6.3 

-5.9 

-4.9 

-3.6 

-2.4 

-1.4 

V 

29.** 

31*4 

33.0 

34.1 

34.8 

35.0 

34.8 

34.1 

33.0 

31*4 


VT 

29.3 

31*2 

32.6 

33.5 

33.8 

33.9 

33.8 

33.5 

32.6 

31.2 

2 9.3 

VA 

1.8 

3.2 

4.9 

6.7 

8. 1 

6.6 

8.1 

6.7 

4.9 

3.2 

1.8 

RADIUS 

3.8 

3.4 

3.1 

2.8 

- Zi7 

” "2.6 

" 2.7 

2.8 

3.1 

3.4 

3.8 

THETA 46.8 

-18.6 

- 1 9 . 8 

-20.3 

-19.3 

-16.4 

-11.5 

-5.1 

1.5 

7.3 

11.5 

13.9 

-"17.9 

DVL 

28.2 

3 2.5 

35.9 

~36.9 

3 3.7 

26.5 

15.8 

”4.7 

8.6 

13.9 

VL 

-19.9 

-20.9 

- 21.0 

-19.5 

-16.0 

-10.6 

-4.0 

2.7 

8.6 

12.8 

1 5 , 2 

VTL 

-17.6 

- l 6 . 9 

-14. B 

-11. 1 

-6.0“ 

OVO ‘ 

“6.0 

11.1 

14.8 

16.9 

|7.6 

VAL 

-2.4 

-4.1 

-6.2 

-8.4 

- 10 . 0 

-10.6 

- 10 . 0 

-8.4 

-6.2 

-4.1 

-2.4 

V 

31*4 

33.4 

34.9 

35.9 

36.4 

36.5“ 

“36.4 

3 5.9 

34.9 

33 • ** 

31.4 

VT 

31.2 

33.0 

33.9 

3*1.0 

33.7 

33.4 

33.7 

34.0 

33. v 

33.0 

31.2 

VA 

3.2 

5.6 

8.5” 

~ 1TV5T 

n.H~ 

lTTV - 

1 3 . 8 

l l .5T 

8.5 

5.6 

3,2 

RADIUS 

3.4 

3.0 

2.6 

2.3 

2.1 

2.1 

2.1 

2.3 

2.6 

3.0 

3.4 

THETA 46.2 

-20.5 

-22.0 

-22.7 

-21 i e 

— 18.7 

1 3 • 4 " 

-6.4 

”0.9 

7.3 

12*0 

14.7 

DVL 

28.8 

32*9 

36. 1 

36.6 

33.7 

26.3 

15.7 

4.6 

8.5 

13*7 

1 7 • 9 

VL 

-23.7 

-25.7 

-26.6 

-25.6 

- 22.0 

-15.9- 

-7.9 

~ 0 * 6- 

8 . 0 

13.4 

l«.s 

VTL 

-20.1 

- 1 9 . 5 

-17.3 

-13.1 

-7.1 

0.0 

7.1 

13.1 

17.3 

19.5 

20.1 

' VAL 

-3.6 

-6.2 

-9.3* 

-12 5T5* 

-- pr . 9~ 

^1579“ 


-12 

“-9.3~ 

~-6*2‘ 

—n* 

V 

33.0 

35.0 

36.2 

36.8 

37.0 

37.0 

37.0 

36.8 

36.2 

35.0 

32.0 

VT 

32.7 

33.9 

33.8 

32.5 

30.7 

29.8 

30*7 

3 2.5 

33*8 

33.9 

32.7 

VA 

5.0 

8.5 

12.9 

17.3 

20.7 

22.0 

20.7 

17.3 

12.9 

8.5 

5.0 

RADIUS 

3.1 

2.6 

2.2 

1.9 

U6 

1 .5 

1.6 

1.9 

2.2 

2 * 6 

3.1 

THETA 4 5.6 

-22.1 

-23*0 

-24.7 

-23.8 

- 20.5 

-14.9 

-7.5 

0.3 

7.2 

12*3 

1&.3 

DVL 

29.4 

33*5 

36.4 

36.8 

33.5* 

2«*. 1 15 .7 


-8,5” 

“13.5 

j7.9 

VL 

-2 7.2 

-30*1 

-31.8 

-31.3 

-27.7 

- 20 . v 

-11.7 

-1.8 

7.0 

13*6 

1*7.4 

VTL 

-22.3 

-21*8 

-19.4 

-14.8 

-8.0 

- OiO 

8.0 

“1 4.8 

19.4 

21.8 

“"22.3 
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-4.9 

-8.3 

-12.4 

- 16.6 

-19.7 

-20.9 

-19.7 

-16.6 

-12,4 

-8.3 

-4.9 

V 

3 4.2 

35*9 

36.6 

3 7 VI” 

37.0 

3**’ 

37.0 

37.1 

” 36.8 

35*9 

3 4.2 

VT 

33.5 

34.0 

32.5 

28.9 

24.6 

22*5 

2 *t . 6 

28.9 

32.5 

34.0 

33.5 

VA 

6.9 
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37.0 
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— 25.6“ 
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17.9 

VL 

-30.6 
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-36.6 

-32.9 
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V 
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37.0 

36.9 

36 • 6 
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34.0 
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29.7 

22.6 

1 2.2 

0.3 

12.2 

22.6 

29.7 

33.3 

34.0 

VA 
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1 5 . 0 

22. 1 

29.2 

-34.5 

" 36.5 * 

™34 i 5“ 

“29.2 

-22.1 

1S.0 
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Fig. B-3 - Example Output 
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The first number following (THETA) is the scan angle THETA in degrees 
for the data until the next (THETA) is encountered. The numbers other than 
the first are the centroid of the weighted averages. 

DVL is the maximum minus the minimum velocities observed within the 
focal volume. Or, if (IVMAX = 1), DVL is the maximum of the magnitude of 
the line -of-sight velocities. 

VL is the line -of-sight velocity at the focal point, 

VTL is the line -of-sight velocity at the focal point due to 
the tangential velocity. 

VAL is the line-of -sight velocity at the focal point due to 
the axial velocity. 

V is the norm of the total velocity vector at the focal point. 

VT is the norm of the total tangential velocity vector. 

VA is the norm of the total axial velocity vector. 

RADIUS is the distance from the axis of the vortex and the 
focal point in meters. 

There are two types of plots. The first plot (Fig. B-4) is velocity versus 
one of the scanning parameters such as PHI. In this type of plot there are a 
number of curves plotted with various letters. Each letter represents a different 
value of another scanning parameter such as THETA In Fig. B-4, A repre- 
sents THETA = 48.0 degrees, B represents THETA = 47.4 degrees and so forth. 

The second type of plot is a velocity contour plot (Fig. B-5). Each letter 
represents a particular contour as indicated at the top of the plot. The contour 
lines between the points must be drawn in by hand. 

B.12 EXAMPLE RUN FOR LASER VELOCIMETER VOLUME SCAN 
SIMULATION PROGRAM 

For an example, take a vortex generated by a Boeing 747 while it is land- 
ing and SCAN it about 1000 meters from the generating aircraft. Approximate 
values for the parameters of this vortex are 
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GAMMAC = 727.0 nj-m-sec/kg 
W8 = 77.3 m/sec 
UNU = 0.0914 m 2 /sec 
Z = 1000 m 
DZERO = 50,000 n 
RHO = 1,192 kg/m 3 

Let the vortex pass directly overhead and parallel to the X axis which 
is on the ground 

OP = (0, 0, 50) 

PQ = (1,0,0) 

The following is the scanning selection 

PHI 1 = 3 deg 
PHI2 = -3 deg 
NPHI = 1 1 
THETA 1 ;= 48 deg 
THETA2 = 42 deg 
NTHETA = 11 
DIST1 = 69 m 
DIST2 = 71 m 
NDIST = 11 

A CO^ gas laser is used, therefore 

FLAMDA = 0.0000106 m 

The value for PSI is a scale factor and, set it to one, 

PSI = 1. 


Let the laser Doppler system be a single, coaxial, continuous wave, 
focused system. Therefore 

ICHOIC * 1 
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The lens of the telescope will be about a half foot in diameter. There- 
fore 


R = 0.0762 m 


When the focal length is 70 meters, half the length of the focal volume is 


Af = 


2\f _ 2*. 0000106*70 

V R 2 tt*(.0762) 2 


= 5.7 meters 


To sample the velocity in the focal volume at least every half meter, 

NWEIGH = 12. 


To print out everything. 


IOUT = 3. 

To plot the centroid of the weighted velocity, 

IPLOT = 1 
IPLOT1 = 1 . 

Let the program compute the absolute maximum of the line -of- sight 
velocity in the focal volume 


IVMAX = 1 

The first input is a comment on the first 72 columns of the first card. 
Then comes the main input which is input by way of FORTRAN namelist. The 
name of the namelist is INPUT and the data are the above described variables. 

The final input to consider is the plotting control. Two sets of plots 
PHI versus THETA and RANGE versus THETA are wanted. Also, in both 
cases, lines will be drawn between the points. 
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B-18 


LOCKHEED - HUNTSVILLE RESEARCH & ENGINEERING CENTER 


LMSC-HREC D225936 




C**H?LE 'OR VOLUME SCANNING PROGRAM n 


CifTANCE - 7.COCOCX1C 


♦ 01 


4.3000X10* 04 


*♦01 


*♦01 


.♦01 


4 .1000 X10* 01 L 


■U rl Ua 1. fii-... B I Q 


0.4000X10 


♦ 00 


■tg-.U tt?. 


i=- u 


o r~ 


t.toooxio * 00 • 

PHI IN DECREES 


-E 0- 


-1 .0000X10 


♦of" 


~ "*~r 

i 


-o — 


4 

-41 


-0.4000X10 


♦ 00 


Fig. B-5 - Example Output 
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Fig. B-6 - Example Output 
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These input cards are 

Card No. 


1 

2 

3 


Col. 10 
1 
1 

blank 


Col. 20 
1 
3 

blank 


Col. 30 
2 
2 

blank 


A sample of the printer output from this example run is shown as Figs. 
B-2, B-3, B-4, B-5, B-6 and B-7. 


B.13 LASER VELOCIMETER VOLUME SCAN SIMULATION PROGRAM 
VARIABLES 


Name 

Units 

Location 

Description 

ABC (20) 

None 

PSCAN 

Plot label. 

AVT(3) 

m/ sec 

COMMON/AO 5/ 

Axial velocity vector at 
illumination point - (rectangular 
coordinate system). 

D 

m/ sec 

COMMON/A 15/ 

D is the distance from illumi- 
nated point of interest to axis 
of vortex. 

DD(20) 

m 

MAIN 

Distance of focus point or 
center of pulse from axis of 
vortex. 

DDIST 

m 

COMMON/A 12/ 

Distance between points in 
range scan. 

DDISTW 

m 

MAIN 

Step size used in sampling 
the illuminated volume. 

DDV 

m/ sec 

PSCAN 

Preliminary result in 
determining delta velocity 
between the curves on the 
graph of velocity versus a 
scanning parameter. 

DISTL 

rn 

MAIN 

Distance from laser to focus 


point. 
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Name 

Units 

Location 

Description 

DISTW 

m 

MAIN 

Distance from laser to point 
of interest in the illuminated 
volume . 

DIST1 

deg 

COMMON/A09/ 

Starting range. 

DIST2 

deg 

common/ A09/ 

Ending range. 

DL2 


MAIN 

A preliminary result for 
the weighting associated 
with a focus system. 

DPHIL 

rad 

MAIN 

DELTA step in PHI scanning. 

DPHIL1 

deg 

common/a 12 / 

Angle between focus points 
when scanning in PHL 

DPTD(3) 

deg ,deg 

,m COMMON/A 12/ 

DPTD(l) is same as DPHIL1 
DPTD{2) is same as DTHEL1 
DPTD(3) is same as DDIST. 

DTHEL1 

deg 

common/a 12 / 

Angle between focus points 
when scanning in THETA. 

DTHETA 

rad 

MAIN 

DELTA step in THETA 
for scanning. 

DUM 

None 

PSCAN 

A dummy variable used in 
calling arguments. 

DV 

None 

PSCAN 

DELTA velocity between the 
curves on the graph of velocity 
vs a scanning parameter. 

DZERO 

None 

common/a 14 / 

DZERO is the profile drag of 
airfoil from generating aircraft. 

EX 


VEL 

Intermediate value for tangential 
velocity. 

FACT1 


MAIN 

A preliminary calculation for 
the weighting associated with 
a focus system. 

FACT 2 


MAIN 

A preliminary calculation for 
the weighting associated with 
a focus system. 

FLAMDA 

m 

COMMON/ A20/ 

FLAMDA = wavelength of light. 

F L.DX (12, 3) 

None 

common/aii/ 

Stores labels for the plots. 

FLDY(12) 

None 

common/aii/ 

Stores labels for the plots. 

GAMMAC 

(N-m-sec 

)/kg COMMON/ A06/ 

Circulation about vortex. 

I 

None 

MAIN 

The PHI do-loop parameter. 
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Name 

I 

LA. 

IABC 

ICHOIC 


IERR 

n 

IIA 

IHNE 

[OUT 


Units 

None 

None 

None 

None 


None 

None 

None 

None 


None 


Location 


Description 


PSCAN I is a subscript used in incre- 

menting the parameter plotted 
on the horizontal scale. 


PSCAN 


PSCAN 

COMMON/ A20/ 


Used in calculation for DELTA 
velocity between the curves 
on the graph of velocity vs 
a scanning parameter. 

X coordinate for label used 
on contour plot. Example (A =) 

Flag for laser type 

ICHOIC = 1 indicates focused 
CW laser 


PSCAN 

MAIN 

PSCAN 


PSCAN 


ICHOIC = 2 indicates pulsed 
laser 

ICHOIC = 3 indicates perfect 
velocity detection. 

IERR is an error flag from a 
plot routine. 

The weighting do -loop parameter. 

Used in calculation for DELTA 
velocity between the curves 
on the graph of velocity vs 
a scanning parameter. 

Determines if lines will be 
drawn between points of SC 
4020 plots. 

ILINE = 0 indicates no lines. 


ILINE > 0 indicates lines. 

COMMON/A20/ Flag for printer output 

IOUT = 1 no velocities are 
printed. 

IOUT = 2 prints only velocity 
specified by flag 
IPLOT1. 

IOUT = 3 prints all available 
velocities 
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Name 

IP 

IPLOT 

IPLOTI 


ISYM(20) 

IV 

IV MAX 


Units Location 


Description 


None 


None 


None 


None 

None 


None 


PSCAN 


COMMON/A20/ 


COMMON/ A20/ 


common/a 11 / 
PSCAN 

COMMON/A 20 / 


Contains address of storage 
location to go to in order to 
accomplish commands of 
variables IX and IY. Used in 
"assigned go to." 

A flag indicating whether or 
not SC 4020 plots will be 
produced. 

IPLOT = 0 indicates no plots. 

IPLOT = 1 indicates plots 
will be made. 


Flag indicating velocities to 
be plotted. 

IPLOTI = 1 Plots centroid of 
weighted velocity. 

IPLOTI = 2 Plots MAX -MIN or 
MAX (ABS(VL)). 


IPLOTI =3 Plots line-of-sight 
velocity. 

IPLOTI = 4 Plots line-of-sight 
tangential velocity 

IPLOTI = 5 Plots line-of-sight 
axial velocity 


IPLOTI = 6 Plots velocity 

IPLOTI = 7 Plots tangential 
velocity 

IPLOTI = 8 Plots axial 
velocity 


Stores plotting symbols 


IV is a flag indicating if current 
velocity is higher or lower than 
contour that is being sought. 


Flag to indicate whether MAX- 
MIN of line of sight velocity in 
illuminated region is to be 
computed or MAX (ABS (line 
of sight in illuminated region)) 
is to be computed. 
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Name 

Units 

Location 

Description 

IX 

None 

PSCAN 

Determines horizontal axis 

IX = 1 horizontal axis is PHI 
IX = 2 horizontal axis is THETA 
IX = 3 horizontal axis is RANGE 

IY 

None 

PSCAN 

Determines vertical axis 

IY = 1 vertical axis is PHI 
IY = 2 vertical axis is THETA 
IY = 3 vertical axis is RANGE. 

IZ 

None 

PSCAN 

Indicates parameter which is 
constant for each plot. 

IZ = 1 PHI is constant 

IZ = 2 THETA is constant 

IZ = 3 RANGE is constant. 

12 

None 

MAIN 

A do-loop parameter. 

13 

None 

MAIN 

Indicates first velocity to be 
printed out when IOUT = 3. 
ICHOIC determines value for 13. 

14 

None 

MAIN 

Indicates second velocity to be 
printed out when IOUT = 3. 
ICHOIC determines value for 14. 

J 

None 

MAIN 

The THETA do-loop parameter. 


None 

PSCAN 

J is a subscript used in incre- 
menting the parameter plotted 
on the vertical scale. 

JABC 

None 

PSCAN 

X coordinate for label used 
on contour plot. Example 
(30) in label (A = 30) 

JCHOIC 

None 

MAIN 

Contains address of storage 
location to go to in order to 
accomplish command of 
variable ICHOIC. Used in 
"assigned go to." 

JOUT 

None 

MAIN 

Contains address of storage 


location to go to in order to 
accomplish command of vari- 
able IOUT . Used in "assigned 
go to." 
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Name 

Units 

Location 

Description 

JVMAX 

None 

MAIN 

Contains address of storage 
location to go to in order to 
accomplish command of 
variable IV MAX. Used in 
"assigned go to." 

K 

None 

MAIN 

The range do-loop parameter 


None 

PSCAN 

Do-loop parameter for plots. 
Two plots are made for each 
value of K. 

KCHOIC 

None 

MAIN 

Contains address of storage 
location to go to in order to 
accomplish command of vari- 
able ICHOIC. Used in "assigned 
go to." 

KK 

None 

PSCAN 

Number of current contour 
program is working on while 
producing contour plot. 

KOUT 

None 

MAIN 

Contains address of storage 
location to go to in order to 
accomplish command of 
variable IOUT. Used in 
"assigned go to. M 

L 

None 

PSCAN 

Used in calculation for DELTA 
velocity between the curves 
on the graph of velocity vs 
a scanning parameter. 

LAB(2, 3 ) 

None 

PSCAN 

Plot label. 

LABLE(8) 

None 

MAIN 

Contain labeling for printer 
output. 

LINE 

None 

MAIN 

Line count for current 
page of line printer. 

NDIST 

None 

COMMON/A09/ 

Number of evenly spaced 
ranges. 

NEW 

None 

MAIN 

Not currently used. 

NNTH 

None 

MAIN 

Number of lines of output 
on printer per range step. 
Determined by IOUT. 

NNX 

None 

PSCAN 

Same as NX except NNX 
is negative if lines are to 
be drawn between points 
on graph. 
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Name 

Units 

Location 

Description 

NPHI 

None 

common/ A09/ 

Number of evenly spaced 
PHIs. 

NTHETA 

None 

COMMON/ AO 9 / 

Number of evenly spaced 
THETAs. 

NV 

None 

PSCAN 

NV is number of positive 
contours to be plotted. 

NW 

None 

MAIN 

Number of sampling points 
in the illuminated volume. 

N WEIGH 

None 

C OMMON/ A 2 0 / 

NWEIGH is half the number 
of evenly spaced points used 
in the weighting of this illum- 
inated volume. 

NX 

None 

PSCAN 

Number of steps in parameter 
which are plotted along hori- 
zontal axis. 

NY 

None 

PSCAN 

Number of steps in parameter 
which are plotted along vertical 
axis. 

OP (3) 

m 

C OMMON/ AO l/ 

Position of a point on the vortex 
axis (rectangular coordinate 
system). 

OS(3) 

m 

common/ao 5 / 

Vector to indicate illumination 
point relative to laser (rectan- 
gular coordinate system) 

PHIA 

rad 

MAIN 

Starting PHI. Same as PHI1 
except for units. 

PHIB 

rad 

MAIN 

Ending PHI. Same as PHI2 
except for units. 

PHIL 

rad 

MAIN 

Current value for the 
angle PHI. 

PHIL1(20) 

deg 

common/ao 4/ 

PHI values to be plotted 
or printed. 

PHIL2 

deg 

MAIN 

Current value for the 
angle PHI. Same as PHI2 
except for units. 

PHI1 

deg 

common/ao 9 / 

Starting PHI. 

PHI2 

deg 

common/ao 9 / 

Ending PHI. 

PI 

None 

COMMON/A07/ 

PI is v = 3.1415926. 
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Name 

Units 

Location 

Description 

PLENGT 

m 

COMMON/ A20/ 

PLENGT is pulse length. 

PQ(3) 

None 

common/ao l/ 

Directional vector for 
vortex (rectangular coor- 
dinate system). 

PSI 

None 

C OMMON/ A2 O/ 

PSI is incident wave 
amplitude. 

R 

m 

C OMMON/ A2 o/ 

R is radius of telescope 
aperture. 

RHO 

kg/ 

COMMON/a 14/ 

RHO is the density of the 
atmosphere. 

R2D 

deg/ rad 

COMMON/AO 7/ 

Conversion factor for 
converting from radians 
to degrees. 

SP(3) 

m 

COMMON/AO 5/ 

Vector from illumination 
point (OS) to point (OP) on 
vortex (rectangular 
coordinate system) 

THETAA 

rad 

MAIN 

Starting THETA. Same 
as THETA1 except for 
units. 

THETAB 

rad 

MAIN 

Ending THETA. Same as 
THETA2 except for units. 

THETAL 

rad 

MAIN 

Current value for the angle 
THETA. 

THETA1 

deg 

COMMON/AO 9/ 

Starting THETA. 

THETA2 

deg 

common/ao 9/ 

Ending THETA. 

THETL1 (20) 

deg 

COMMON/A04/ 

THETA values to be plotted 
or printed. 

TITLE (12) 


COMMON/AO 2/ 

Contains inputted description 
of run. 

UNU 

2/ 

m / sec 

common/ao 6/ 

Eddy viscosity. 

V (20, 20, 20) 

m/ sec 

COMMON/A03/ 

Velocity values to be plotted: 
First index varies PHI. 
Second index varies THETA. 
Third index varies range. 

VC ON 

m/ sec 

PSCAN 

VCON is velocity of contour 
to be plotted. 
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Name 

Units 

Location 

Description 

VC ONI 

m/ sec 

PSCAN 

Velocity of maximum 
contour to be plotted. 

VL 

m/ sec 

MAIN 

Line -of- sight velocity at 
point of interest. 

V LINE 1 

m/ sec 

C OMMON/A 15/ 

VLINE1 is the line-of-sight 
component of the tangential 
velocity. 

V LINE 2 

m/ sec 

COMMON/A 15/ 

VLINE2 is the line-of-sight 
component of the axial 
velocity. 

VLW 

None 

MAIN 

Sum of line-of-sight velocity 
multiplied by weighing factor 
over the illuminated volume. 

VMAVT 

m/ sec 

C OMMON/A 15/ 

VMAVT is the norm of the 
tangential velocity. 

VMAVTD 

l/ sec 

VEL 

Scalar used to produce the 
axial velocity vector AVT(3) 
from vector PQ(3). 

VMAX 

m/ sec 

PSCAN 

Maximum value for velocity 
of current plot being made. 

VMAXX 

m/ sec 

PSCAN 

The maximum velocity on 
the graph of velocity vs a 
scanning parameter. 

VMIN 

m/ sec 

PSCAN 

Minimum value for velocity 
of current plot being made. 

VMPQ 

m 

VEL 

Norm of vector PQ(3). 

VMVT 

m/ sec 

COMMON/A 15/ 

VMVT is the norm of the 
tangential velocity. 

VMVTD 

l/ sec 

VEL 

Scalar used to produce the 
tangential velocity vector 
VT(3) from vector VV(3). 

VMVV 

m 

VEL 

Norm of vector VV(3) 

VP 

None 

PSCAN 

VP is used in computing 
a linear interpolation 
between data points, 
for the contour. 

VT(3) 

m/ sec 

C OMMON/A 05/ 

Tangential velocity vector 
at illumination point (rec- 
tangular coordinate system). 
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Name 

Units 

Location 

Description 

VTOL 

m/ sec 

MAIN 

Norm of velocity vector at 
focus point. 

VV(3) 

None 

COMMON/ AO 5/ 

Vector pointing in direction 
of tangential flow at illum- 
ination point (rectangular 
coordinate system). 

VV (20, 20) 

m/ sec 

common/a 10 / 

Values of VV are obtained 
from V and stored in proper 
order for plotting. 

VVMAX 

m/ sec 

MAIN 

Maximum line -of- sight 
velocity in the illuminated 
volume. 

VVMIN 

m/ sec 

MAIN 

Minimum line -of- sight 
velocity in the illuminated 
volume. 

VVP 

deg or m 

PSCAN 

VVP is one coordinate on 
a contour. 

VVV(20, 20) 

m/ sec 

common/a 10 / 

Values of VVV computed 
from VV which has had a 
DELTA velocity added to it 
in order to separate the 
graphs on the velocity plot. 

V 1 (20, 8) 

m/sec 

MAIN 

Stores velocities for one 


sweep in PHI. 

Vl( , 1) is centroid of 
weighted line of sight 
velocity. Vl( ,2) is MAX- 
MIN or MAX (ABS(VL)). 

Vl( , 3) is line of sight 
velocity at focus point. 

VI ( , 4) is line of sight of 
tangential velocity at focus 
point. VI ( ,5) is line of 
sight of axial velocity at 
focus point. Vl( ,6) is 
norm of total velocity. 

VI ( , 7) is norm of tan- 
gential velocity. Vl( ,8) 
is norm of axial velocity. 

WAT None MAIN Intensity of pulse. Presently 

set to one. 


G 
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Name 

Units 

Location 

Description 

WF 

None 

MAIN 

Weighting factor for point 
of interest. 

WFT 

None 

MAIN 

Sum of weighting factors 
over the illuminated volume . 

W8 

m/ sec 

COMMON/AO 6/ 

Freestream velocity. 

2 

m 

COMMON/AO 6/ 

Vortex axial coordinate from 
source. 
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Fig. B-8 - Flow Chart of Laser Velocimeter Volume Scan 
Simulation Program 
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B. 14 LASER VEL.OCIMETER VOLUME SCAN SIMULATION PROGRAM LISTING 


*JOP 

*SETUP LR ? 
PAPE I GN 
*EXECUTE 
* t R JOE LMSG 
TIPFTC MAIm 
C 


TAPF,SC4O?0 

SYSLB? 

IB JOB 

MAP 

DFCK 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

r 

c 

c 

c 

c 

c 

c 


c 

r 

r 

r 


LASER velocimftfr volumf scan simulation program 


this program simulatfs one vortex and 

VELOr I METER • OUTPUT ON SC40?0 PLOTS ARE 
OF POSITION* 


SCANS IT WITH ONE LASER 
AIR VELOCITY AS A FUNCTION 


this program computes line oe sight 
WFIGhTFO lime of sight velocities. 

SC40 ?r\ plotting IS PERFORMED 
VELOCITY CONTOURS 


NAMELIST (INPUT) 


velocities and 


AAMM AC 

CIRCULATION 


W8 

FRFFSTREAM VELOCITY 


UNI 1 

FDDY viscosity 


7 

VOPTFX AXIAL COORDINATF 

FROM SOURCE 

OP 

POSITION OF A POINT ON 

THE VORTEX AXIS 

PO 

DIRECTIONAL VICTOR FOR 

VOPTFX 

PHI 1 

STARTING PHI 


D H I p 

ENDING PHI 


N°H I 

NUMRPP OE EVENLY SPACED 

PHI S 

THETA 1 

starting theta 


THETA P 

ENDING THETA 


K'THFtA 

NUMRER of FVENLY EPACED 

thft a s 

01 ST 1 

STARTING RANGE 


0 ISTO 

fnding range 


NO I ST 

NUMBER OF EVENLY SPACED 

PAnGFS 

PS I 

INCIOFNT <VAVF AMPLITUDE 


R 

RADIUS OF TELESCOPE APFRTURF 

FLA^pfl 

WAyF LFNGTH OF LIGHT 


PLENPT 

PUL SF LFNGTH 


r>7ER^ 

PROFILE DRAG OF THE AEROFOIL 

RHO 

DENSITY OF ATMOSPHERE 


ICHn T c 

FLAG 



NWE I r.H 


I Ot IT 


T PL^T 


I PL 


ichoic = 1 

ICHOIC = * 

ichoic = o 

HALF THF MllMPFR 


continuous WAN/F laser 

PULSFn LAcrTQ 

PFPFFCT VELOCITY DFTF C T I ON 
OF EVENLY SPACFD POINTS ifcrn 


I N 


THF WEIGHT l MG 
FLAG 

I OUT = 1 

I OUT = 7 

I OUT = 0 

FLAG 

I PLOT = P 
1 PLOT = 1 

FL A G 

I PLOT! = 1 


NO VFLOCITYS APF PR I NTEP 
PRINTS ONLY VELOCITY SP<=T I F I FO 
RY I PLOT 1 

PRINTS ALL VELOCITIES 

NO SC4020 PL OT c 
cr 40 pn plot c 

PLOTS CENTRO I p OF uiEI^Hjfn VELOP T TY 
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r 

I PLOT 1 

= ? 

r* 

I PLOT 1 

rr 3 

r 

I PLOT 1 

= 4 

r 

I PLOT 1 

= 5 

r 

I PLOT 1 

= 6 

r 

I PLOT 1 

= 7 

r 

I PLOT 1 

= 8 

r 

IV/MAv flag 


r 

IVMAX 

= 0 

r 

IVMAX 

= 1 


r 


PLOTS M AX— M I N OR MAX(ABS(VL)) 

PLOTS LINF OF SIGHT VFLOC I TY • 

PLOTS LINF OF SIGHT TANGFNTIAL VFlOC IT' 
PLOTS LINF OF SIGHT AXIAL VFLOC I TY 
PLOTS VELOCITY 
PLOTS TANGFNTIAL VELOCITY 
PLOTS AXIAL VELOCITY 

COMPUTE MAX-MIN OF VL 
COMPUTE VFLOC I TY MAX ( ARS ( VL ) ) 


C 

r 

r 

r 

C 

C 

C 

C 

r 

C 

C 

c 


c 

r 


PLOTING INPUT 

COLUMN 10 O NO LINFS RFTWFFN POINTS 

1 LINFS RFTWFFN POINTS 

COLUMN 20 1 PHI PLOTFD ON HORIZONTAL AXIS 

2 THFTA. PLOTFD ON HORIZONTAL AXIS 

3 DISTANCF PLOTFD ON HORIZONTAL AX I A 

COLUMN sc 1 PHI PLOTFD ON VERTICAL AXIS 

2 THFTA PLOTFD ON VFRTICAL AXIS 
s DISTANCF PLOTFD ON VFRTICAL AXIS 

COMMON / AOl / OP<3> , POO) 

COMMON / AC? / TITLF(IP) 

COMMON / AOS / V(?0,20,?0) 

COMMON / A04 / PH I L 1 ( 2.0 ) , THFTL1(20>, D I STL 1(20) 

COMMON / A06 / GAMMAC* U/8, UNU, Z 
COMMON / ACT / PI, R?D 

COMMON / AOO / PH II, PHI 2* NPHI, THFT A 1 , THFTA?, NTHFT A , OIST1 

* nic. T ?, NO 1ST 

COMMON / Ais / DPH I L 1 , DTHFL 1 , DDlsj 

COMMON / A14 / DZFDO, RHO 

COMMON / A IF / VLINF1, VLINF?, VMVT, VMAVT , D 

COMMON / A20 / PS I , R. FLAMDA', PLFNGT, ICHOIC, NWEIGH, IOUT, 

T I PLOT* I PLOT 1, IVMAX 
D IMSms I ON VI (?0,R) ,DD (20) ,LARLF( 8) 

DATA LABLF / 6HWFIGHT,SHDVLOHVL *3HVTL« 3HVALOHV OHVT *3HVA 
DATA PI, P?D / 3.141FR?6, F7.29F78 / 

DATA GAMMAC » W8, UNU, Z» OP, PQ , PHIU PH 1 2 , THFTA1, THSTA2, 

* DISTl, D I S T 2 » NPHI, NTHFT A , ND I ST , NWEIGH, IOUT, I PLOT 

S / 18*0,, 3*0, 1,3,0 / 

n.ATA PSI, R, FLAMD A / 3*0. / 

DATA PLFNGT, OZFRO, RHO, ICHOIC, IPL0T1, IVMAX / 3*0., 3*1 / 

NAMFt 1ST / INPUT / GAMMAC, W8 , UNU, Z, OP, PO , 

1 PH I 1 , PHI?, NPHI, THFT A 1 , THFTA2* NTHFT A , DISTl, DIST2, ND 1ST 

2 °S I , P, FLAMDA, 

3 PLFnGT, D7SR0, RHO, 

4 ICHOIC, NWEIGH, IOUT, I PLOT , I PLOT 1, IVMAX 
1C RFAp <F,F00) titlf 

pCC FORMAT ( 1 2A6 ) 

WRITS (8, 800) TITLF 
800 FORMAT ( 1H1 ,30X, 1 2A8) 

I NPl IT 

PFAO (F, INPUT) 

"IP ITS (8, INPUT) 


/ 


r 
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C INITIALIZE 

NEW - 1 

PHI A = PHI 1 /RPO 
PH IB = PHIP/RPD 
THETAA = THETA 1/R2D 
THETAR = THETA£/R?D 
GO TO ( 1 1 . IP* 1 3) , I OHO I C 
11 A?dr,M 4 1 TO JCHOIC 
AES t ON 44 TO KCHOI0 

r- a to 14 

1? Accjr-N 4 P TO JCHOIC 
A SS y ON 4« TO KCHOIC 
14 13 = 1 

14 = p 

00 TO 1 E 

13 ASS I r.N BO TO JCHOIC 

13 = 3 

14 = 4 

1* OO TO (PI* PP ♦ 1 6 > ♦ I OUT 
PI ASS Toni 31 TO JOUT 
AES I CM 110 TO KOUT 
NNTH = O 
OO TO 17 

PP ASCIAM 4B TO JOUT 
ASS I ON BE TO KOut 
NINTH = ( 1 O- I 3 ) *NTHETA 
00 TO 17 

10 ASc; I ON 4 B TO JOUT 
ASSIGN B 6 TO KOUT 
NNTH = M THETA 

17 IF ( I \/VAX .OT.O) 00 TO 18 
ABE I ON PI TO JV/VIAX 
op rr 1 9 

IB ASS I ON 9 P TO JVNAX 
lo CONTI Ml IE 

PPH 1 1 = (PHTB - PH I A ) /FLOAT ( MPH I - 1) 

PTHFtA = (THETAR - THETA A ) /FLOAT ( NTHETA - 1) 

ODIST = (DISTP - D I ST 1 ) /FLOAT (NO I ST - 1) 

DIOTL = DIET1 - DPI ST 

LINE = lOO 

0PHIL1 = DPH I L*RPD 

DTHEL 1 = DTHFT A*R?p 

THETAL = THFTAA - OTHETA 

THETAP = TH r T AL*PPD 

PHIL = PH I A - DPHIL 

PHILP = PHTL*RPD 

e ACT 1 = ( PS I * D I *P**?1 **P/EL AMDA**? 

FACT? = E|_AKDA/(PI*R**3) 

DO Pn 13=1 , NPH I 
PHIL’ = PHI! P + DPHIL 1 
po PHIL 1 ( IP) = D HIL2 
DO I P=l, NTHETA 

THETAP = THETAP + DTHEL 1 
PE THETL. 1 ( I P > = THETAP 

r 

r range do loop 

DO 1 po K=1,NDIET 

THETAL = THETAA - DTHEtA 
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OISTL = OISTL + OD 1 ST 
nicri. 1 (K)=0TGTL 
GO Tn JOUT, (4R«31 ) 

4R IF (L IMF + NINTH .LT. 64) GO TO 30 
WP I Tc (6. 600) TITLF 
LINF = 1 

30 WR I Tf ( 6* 6?0 ) OISTL. (PHIL1 ( I?) . 13=1 *NPHI ) 

630 FORMAT < //61X. 1 OHDI STANCE =.F7.0/12H PHI = .20F6.1) 

LINF = LINF + 4 

r 

O THFTA OO LOOP 

3 1 no 1 1 n J=1 ( NTHFT A 

THFtal = TH C T AL + OTHFTA 
PHIL = PHIA - OPHIL 
C 

c phi oo loop 

OO 1 no 1 = 1 « NPH I 

WMAy = -1 * F30 

W/MIm = 1 .F30 

PHIL = PHIL + OPHIL 

GO Tn JCHO I C » (41 . 43 . 30 ) 

4 1 OOlSTW = 4 • *F ACT 3 *n I STL**2/FL0AT ( 2*NWE I GH ) 

GO TO 43 

4 P DDISTW = PLFNGT /FLOAT ( 2*NWE I GH ) 

4*5 OISTw = OISTL - OO I STW*FLOAT ( NWE I GH + 1) 

NW = 2* NWE I GH + 1 

VLw = O . 

vft =o. 

r 

0 RAMGF OO LOOP FOR I LLUM I N ATFD VOLUME 

OO on I 1 = 1 ,NW 
PI ST'.' = O 1ST W + DniSTU 
r 

c compute weighting factop 

GO Tn KCHOTC. ( 44,43) 

44 OL3 = ( fact?*o I STL*OI STW )**?. 

WF = F ACT 1 / ( 1 • + (OISTL - DISTW)**2/DL2)/OISTW**2 

GO m 4 6 
4 3 WAT = 1 . 

T F ( I T .GT. NW^IGH) GO TO 60 

WF = WAT *FLOAT ( I I — 1 ) /FLOAT ( NWE I GH ) 

GO TO 46 

60 WF = WAT*FLOAT ( NW- I I ) /FLOAT ( NWE I GH) 

r. 

r sum weighting factor 

46 WET s WFT + WF 
r 

r CALL FOP V/ELOCITV AT POINT OF INTEREST 

CALL V/FL ( VL.PHTL , THFTAL » D I STW) 
r 

r LOOK FOR FOCUS POINT AND SAVE VELOCITIES 

I F ( I T .NF. NWE I GH + 1) GO TO 80 

V 1 ( I , 3 ) = VL 

V 1 ( I , 4 ) = VL I NE 1 

V 1 ( T , 3 > = VL I NE2 

VTOI... = SORT ( VMVT**? +VMAVT **? ) 

V 1 ( I , 6 ) = VTOL 

V 1 < I , 7 ) = VMVT 
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c. 

r 

r 

c 

on 

r 

C 


r 

r 


n i 


C 

9 ? 


ro 


l no 
r 
o 

re 

ftin 

R6 


1 0 4 
61 1 

61 ? 
1 1 n 

1 ?n 


Vl(I,R) = VMAV/T 
00 ( 1 ) = n 

FIND iv AX I Ml )M AND MINIMUM OF LINE OF SIGHT 
I F ( V|_ .GT. VVMAX) VVMAX = VL 
IF (VI .LT. VVMIN) VVMIN = VL 

SUM WEIGHTING FACTOR TIMFS LINE OF SIGHT VFLOC I TY 
VL'" = VL W + VL*V.'F 

COMPUTE CFNTPOIO OF VELOCITY CURRENT IN DETECTOR 
VI (1,1) = VLW/WFT 

GO To JVMAX , ( R 1 , 9? ) 

COMPUTE MAXIMUM M t NU C MINIMUM 

V 1 ( I , ? ) = VVMAX - VVMIN 

GO to ion 

COMPUTE ARSOLUTE MAXIMUM 
VHIi?) = AMAX1 (VVMAX, -VVMIN) 

co To inn 

CALL VFL(V1 ( I ,0) , PH I L , THFT AL , D I STL ) 

VI ( I ,4 ) = VL INE1 

VI ( I , E) = VLINFF 

VTOL = RORT ( VMVT**F +VMAA/T**?) 

VI ( I , 6 ) = VTOL 

Vl ( I ,7) a VMVT 

V 1 ( I , R ) = VMAVT 

DD ( I ) = O 

Vtl.J,*) = V1(I,IPL0T1) 

PRINT OUTPUT ON LINE PRlNTEP 
GO TO KOlIT, (RR,F 6 , 1 10) 

1*1 R I Tr ( 6 , 6 1 O ) THETL 1(J),(V(I1,J»K)*T1=1, NPH I ) 

FORMAT ( 6 H THFTA ,F 6 . 1 ,F0F6. 1 ) 

GO TO MO 

HR I Tp ( 6, 61 o ) THRTL 1 ( J ) , ( Vl ( I 1 , I R ) , I 1 = 1 , NPH I ) 

DO 1 04 I ? = I 4 , R 

HR I Tr ( 6 . 6 1 1 ) LARLE (IF), (\/ 1(11, IF), 11 = 1, NPH I ) 

FORMAT ( 6 X, A6,F0E6. 1 ) 

HP I TR (6,61?) ( OD ( I 1 ) , I 1 = 1 , NPH I ) 

FORMAT ( 6X«6HRA0IUS,?0F6. 1 ) 

C.ONT 1 NUF 

LINE = LINE + NNTH 
CONTINUE 


C CALL PLOT I NO PROGRAM 

I E ( I r>L OT .GT. O) CALL °SCAM(NFW) 

NFni = n 

r 

r RETURN TO FTART FOP NEXT CAFF 

GO TO in 
END 
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* ! PPTr. \/EL nFCM 

FUBPPUT I MF VFL ( VL I NF % PH I L » THFT AL « D I STL ) 
r 

r. THIS POUTTNF COMPUTE TANGENTIAL AND AXIAL VELOC I TYS OF A 

G V/OPTFy AT A POINT OF INTRFST. ALSO THE LINE OF SIGHT COMPONENT 0 ' 

r rPOM THF LASTFP APF COMPUTED 


r 

r 


r 

r 


r 

r. 


r 

r 


r 

c 


r 

C 


r 

C 


c 

c 

r 

r 


COMMON / AD] 
COMMON / ftfiR 
COWWOM / AOS 
COMMON / A P 7 
COMMON / A 1 h 
COMMON / A 1 s 


/ OP ( 3 ) « PO ( 3 ) 

/ OF ( 3 > » SP ( 3 ) » VV ( 3 ) » VT ( 3 ) 
/ GAMMAC* "IRi UNU» 7 
/ PI* RPD 
/ D7FP0, OHO 

/ VL I NF 1 « VL I NF3 « VMVT * VM A VT < O 


* AVT( 3) 


os IF VECTOR FROM LASER TO POINT OF INTRFST 
OS< n = DISTL*FTN(THfTAL)*COS(PHIL) 

0°(?| = niSTL*FIN(THFTAL)*SIN(PHIL) 

OS ( 3 ) = ni STL*COF ( THFTAL ) 


SP IS VECTOR FROM POINT OF INTRFST TO POINT ON VORTEX 
FP( 1 ) = OP( 1 ) - OS( 1 ) 

SP( ?) = OP(?) - o° ( ? ) 

= OP (3) - OS (3) 


VV IF VFCTOR IN DIRECTION OE TANGENTIAL VELOCITY 
CALL CPO°F(SP,PO,VV) 

V/MPO = V M AG ( PO ) 

V/MVV = VMAG(VV) 

D IS DISTANCE from VORTEX AXIS TO POINT OF INTRFST 
n = \/m\/V/VMPO 

PV = EXP ( - ( 4 .*1 INU*7 ) ) 


VMVT IS MAGNITUDE OF THE TANGENTIAL VELOCITY AT THE POINT 
VMVT = GAMMAC*(1. - EX ) / ( ? •*P I *n ) 

VMVTP = VMVT/VMVV 


VT IS THE TANGENTIAL VELOCITY VECTOR AT POINT OE INTREST 
VT < 1 ) =VMVTD*W(1) 

VT (?) =\/M\/TD*VV ( .? ) 

VT ( 3 ) = \/MVTD*VV ( ) 


VMA\/r 

V/MAy/Tpl 

AVT ( 1 ) 
AVT ( P ) 
AVT ( S) 


VMAVT IS MAGNITUDE OF AXIAL VELOCITY 
O7EP0*fx/ ( A . *P I *RH0*UNU*7 ) 

= VMAVT /VMPO 
= \/MA VTD*-PO ( 1 ) 

= VMAVTD->TPO ( ? ) 

= VMAVTD*P0(3) 


AT POINT OE 


INTREST 


VL I NF 1 IS LINF OF SIGHT COMPONENT OF TANGFNTIAL VELOCITY 
VLINFt = -DOT ( VT ,OS) /D! FTL 


VL I NF? IS LINF OF SIGHT COMPONENT OF AXIAL V C LPC I T Y 
VL TNfo = -HOT ( AVT , OS ) /n I FTL 


r VLINF 1° TOTAL LINE OF F I GTH V F LOClTV COMPONENT 

VLINF = \/L I NF 1 + VLINF? 

PFTl |ON 
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*• i n F r rc per am npox 

Cl iRpnt IT I NP per AN ( MCW ) 

r 

r THIS ROUT 1 MP HANDLES THF PLOTING OF THF VELOCITY? 

r 

rnMMAM / APR / TITL.FU?) 

CPMMPN / APR / V(?O,?n,?0) 

COMMON / AP4 / PTDl(pO,3) 

COMMON / APO / PTD(3*3) 

COMMON / A 10 / VV (?0» RO ) * VVV(?0»2C') 

COMMON / All / IFYM(RO), FLDX ( l?i3) » FLDY ( 1 R ) 

COMMON / AIR / DPT0C3) 
n 1 MpMC I ON LAR ( R * 3 ) 

O I MPMC ! ON NPTO (3,3) 

01 MPMC ION ARC(RO) 

FOU I vALFNCF (DTP, NPTO) 

OAT A ( L ARC 1,1 ),T=1,2) / 6H P,6HHI = / 

PATA (l„AR( I ,R) , I =1 ,?) / 6H THP , 6HT A = / 

OATA (L*R( T ,3) ,1=1 ,?) / 6HDTPTAN,6HCF = / 

OAT A I CYM / 17,18,19,20 ,21 , R2 , 23 , R4 , 25 , 33 ,34,38,36,37,38,39, 

1 40,41 , R O , R 1 / 

DATA <FLDX( I , 1 > « 1=1 , 12) / 6HPHI 1N,6H DFGRF , 6HFS ,9*1H / 

DATA ( FLDX (I, R), 1=1,1?) / 6HTHETA ,6HIN DFG , 6HRFFS ,°*1H / 

DATA ( FLDX (I,3),T = 1*1?) / 6HD I ST AN , 6HCF IN , 6HMFTPRS , 9* 1 H / 

DATA CLOY / 6HVEL0C I , 6HTY IN , 6HD I RFCT , 6H I ON 0F,6H LASER, 

* 6H MF , 6HTFRP/P , 6HF0 ,4*1H / 

PAT A ARC /?HA= ,2HB= , RHC= , 2HD= , ?HF= , ?HF= , ?HG= , 2HH = , 2H I = , 2HJ= , RHK= , 
R 2HL= , RHM= , 2HN= , 2HO= , 2«Pr , 2H=0, 2HR= , 2HP= , 2HT=/ 

CALL CAMPAV(9 ) 
r 

r CALL FOR PLOTING OF INPUT 

CALL P INPUT 

r 

r PFAD PLOTING INSTRUCTION 

0 ILINF = 0 INDICATFP NO LINFP BETWEEN PLOTED POINTP 

o TLINE = 1 INDICATES LINES BETWEEN PLOTED POINTS 

C IX = 1 INDICATFP PHI TO RE PLOTED ON HORIZONTAL AXIc 

C IX = 2 INDICATFP THFT A TO RF PLOTED ON HORIZONTAL AXIS 

O IX = 3 INDICATFP RANGE TO RF PLOTED ON HORIZONTAL AXIS 

r T v = 1 INDICATES PHI TO fr PLOTED ON VERTICAL AXIS 

C TY = 2 INDICATFP THFjA TO RF PLOTFQ ON VERTICAL AXIS 

C TY = 3 INDICATFP RANGE TO RF PLOTED ON VFRTICAL AXIS 

IP °F An ( f , fpp ) I L I h'F , IX, IV 
= 00 pORMAT (RJIO) 

I F ( I x .LF. p .OR. IY .LF. O) GO TO 400 
IZ = 6 - IX - IY 
NX = NPTD ( 3 , I X ) 

MY = NPTD ( 3 , I Y ) 

MZ = NPTO ( 3 , I Z ) 

NNX = NX 

IF( ILINF .GT. 0) NNX = -MX 
no R nn K=l, N 7 
VMAy =- i ,c?o 
\/M I M = 1 . PRO 

CO jo (1R, 13, 14), IX 
1 R on TO (IF, 16, 17), IY 

13 GO TO (18, IF, ip), iv 

14 GO TO ( RO ,R1,1R),TY 


B-41 



LMSC-HREC D225936 


1 c '<<9 T T*" ( 6 * 6 1 o ) 


6 1 o F 0 PMAT( 99 H #** FPPOP *** IV = I Y) 



ro Tr )r> 




1 * 

A<;cjrM pr, 
ro To 40 

TO fo 



1 y 

APPTON ?y 

ro TO 40 

To JP 



IP 

A^cirN* PO 
(to TO 40 

TO TP 



lo 

A P c TON? PO 

OO TO 40 

TO IP 



pn 

A PC T on or 

ro to 40 

TO TP 



P\ 

ASP I ON PI 

TO IP 



40 

OC 41 J= ! t 

MV 




OO 4 1 1 = 1 1 

MV 




ro t^ T P ♦ ( « p 8 « P° 

,30, 

31 ) 

.?p 

VV ( 1 , J ) = 

90 TO 9? 

Vt I , J»Y) 



?7 

VV ( I , J ) = 

90 TO op 

v ( i , y , j ) 



PP 

VV ( T , J ) = 

90 TO 9? 

V ( J , I , K ) 




W< I , J) = 
90 TO 9? 

V(K, I , J) 



S" 

vv < r , j ) = 

90 Tn 3 p 

V(J,Y« I ) 



PI 

VV ( I , J ) = 

V ( K . J , I ) 



PP 

tp< vv( r , J) 

• 9 t # V/MA.y ) 

V M AX = 

* 1 

|F< VV( T . J) 

.LT. V M T N ) 

vyin = 


TF(VWAV .LF, VWIN) GO 

TO 

POO 


nr>\/= ( \/v» A y 
L = 0 

- VIYIN)*.?8 




T I A - oov# 1 00 no, 




I A= I T A 




5^ 

L = L+ 1 

I A= I A/ 1 0 





I F ( I A , NF . 0 ) GO TO 80 





VV ( IiJ) 
VV ( I * J ) 


OV = ( T I A/1 O** ( L- 1) ) * 1 0** ( L- 1 ) 


nv = nv#,nooi 
OO Cr Jap, MY 
no Fc i=i, My 

VVV(T «J) = VV( I . J) + OV/*FLOAT( J- 1 ) 

\/MAyv = \/M Ay + OV*FL.OAT ( NV- 1 ) 


PLOT 8 VAPIARLF OF HORIZONTAL avis V^PGUF VFLOCITV 
FALL OU I l<3L ( — 1 , PTD( 1 * IX) , PTD(?, IX) , VMIN,VMAXX, ISYM< 1 ) , FLOX( 1 , IX) 
* FLnv , NNy , PTO 1 ( 1 , T X ) , VV (1,1)) 
no inn J = .?,MY 

1 00 CALL OU IK.3L ( 0 , DUN , OUW ,OUM,DUM, I SYM ( J ) , DUN , HUM , NNX , PTD 1 ( 1 , I X ) « 
ff >VVV ( 1 , J ) ) 

CALL DRINTV(7?,TITLF, 38,1008) 

CALL PR I MTV ( 1 0»LAR( 1 , 17) ,61 6 ,988) 

FALL L A9LV ( PTD 1 ( K , T7) .904,088,-6, 1,6) 

CALL n P I NTV ( A, 4 Hnv =,840, 988) 

CALL L A PL V ( OV ,88o, 988,-6,1,6) 

CALL PPAA'FV 

CALL PPIMTV(79,TTTLF, 1 24, loop) 
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TALI., np I MTV ( 1 0,LAR( l , 17) , 740 , 1008) 

C ALL L APLV ( PTO 1 ( K , 17) ,P?8i 1 OOR, -6, 1 .6) 

CALL PP INTV <18, FLOY < 1 «IX)«440, 16) 

CALI. AFPMTV < 0 . -1 R, 1 R , FLOX ( 1 , 1 Y ) , 1 6 « 6E6) 

call op i oi v < ? , PTn< i , i x ) , pt d< 2 ♦ I x ) , pto< i , i y ) , ptr < ? , i y ) , opto ( i x ) 

T nPTn(!Y), c P,PO,-?,-?,-P,-=;) 

T ARC = 144 

J ARC = 1 6R 
KK= 1 

MV = V/Mflv/R. 

VCON = NV*R 

I F ( v m ,A X .LT. O.) VC ON = VCON — R • 

VCON 1 = VCON 

Pio IF ( Vo ON .LT. VMIM) GO TO 22E 

CALL PP I NT V ( 2 » ARC ( X'K ) » I ABC « 988 ) 

CALL L ARLV ( VCON « JA°C . PR 8 , 3 , 1,3) 
no 7 ?o T = 1 , NX 
I V = ? 

T F (VC ON .AT. VV (1*1)) IV=1 
no 7 no J-Pt NY 
CO T n ( 777, ??R) , IV 
07? I F ( ycON .OT . VV ( I , J ) ) CO' TO 770 
TV = 7 

CO TO 774 

773 t F ( VCON .LT. VV ( I , J ) ) CO TO 220 
I V = 1 

??4 VP = ( VCON - VV ( I * J— 1 ) ) / ( VV ( I * J ) — VV ( I * J— 1 ) ) 

VVP = VP*(PTD1(J , IY)-PTD1 ( J-l , IY) ) + PTOl(J-l.lY) 

PLOTR VELOCITY CONTOUR 

CALL A PLOT V ( 1 * D Tn 1 ( I , ly) ,VVP, 1*1,1, ISYMCKK) , I FRR ) 

??o CONTINUE 
kic= 1 

VCON = VCON - F. 

I ARC = I ARC + R6 
J A BC = J ARC + R6 
CO TO 710 

7?f Kl<r=l 

VCON =VCOM 1 

?oo I F ( VC ON .LT. VM IN) CO TO 200 
OO ??6 J=1 , NY 

I V = 7 

IF (VCON .CT. VV(1»J>) IV = 1 

OO 2?6 I =7, NX 
CO TO (777,778), TV 
777 IF (VCON .CT. VV ( I , J ) ) CO TO 276 
IV = 7 
CO TO 770 

77° I F ( VoON .LT * VV ( T , J ) ) CO TO 2 26 
I V = 1 

??o v D = (VCON - VV ( I — 1 , J ) ) / ( VV ( I , J ) - VV( T — 1 ,J) ) 

VVP = VP*(PTD1 ( I 4 IX) - PTO 1 ( I - 1 , I X ) ) + PTD 1 ( I — 1 , I X ) 

PI_OTR VELOCITY CONTOUR 

CALL A PLOT V ( 1 , VVP, PTD 1 ( J , I Y ) , 1 , 1 » 1 , I RYV ( XX ) , I ERO ) 

776 CONTINUE 

KK = ICY + ] 
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\/com = \/roM -o. 

CO T<o POO 
On C ONT T Nt JF 
GO T n in 
on CALL CL FAN 
pcr Tt JON 
C MD 
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l p FTr P I mpi it nrcK 

eilROoUTlNF PJMPUT 

rOMMnM / API / nDH) , PO(3) 

COMMON / AO? / TTTLFMP) 

COMm^m / A 06 / GAMWAC, HR, t.INU . Z 

COMMON / AOP / PH I 1 » PH I ? t NPH I * T HE T A 1 « THETA;?, NTHET A , 0IST1, 
«. Olcr?, K'DI^T 
COM^pm / A ! A / D7FP0* PHO 

COMMON / APO / PFI, R ♦ FLAMOA, PLFNGT, ICHOIC, WEIGH, IOUT, 

* IPLPT, I PLOT 1 * IVMAX 

call FCr>UTV(P ) 

wQJ T=- ( 1 6,6 PC ) 

uiP J tp (16.6PP) TITLE 

6PC FORMAT C 1H1 1 PA 6) 

WPITp <16,610 GAMMAC, WR» UNU , Z, OP, PO, PH I 1 , PH IP, NPH I, 

1 THFTA1, TH^TAP, NTHFTA, DIST1, OISTP, NOIST, PSI, R, FL AMO A , 

P PLFMGT, OZFRO, RHO, ICHOIC, NWEIGH* I OUT, I PLOT, I PLOT 1, IVMAX 
6JP FORM/nT ( //* 1 X, 1 6HNAMFL I FT (INPUT),// 


1 

1 PH GAMMAC 

= 

,fj 6 .B, 

3 X* 

1 1 HWR 

= . F 1 6 • 

8 . 3 X « 

p 

1 1 HUNt l 

= 

,F 1 6 .B. 

3 x. 

1 1 HZ 

=*F 16 . 

8 . / 

3 

1 2 H OP( 1 ) 

= 

, F 1 6 • R , 

3 X * 

1 1 HOP < ?) 

= .F! 6 . 

8 . 3 X. 

A 

1 i H OP( 0 ) 


, F l 6 .R, 

/ 





1 PH PO< 1 ) 

- 

, f 1 6 * B , 

3 X. 

1 lHPOt P) 

“.FI 6 • 

8 . 3 X% 

6 

1 1 MPO ( 7 ) 


, F 1 6 • B, 

/ 




7 

1 PH PHI 1 

= 

, F 1 6 • B , 

3 Xi 

1 1 HPH I P 

= . F 16 . 

8 . 3 X « 

P 

1 1 UMPHI 

SS 

, IF,/ 





o 

1 PH t HF T A ) 

= 

, F 1 6 • B , 

3 X . 

1 1 hthftap 

= . F 1 6 • 

8 . 3 X . 

A 

1 1 HNTHFTA 

= 

, IB./ 





R 

1 ?h niFTi 

= 

, F 1 6 .B, 

3 X . 

1 1 HO I FTP 

=.F 16 « 

8 . 3 X ♦ 

r. 

1 1 h mo I ft 

= 

.IF./ 





n 

1 ?H PFI 

= 

, F 1 6 • B , 

3 X % 

1 1 HR 

= . F 1 6 • 

8 . 3 X . 

f 

1 1 H^LAMDA 


. F 1 6 • R , 

3 X. 

1 1 HPLFNGT 

= . F 1 6 • 

Bi/ 

F 

1 PH 07^00 

— 

, F 1 6 • B , 

3 X . 

1 1 HRHO 

=iF 16 . 

8 . 3 X . 

r, 

1 1 m I CHOI C 

= 

, IB, 1 IX 

. 

1 ihnwfigh 

= . IP./ 


H 

1 PH I OUT 

= 

, IB. 1 IX 

. 

1 1 H I PLOT 

= . I 8 . 1 

IX. 

I 

1 1 MT PLOT 1 

5S 

, IB. 1 IX 

. 

1 1 H IVMAX 

= ♦ IP) 



PFTt tPN 
FNO 
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tTpFtC C.POcc n^CK 

SUBROUTINE C POSS ( VA , VR» VO ) CROon 

r FuRROUTINF TO PFRFORM CROSS PROOUCT CROno. 

VFCTOR VA is CROSSED INTO VECTOR VR TO FORM VFCTOR VC CROnn 

n i mfms i on va(s) ,vr(S) «vC(R) r.Rono. 

VC ( 1 )sVA(?) * VR ( 3 ) — VR ( P ) *VA ( 3 ) CROPP' 

VC ( 2 ) = V A ( 3 ) *VR ( 1 ) — VR ( 3 )*VA ( 1 ) CROnn, 

VC ( *» > = VA { 1 ) *vp ( P )-VBC 1 ) *VA ( P ) CRQRP' 

R^Tl lr*>N CRO nr " 

fmd rnorr- 1 

STRFTC POT dfck 

F! INCTJON OPT (\./ A, VP) nOTPP 

c FUNCTION TO PERFORM POT PROOUCT DOTop 

O I MFMS I ON VA ( 3 VR < 3 ) OOTOO- 

DOT =\/A ( 1 ) *VR( 1 )+VA(P) *VR ( 2 ) + VA ( 3 ) *VB ( 3) DO TOP 1 

C VFCTOR VA IS DOTTED INTO VECTOR VB TO VORM THE SCALAR DOT PRODUCT DOT 0 r 

RETURN DOTPPt 

FMD DOT pp' 

TIPFTC VWAr. DECK 

FUNCTION VMAG(VA) \/ a* A p P 

F FUNCTION TO DFTFRMINt Victor MACNITUPF \/MAop. 

D I MFstQ I ON \/A(R) VMAPP 

V MAc = SORT ( VA ( 1 )*\/A ( t )+VA < 2)*VA < ,?)+VA( 3 )*VA < 3) ) VMA.no 

DFTURN \/ m A P p ‘ 

FMD \/MA on 
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"■•nAK 

FyAMPLF FOP VOLUMF FC.ANNINP PROGRAM 
I N^t IT 

GAMMAC = 720, « 

t.iP - 77 , 7 , 

I IN.ll .1 = ,AO|A, 

7 = 1 nop, , 

r>7rpn = ernonn , , 

Rho = 1 ,10?, 

On — o,, n,, F o , , 

PO — 1 . , n . , n,, 

PHI l = 7. , 

PHI 7 = -7. , 

NPHT = 11, 

TH C T A 1 = 4 2.* 

THFTA2 = 42., 

mth^ta = H, 

0 T FT 1 = FP., 

O I c T2 = 71 . , 

NO I FT = 11, 

F L A m n A - , nooning, 

PFT = 1,, 

I CHOIC = 1 , 

R = ,07 F 2 , 

NHF1GH =1?, 

1 Ol IT = 7 , 

I PLOT = 1 , 

I PLOT 1 = 1 , 

I v/MAV = 1 , 

<r- 

1 1 2 

1 7 2 
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Appendix C 

LASER VELOCIMETER LINE SCAN SIMULATION PROGRAM 


C - 1 c° 
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Appendix C 


C.l DESCRIPTION 

Program simulates one vortex and scans it in a straight line with one 
or two laser velocimeter s. A sequence of vortex positions, scan positions 
and laser parameters is automatically run for comparisons. Outputs are 

plots of velocity versus position. 

C.2 CAPABILITIES 

1. Vortex can be in any position. 

2. The tangential and axial velocity profile of the vortex can be 
input with different values. 

3. Constant winds can be added. 

4. The vortex positions, scan positions and laser parameters 
of the automatic sequence can be altered. 

5. For two lasers the separations may be altered. 

C.3 RESTRICTIONS 

1. Infinite S/N is assumed. 

2. Aerosol is assumed to be constant throughout the vortex. 

3. The vortex is assumed to be uniform along its length. 

4. The vortex axis is assumed to be a straight line. 

5. No change in vortex position or strength is allowed with change in 
time. 

c-i CP) 

LOCKHEED - HUNTSVILLE RESEARCH & ENGINEERING CENTER 



LMSC-HREC D225936 


6. Scanning is restricted, for one laser scanning is either in line 
with the laser (range scan) or perpendicular to that line (lateral 
scan). For two lasers, the scanning is the same as if there were 
one laser directly in between the two and the single laser con- 
trolled the scanning. 

C.4 POLAR COORDINATE SYSTEM OF LASER 


The laser's focal point orientation is represented in a polar coordinated 
system in the XZ plane of the rectangular coordinate system. 



The transformation from the polar coordinate system to the laser's rectangular 
coordinate system is 

X = r cos 0 
s s s 


Z = r sin 0 
s s s 

C. 5 SCANNING OF LASER 

With one laser the scanning is done with r constant (LATERAL SCAN) 

s 

or with 0 g constant (RANGE SCAN). With two lasers the scanning is treated as 
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in the one laser case but using a reference point directly between the two lasers. 
A sample of these is shown in Fig. C-l. In this figure the vortex is shown in its 
normal position with its axis parallel to the Y axis of the rectangular coordinate 
system. 

C.6 COMPUTATION OF TANGENTIAL VELOCITY COMPONENT OF VORTEX 
IN DIRECTION OF LASER 


This is computed the same way as shown in Appendix B. 


C.7 COMPUTATION OF AXIAL VELOCITY COMPONENT OF VORTEX IN 
DIRECTION OF LASER 


This is computed the same way as shown in Appendix B. 


C.8 COMPUTATION OF VELOCITY INTERPRETED BY A COAXIAL FOCUSED 
LASER DOPPLER SYSTEM 


This is computed the same way as shown in Appendix B. 


C.9 COMPUTATION OF VELOCITY INTERPRETED BY A COAXIAL PULSED 
LASER DOPPLER SYSTEM 


This is computed the same way as shown in Appendix B. 


C.10 COMPUTATION OF VELOCITY INTERPRETED BY A COAXIAL FOCUSED 
CONTINUOUS WAVE LASER OR A COAXIAL PULSED LASER USING THE 
HIGHEST ACTIVATED FILTER TECHNIQUE 

In addition to the centroid method described in Appendix B this program 
can use the technique of highest activated filter. The same weighting over the 
illuminated volume is used but in this case the weighting value for each sampled 
location is added to a slot representing a filter sensitive to a particular velocity 
range. After the illuminated volume has been sampled, the line of sight velocity 
representing the volume is chosen to be the mid point of the filter, representing 
the highest velocities, which has a predetermined fraction of the signal or more. 
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C.ll INPUT FOR LASER VELOCIMETER LINE SCAN SIMULATION 
PROGRAM 


The first input is a title card (columns 1 through 72). Then three name- 
lists comprise the remainder of the input. Their names are VORTEX, SYSTEM, 
and SCAN. They must be input in that order. All input variables have preset 
values. The values describing the vortex velocity flow field are representative 
of that produced by a Boeing 747 while landing. 


Namelist VORTEX has the following variables. 


Name 

Preset Value 

Units 

Description 

GAMMAC 

729.0 

n-m-sec/kg 

Circulation of vortex from 
generating aircraft. 

W 8 

77.3 

m/sec 

Freestream velocity 

UNU 

.0914 

m /sec 

Eddy viscosity 

Z 

1000.0 

m 

Down stream distance of 
vortex from generating 
aircraft . 

DZERO 

50000.0 

N 

Profile drag of the aerofoil 
(used for axial velocity. If 
no axial velocity wanted set 
DZERO to zero). 

RHO 

1.192 

kg/m 3 

Density of atmosphere 

VROTH 

0 . 

deg 

Vortex rotation about 
reference sweep line 

VROTHC 

10. 

deg 

Change in vortex rotation 
about reference sweep line. 

V ROT V 

0 . 

deg 

Vortex rotation about site 
line. 

VROTVC 

10. 

deg 

Change in vortex rotation 
about site line 

NRANGE 

2.0 

none 

Indicates nominal vortex 
range 

VRANGE(8) 

50.0 

m 

Vortex ranges 


100.0 

300.0 

1000.0 

C -5 
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Name 

Preset Value 

3000.0 

5000.0 
-1.0 
-1.0 

Units 

Description 


VMAX 

34.0 

m/ sec 

Upper velocity bound 
scale 

on plotting 

VMIN 

0 

m/ sec 

Lower velocity bound 
scale 

on plotting 


Namelist SYSTEM has the following variables. 


IUNITS 

1.0 

none 

FLAG 

IUNITS = 1 indicates one laser 
IUNITS -Z indicates two lasers 

ICHOIC 

1.0 

none 

FLAG 

ICHOIC = 1 indicates continuous 
wave focused laser 
ICHOIC -Z indicates pulsed 
laser 

PSI 

1.0 

— 

Controls magnitude of weighting 
factor. Leave set to 1. 

FLAMDA 

.0000106 

m 

Wave length of laser light 

NR 

2.0 

none 

Indicates nominal radius of 
telescope aperture 

R(8) 

0.1 

0.2 

0.3 

0.5 

1.0 

2.0 

4.0 

6.0 

m 

Radius of telescope aperture 

NPLENG 

2.0 

none 

Indicates nominal pulse length 


C-6 


LOCKHEED - HUNTSVILLE RESEARCH & ENGINEERING CENTER 



LMSC-HREC D225936 


Name 

PLENGT(8) 

WEIGHL 

IABS 

SLASER 

DLASER 

EIjASER 

Namelist 

IDIREC 

SCAWL 

SCANR 

NSWEEP 

SWEEP(5) 


Preset Value 

Units 

Description 

30.0 

90.0 

150.0 

300.0 

900.0 

m 

Pulse length 

-1.0 



-1.0 



-1.0 




0.5 

m 

Distance between sample points used 
in weighing 

1.0 

none 

Flag 

IABS = 1 indicates no distinction be- 
tween negative and positive velocities 
are detected 

IABS -2 indicates positive and nega- 
tive velocities are detected. 

0.0 

m 

Distance of half laser* s separation 
along site line 

0.0 

m 

Distance of laser from sweep plane 

100.0 

m 

Distance of lasers from site line in 
sweep plane 


SCAN has the following variables. 


1.0 

none 

Flag 

IDIREC = 1 indicates lateral scan 
IDIREC ~Z indicates range scan 

30.0 

m 

Scan length 

0.5 

m 

Scan resolution 

3.0 

none 

Indicates nominal sweep position 

-15.0 

-7.0 

0 

7.0 

15.0 

m 

Sweep positions 
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Description 

Flag 

IVMAX = 1 indicates the maximum 
minus the minimum line-of-sight 
velocities observed in the illuminated 
volume will be plotted as X. 

IVMAX = 2 indicates the maximum of 
the absolute line-of-sight velocities 
observed in the illuminated volume 
will be plotted as X. 

Flag 

IWEIGH = 1 indicates that the highest 
activated filter will be plotted as W 
IWEIGH = 2 indicates that the centroid 
of the weighted velocities will be 
plotted as W 


WFILT 

0.94 

megahertz 

Band width of filters 

DETLEV 

0.1 

none 

Detection level in fraction of signal 
for filters 

OFSETX(2) 

0 

m 

Offset of laser focus point along site 


0 


line due to misalignment 

OFSETY(2) 

0 

m 

Offset of laser focus point from scan 


0 


plane due to misalignment 

OFSETZ (2) 

0 

m 

Offset of laser focus point along sweep 


0 


line due to misalignment 

WIND(3) 

0 

m/ sec 

Constant WIND vector 


0 

0 


Name 


Preset Value Units 


IVMAX 


none 


IWEIGH 


none 


After case is completed, control returns to start of program. 
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C.1Z OUTPUT FOR LASER VELOCIMETER LINE SCAN SIMULATION 
PROGRAM 


The input is outputted on both the printer and the plotter. See Fig. C-2. 

The remainder of the output consists of five plots each containing three 
graphs with four curves. Examples of these are shown in Figs. C-3 through C-7. 
The four curves that appear on each graph are identified by four letters W, 

X, T and L. They are all velocities plotted as a function of position. 

W is velocity of the highest activated filter or as indicated by the 
input IWEIGH, or in other cases W is the centroid of the weighted 
velocities. 

X is the maximum of the absolute line-of-sight velocities observed 
in the illuminated volume, or as indicated by the input IVMAX. In 
other cases X is the maximum minus the minimum line-of-sight 
velocities observed in the illuminated volume. 

T is the norm of the total velocity vector at the focal point. 

Lis the line of sight velocity at the focal point. 

The first plot Fig. C-3 contains three graphs in which each represents 
a different sweep position. For the lateral scan, the sweep position represents 
a difference in range between the vortex axis and the sweep. For the range 
scan, the sweep position represents the distance that the sweep has missed the 
axis . 


The second plot Fig* C-4 contains three graphs in which each one represents 
a different vortex range. The center graph is identical to the center graph of 
Fig. C-3. This graph is called the nominal. It appears in the center of each plot 
and only one parameter for each of the plots will vary with the graphs. 

The third plot, Fig . C - 5, contains three graphs in which the vortex is ro- 
tated in different amounts about the reference sweep line. 

The fourth plot, Fig.C-6, contains three graphs in which the vortex is 
rotated in different amounts about the site line. 
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The fifth plot, Fig. C-7, contains three graphs in which the telescope 
aperture radius is different for each. Or in the case where a pulsed laser is 
used, the three graphs represent different pulse lengths. 

For two lasers W, X and L are the resultants from the components seen 
by the two lasers. However, since X has no sign, it often badly represents the 
velocity. T is still the norm of the total velocity vector at the focal point. In 
the event of misalignment the focal point that T represents is that of the second 
laser. 
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CA#C MO 24 1 LASER RA M6€ MO LATERAL SCAM FILTERS 

NAMELIST (VORTEX) 

6*IWAC » 0.T29OOOOOE 09 Ml * O.7T30OOOOE 02 LMu * 0 .914000O0E-0I * * O.IOOOOOOOE 04 

OzCRO « O. OOOOOOOOE- 36 RhO a 0. It 920000 £ 01 VROTM » Q . OOOOOOOOE- 36 VROTMC ■ 0,90000000 C 02 

VROTV * O. OOOOOOOOE- 96 VROTVC * 0.90000000E 02 K*ANGE * 9 

VRAMCC(t) « 0 . 30000000 E 02 VRAM6C (2) * O.IOOOOOOOE 09 VRANGC <3) s 0.90000000E 09 VRANCC(4) a O.IOOOOOOOE 04 

VRAMtE<3) a 0.90000000C 04 VRAM6C 1 6) a 0.5000000QE 04 VRANGC(7) a -O.IOOOOOOOE 01 MtANCC(S) « -O.IOOOOOOOE 01 

VHAX a 0.94000000C 02 V*IN s O. OOOOOOOOE -3© 


NAMELIST (SYSTEM) 


IUNIT6 

a 

t 












lOtOf c 

a 

1 


PSI 

a 

O.IOOOOOOOE 

01 

FLAMDA 

3 

0. 1 0 600000 E- 04 

m 

a 

r 

R(l) 

a 

O.IOOOOOOOE 

00 

R(2) 

B 

0.20000000E 

OQ 

R (9) 

= 

0.90000000E 00 

R (4) 

a 

0 . 90000000E 00 

R(9) 

a 

O.IOOOOOOOE 

01 

R<6> 

a 

0.20000000E 

01 

R (T) 

s 

0.40000000E 01 

R(8) 

a 

0.600CIOOOOE 01 

NPLCNS 

a 

' 2 












plcn*t(1 > 

a 

O . 30000000C 

02 

PLEMCT<2) 

3 

0.90000000E 

02 

PLEWC-T (51 

a 

0.15000000E 05 

PLENST (4) 

a 

0.30000000E 09 

PLEM*T(9> 

a 

0.90000000C 

09 

PLEN9T (6) 

= 

-O.IOOOOOOOE 

01 

PLEWGT(T) 

= 

-O.IOOOOOOOE 01 

PLENOT (9 ) 

a 

-O.IOOOOOOOE 01 

vcichl 

a 

0 . 90000000E 

00 

IABS 

3 

1 


9 LA8ER 

= 

0. OOOOOOOO E- 96 

olascr 

a 

0. OOOOOOOOE- 96 

CLASER 

a 

O.IOOOOOOOE 

09 












NAMELIST (SCAN) 


IOIRCC 

a 

1 

9CANL 

3 

0.30000000E 02 

SCAN* 

= 

0.5 OOOOOOOE 00 

NShCEP 

a 

3 

»wCCPP cl ) 

a 

-0 . TOOQOOOQC 01 

SWEEPP (2) 

3 

-0.20000000E 01 

SWEEPP (3) 

- 

0 . OOOOOOOOE- 36 

SWEEP P (4) 

a 

O.TOOOOOOOE 01 

•MCCPPCS) 

a 

O.ISOOOOOOC 02 










ivwax 

a 

c 

IV£l«tt 

S 

1 

fc#*ILT 

3 

0 .9400 OOOO E 00 

OCT LEV 

3 

O.IOOOOOOOE OO 

OFSCTX (1) 

s 

0 . 0000 000 OE -3 0 

OFSETX (2) 

= 

0*0000000 OE -38 

OFSETY (1) 

= 

0. OOOOOOOOE -38 

OFSETY (2) 

= 

0. OOOOOOOOE -38 

orscTZd) 

3 

0. OOOOOOOOE -30 

OF SC TZ (2) 

3 

0.0000000 OE -3 8 







WIND (1) 

a 

0 . OOOOOOOOE —9 6 

WIND (2) 

3 

0. 0000000 OE -38 

WIND (3) 

3 

0. OOOOOOOOE -38 





Fig. C-2 - Output for Laser Velocity Line Scan Simulation Program 
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V 

K 

L 

0 

c 

1 

? 

V 


8 

M 

/ 

d 

c 

c 


data point position with respect to site line in heters 



i 


Fig. C-3 - Example Output 


C- 12 


LOCKHEED - HUNTSVILLE RESEARCH & ENGINEERING CENTER 


A«® % * 


LMSC-HREC D225936 




s 

* 

V 

s 

M 

b 

s 

c 

S 

C 



Fig. C-4 - Example Output 
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OAT A POINT POSITION WITH RESPECT TO SITE LINE IN HE TER 8 



Fig. C-5 - Example Output 
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Fig. C-6 - Example Output 
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E 900 LATERAL SCAN FIllOl 
TCLCSCOSC RADIUS IS 3.000 MO 


Norm. 
t * Tota" 


X » NAX(ASS(Vl>> U " J 
L * IIIC OF SICNT W * HI«H ACffVATCO FflTC* 


| BB ■■■■ ■■■■ ■■■■ ■■■■ BBBB K'<iBBtMW3 ■■BT/iBBfttMBBBB BBBB BBI 

iiiniiiiiiinniiiii/iiiiiiiiiiiv JiiatAVinainaa 

|BBBBBBBBBBBBBBBBBBBBr'BBaBBBMBBBBrJBBBBBBMBBBBBBBBI 

|BiBBBBBBBBBBBBBBBflr.:BBBBBBBBBtaBBBBBBBBBBKJBBBBBBI 

lii5 B B BBBBBBBBBBBBI &^ BBBBBBBBB t» B VABBBBBB BB K^«MBBI 

5i9 B S BBBBBBBBB *^^ BBBBBBBBBB ^ BB ^ BBBBBBBBBB;fr ^ 1BBI 

■■■■■■■■M-^'flBBBBBBBBBBBBBB'lBIBBBBBBflBHBBBB^I 


■■■■■■■I 

BBBBBBBI 

■■■■■■■I 

■■■■■■■I 

BBBBBBBI 


iwysjatiimmmmMHiMifMMiiHm— imi WHH 

U«*»*«BBBBBBBBBBBBBBBBBBBBBBBUMBBBBBBBBBBBBBBBBBBBBIH 
BBBBBBBBBBBBBBBBBBBBBBBBBBBBBIUIBBBBBBBBBBBBBBBBBBBBBBB 
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C . 1 3 EXAMPLE RUN FOR LASER VELOCIMETER LINE SCAN 
SIMULATION PROGRAM 


As an example, the axial component of the vortex 1000 meters downrange 
from a Boeing 747, will be examined* For the nominal case, let the vortex be 
parallel to the Y axis at a range of 300 meters* For the laser system, one con- 
tinuous wave focused, CO^ laser will be simulated. Assume that the velocity sign 
cannot be determined. Let the nominal telescope aperture radius be 0*2 meters. 

A lateral scan of 30 meters with 0.5 meter increment is used. Let the 
nominal sweep portion go through the center of the vortex. On the W curve 
the filter response is observed, and on the X curve the maximum velocity in 
the focal volume is observed. 


The first input must be the title card. CASE NO. 24 1 LASER RANGE 
300 LATERAL SCAN FILTERS. 


The next input is the namelist VORTEX for its variables. 


GAMMAC 

leave 

preset 

to 

729.0 

W8 

leave 

preset 

to 

77.3 

UNU 

leave preset 

to 

0.0914 

Z 

leave 

preset 

to 

1000 

DZERO 

set to 

0 



RHO 

leave 

preset 

to 

1.192 

VROTH 

leave 

preset 

to 

0 

VROTHC 

set to 

30 



VROTV 

leave 

preset 

to 

0 

VROTVC 

set to 

30 



N RANGE 

set to 

3 



VRANGE 

leave 

as preset 



Set for Boeing 747 


Since no axial velocity 
is desired 

Typical low altitude density 


To look at larger off- 
nominal vortex rotations 


To look at larger off- 
nominal vortex rotations 

300 meters is the third 
preset value for VRANGE 
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VMAX 

VMIN 


The next 

IUNITS 

ICHOIC 

PSI 

FLAMDA 

NR 

R 

NPLENG 

PLENGT 

WEIGHL 

IABS 

SLASER 

DLASER 

EEASER 

The last 

IDIREC 

SCANL 

SCANR 

NSWEEP 


leave preset to 34 


leave preset to 0 


Good value for selected 
vortex 

Since all velocity will be 
positive, zero is a good 
lower limit 


namelist to be inputted is SYSTEM, For its variables: 


leave preset to 1 
leave preset to 1 
leave preset to 1 
leave preset to 0,0000106 

leave preset to 2 

leave as preset 
leave as preset 

leave as preset 
leave preset to 0.5 

leave preset to 1 
leave as preset 1 
leave as preset \ 
leave as preset | 


This is the wavelength of 
light emitted from a CO 
laser 

0.2 meters is the second 
preset value for R 


This is of no significance 
since pulses are not being 
used 


This is a reasonable interval 
to sample the illuminated 
volume 


No significance, since only 
one laser is needed 


namelist to be inputted is SCAN. For its variables: 


leave preset to 1 
leave preset to 30.0 
leave preset to 0.5 
leave preset to 3 


0 is the third preset value 
for SWEEP 
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SWEEP 

set SWEEP ( 1 ) to -7 
and SWEEP(2) to -2 



IVMAX 

leave 

as preset to 2 



IWEIGH 

leave 

as preset to 1 



WFILT 

leave 

as preset to 0.94 

Filters are to have an 





effective 5 m/sec velocity 
width 

DETLEV 

leave 

preset to 0.1 


This is a reasonable signal 
to-noise fraction for the 





filter 

ofsetx 

leave 

preset to 0., 0., ] 



OFSETY 

leave 

preset to 0., 0., 

! 

Perfect alignment is 

OFSETZ 

leave 

preset to 0., 0., 1 


as sumed 

WIND 

leave 

preset to 0., 0., 0., 

No wind 


The output from this sample run is shown in Figs. C-2 through C-7. 
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C.14 LASER VELOCIMETER LINE SCAN SIMULATION PROGRAM 
COMMON VARIABLES 


Name 

Units 

Location 

Description 

AVT(3) 

m/sec 

COMMON /AO 5/ 

Axial velocity vector 
at illumination point — 
(rectangular coordinate 
system). 

D 

m/ sec 

common/a 15/ 

D is the distance from 
illuminated point of 
interest to axis 
of vortex. 

DD( 100, 2) 

m 

COMMON/A03/ 

Contains distance from 
vortex axis to focus 
point. 

DD( , 1) is for other than 
nominal case. 

DD( , 2) is for nominal 
case. 

DETLEV 

— 

common/ A22/ 

Detection level of filters 
in fraction of total return 
signal. 

DISPLA(3) 

m 

common/ao 1 / 

For use with two lasers. 

Is dispacement vector for 
laser from reference point. 

DLASER 

m 

COMMON/ A2 l/ 

Determines laser separa- 
tion for two lasers, DLASER 
is distance of lasers from 
sweep plane. 

DZERO 

n 

COMMON/A 14/ 

DZERO is the profile drag 
of airfoil from generating 
aircraft. 

ELASER 

m 

COMMON/A2 l/ 

Determine laser separa- 
tion for two lasers. 
ELASER is distance of 
lasers from site line in 
sweep plane. 

FLAMDA 

m 

COMMON /A2l/ 

FLAMDA = wavelength 
of light. 

GAMMAC 

(N -m- s ec)/kg 

COMMON/A06/ 

Initial circulation 


about vortex. 
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Name 

Units 

Location 

Description 

IABS 


COMMON/ A2 1/ 

A flag where IABS = 1 
indicates that the laser 
cannot detect sign of line - 
of-sight velocity. IABS = 

2 indicates that laser can 
detect sign of the line-oL 
sight velocity. 

ICHOIC 


COMMON/A2 1/ 

A flag where ICHOIC = 1 
indicates a CW focused 
laser and ICHOIC = 2 
indicates a pulsed laser. 

IDIREC 


COMMON/A22/ 

A flag where IDIREC = 1 
indicates lateral scan and 
IDIREC = 2 indicates range 
scan. 

IUNITS 


COMMON/A2 l/ 

A flag where IUNITS = 1 
indicates one laser and 
IUNITS = 2 indicates two 
lasers. 

IVMAX 


common/ A22/ 

A flag where IUMAX = 1 
indicates the maximum 
minus the minimum line- 
of-sight velocity is to be 
plotted as X. IVMAX =2 
indicates the maximum 
absolute line-of-sight velocity 
is to be plotted as X. 

I WEIGH 


COMMON/A22/ 

A flag where IWEIGH = 1 
indicates that the velocity for 
the highest activated filter 
will be plotted for W. IWEIGH 
= 2 indicates that the velocity 
from the centroid of the 
illuminated volume will be 
plotted for W. 

JJ 

— 

COMMON/ A08/ 

Indicates which plot is being 
made. 

NOM 

— 

COMMON/A08/ 

NOM = 1 indicates non- 


nominal scan being made. 
NOM =2 indicates nominal 
scan being made. 
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Name 

Units 

Location 

Description 

NPLENG 

~ 

COMMON/A2 l/ 

Indicates nominal pulse 
length. 

NR 

— 

COMMON /A2l/ 

Indicates nominal radius 
of telescope aperture. 

NRANGE 

— 

COMMON /A20/ 

Indicates nominal vortex 
range. 

NSWEEP 

— 

C OMMON / A2 2 / 

Indicates nominal sweep 
position. 

NTHETA 

— 

COMMON/AO 8/ 

Number of points along 
sweep. 

OFSETX(2) 

m 

COMMON/A30/ 

Offset of laser focus point 
along site line due to mis- 
alignment. OFSETX(l) is 
for first laser and OFSETX 
(2) is for second laser. 

OFSET Y(2) 

m 

COMMON/A30/ 

Offset of laser focus point 
from scan plane due to mis- 
alignment. OFSETY(l) is 
for first laser and OFSETY 
(2) is for second laser. 

OFSETZ(2) 

m 

COMMON/A30/ 

Offset of laser focus point 
along sweep line due to mis- 
alignment OFSETZ(l) is for 
first laser and OFSETZ(2) is 
for second laser. 

OP(3) 

m 

COMMON /AOl/ 

Position of a point on the 
vortex axis (rectangular co- 
ordinate system). 

OS(3) 

m 

COMMON/A05/ 

Vector to indicate illumina- 
tion point relative to laser 
(rectangular coordinate 
system). 

OSS(3) 

m 

COMMON/AO 5/ 

Vector position of focus 
point for IASEL. 

PI 

— 

COMMON/A07/ 

PI is 7T = 3.1415926. 

PLENGTH(8) 

m 

COMMON/A2 l/ 

Pulse length. 
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Name 

Units 

Location 

Description 

PQ(3) 

— 

common/ao 1/ 

Directional vector for 
vortex (rectangular co- 
ordinate system). 

PSI 

— 

COMMON /A2l/ 

Is incident wave amplitude. 

R(8) 

m 

COMMON/ A2 1/ 

Radius of telescope aper- 
ture. 

RHO 

kg/m 3 

COMMON/A 14/ 

RHO is the density of the 
atmosphere. 

R2D 

dig/ rad 

COMMON/ A07/ 

Conversion factor for 
converting from radians 
to degrees. 

SCANL 

m 

COMMON/A22/ 

Scan length. 

SCANP(IOO) 

m 

COMMON/A03/ 

Is distance from focus 
point to site line. Also 
is variable that runs on 
horizontal axis of plots. 

SCANR 

m 

COMMON/A22/ 

Scan resolution. 

SP(3) 

m 

COMMON/AO 5/ 

Vector from illumination 
point (OS) to point OP on 
vortex (rectangular co- 
ordinate system). 

SLASER 

m 

COMMON/A21/ 

Determines laser separation 
for two lasers. SLASER is 
distance of lasers along 
site line from reference 
point. 

SMAX 

m 

COMMON/AO 8/ 

Maximum position along 
sweep. 

SMIN 

m 

COMMON/AO 8/ 

Minimum position along 
sweep. 

SWEEPP(5) 

m 

COMMON /A2 2/ 

Sweep position. 

TITLE( 12) 

aphanumeric 

data 

COMMON/ A02/ 

Contains INPUT descrip- 
tion of run. 

UNU 

2 / 

m /sec 

COMMON/AO 6 

Eddy viscosity. 
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Name 

Units 

Location 

Description 

VAL 

m 

COMMON/AO 8/ 

Distance of sweep from 
reference sweep line (for 
lateral scan). 

Distance of sweep from 
site line (for range scan). 

VLINI 

m/ sec 

COMMON/A 15/ 

Line -of-sight component 
of the tangential velocity. 

VL.IN2 

m/ sec 

COMMON/A 15/ 

Line-of-sight component 
of the axial velocity. 

VMAVT 

m/ sec 

COMMON/A 15/ 

VMAVT is the norm of the 
tangential velocity. 

YMAX 

m/ sec 

COMMON/A08/ 

Upper velocity bound on 
plotting scale. 

VMOSS 

m 

COMMON/AO 5/ 

Norm of vector (OSS). 

VMIN 

m/ sec 

COMMON /AO 8/ 

Lower velocity bound on 
plotting scale. 

VMVT 

m/ sec 

COMMON /A 15/ 

VMVT is the norm of the 
tangential velocity. 

VRANGE(8) 

m 

COMMON/ A20/ 

Vortex ranges. 

VROTH 

deg 

COMMON/A20/ 

Vortex rotation about 
reference sweep line. 

VROTHC 

deg 

COMMON/A20/ 

Change in vortex rotation 
about reference sweep line. 

YROTV 

deg 

COMMON/A20/ 

Vortex rotation about site 
line. 

VROTVC 

deg 

COMMON/A20/ 

Change in vortex rotation 
about site line. 

VT(3) 

m/ sec 

COMMON/AO 5/ 

Tangetial velocity vector 
at illumination point 
(rectangular coordinate 
system). 

VTOL 

m/ sec 

COMMON/A 15/ 

Norm of velocity vector at 
focus point. 


C 
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Name 

VV(3) 

VI(100, 4, 


WEIGHL 

WFILT 

WIND(3) 

W8 

Z 


Units Location 

- COMMON/AO 5/ 

2) m/ sec COMMON/ A03/ 


m COMMON/ A2 1/ 

megahertz COMMON/ AZZ/ 

m/ sec COMMON/A31/ 

m/ sec COMMON/ A06/ 

m COMMON/A06/ 


Description 

Vector pointing in direction 
of tangential flow at illumina- 
tion point (rectangular co- 
ordinate system). 

Velocities to be plotted: 

VI( ,1, ) is the weighted 

average or filter output. 

VI( ,2, ) is VMAX-VMIN 

or MAX(ABS(VL)). 

VI( ,3, ) is norm of total 

velocity vector. 

VI( ,4, ) is line of sight 

velocity at focus point. 

VI( , , 1 ) is for other than 

nominal case. 

VI( , , 2 ) is for nominal 

case. 

Distance between sample 
points in weighing. 

Band width of filters. 

Constant wind. 

Freestream velocity. 

Vortex distance from 
generating aircraft. 
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C.15 LASER VELOCIMETER LINE SCAN SIMULATION PROGRAM LISTING 


*CFtijd lb? 

G Acc | GN 
T ATFMP 

^yfci its 
MPJop l _vsr 

CTPFTC m A T m 
C 


TAper, 

cycL ra ? 

OPnro , 77777, 1 , O' IMP 
I BJOR 

DFCK 


C 

c 


THIS PROGRAM is FOP DETERMINING thf usefulness OF VARIOUS 2D 
L A S F p popPLsp SVFTFMC pop USE IN SCANlNG VORTTCFS 



r 


r 

r 

r 

r 

r 

r 

c 

c 

c 

c 

r 

r 

r 

r 

r 

r 


MAVFL I ST ( VOPTFV ) 


G .A M M A f 
»*>R 

UNU 

7 

H7FPP 

PHO 

VPOTH 

VPOTHC 

\/ROTw 

VROTn/C 

MRANpcr 

VpAwrr ( S ) 
W MAy 
\/M I M 


CIRCULATION 
FRppcjRFAM VFLOCITV 
EDPY VISCOSITY 


VOPTFX AVTAL COQPDINATF FROM SOURCE 
PROFILE OP AG OF THE AEROFOIL 


DENSITY OF ATMOS.PHFPF 

VOPTFX POTATION ABOUT REFERENCE SWEEP L I Ns 

CHANGE IN VOPTFX ROTATION ABOUT REFERENCE SWEEP LINE 


VOPTFX POTATION ABOUT S I TF LINE- 
CHANGE IN VOPTFX POTATION ABOUT S I T r LINE 
I NO I o A TEs NOMINAL VORTEX RANGE 
VQPTcry pam or 

UPPER VELOCITY ROUND ON PLOT I MG SCALE 
LOi'Fp VELOCITY POUND ON plot I mg SCALE 


r 

r 

r 

r 

r 

r 

C 


r 


■ r 


r 

r 


c 

r 

r 


MAMFLlFT ( SYSTEM ) 


iun i ts 


T CHO T r 


PS T 

PL Avn A 
NR 
PCS) 
NPLEMG 

plengt ( R ) 

WE I GmL 

I ARC 


cl Ac.ro 
PL Acer D 
FLActrD 


FLAG 


I UN ITS = 1 ONE LASERS 

HIM ITS = 2 T" f O LASFPS 

P*L A G 


I OHO I G = 1 CONTINOtjs wavF LASER 

I CHO I G = 7 PULSFD LASFP 

T NO T DENT A\/ c Ampl T TUDF 
WA\/F LENGTH os L I OHT 

INDICATES NOMINAL RADIUS OF TELESCOPE APERTURE 
RADIUS Os TELESCOPE APERTURE 
INDICATES NOMINAL PULSE LENGTH 


PUL c F LFMGTH 

DISTANGS BETWEEN SAMPLE POINTS USED IN WEIGHING 
FL A G 


I APS s 1 Ans(\/L) is ucpn 

I ABC - ? vL IS USED 

DICTAMGS OF L A s c*p c; ALONG S T TF L IMF 
DISTANCF OF LASFRC FROM F'^FH PLAMc 

DISTANCE OF LASFRC FROM S I TF LINE IN SWEEP PLANE 


NAME- L 1ST ( CO AN) 


r jninrr 


SLAG 
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r 


inlPrc = 1 

LATERAL SCAN 

r 


I 0 I P r C = P 

range scan 

r 

cr AN’i 

PC AM L r MATH 


r 

^CAm.d 

PC AM PPPOLUTTOM 


c 

Mc;t,fcrcrp 

INDICATES NOMINAL S'"FFP d 0F I T I ON 

r- 

Ci/c’pnD ( e; ) 

SWEEP POSITION 


r 

TN/MAV 

EL A. C 


r 


IVMAX = 1 

VMAX - VMIN IS PLOTFD 

r 


IVMAX = 2 

MAX(AF?S(VL)) IS PLOTFD 

r 


IV M AX = 3 


C 

T W F I CH 

FLAG 


r 


I WE I GH = 1 

PLOT HIGHEST ACTIVATED FILTER 

c 


IWEIGH = 2 

PLOT CENTROID 

r 

WF I LT 

BAND WIDTH OF El 

LTERs tn ME GAHFRTZ 

r 

nFTl_r~\/ 

DEJECTION level 

IN FILTERS 

r 

OF^Ftx ( P ) 

OFESFT OF LASER 

FOCUS POINT ALONG SITE LIME DO TO 

r 


N I SAL I GMMEnt 


r 

OF-’^Ftv ( P ) 

offset of laser 

FOCUS POINT FROM SCAN PLANE DO TO 

r 


m I c, A l_ I GN'Mc-mt 


C 

OF SETZ ( P ) 

0FFSET OF LASER 

FOCUS POINT ALONG SWEEP LINE DO TO 

r 


MISALIGNMENT 


r 

r 

W IN 0(3) 

CONSTANT '</ 1 NO 


r~ 

COMMON / AO I 

/ OP ( 3 ) , PO( 3 ) , D 

I SPLA ( 3 ) 


COMMON / A OP 

/ T I TL C ( 12) 



COMMON / A OP 

/ VI ( 100,4,2) * SCANP(IOO), not loo,^) 


COMMON / A Op 

/ OS(P), SR (3), VV ( 3 ) , VT ( 3 ) , AVT(3), OSS C 3 ) * VMOSS 


COMMON / A o/S 

/ GAMMAC, ’..<8, UNU 

« 7 


COMMON / A 07 

/ 91 , P2D 



COMMON / flOB / VM AX , VMIN, SMAX* .CM IN, NTHETA, NOM , J J , VAL 
COMMON / A14 / nzEPO, PHO 

common / A 1 5 / VL INI, VLIN2, VMVT, VMAVT « D, VTOL 

COMMON / A20 / VROTH, VROTHC* VROTV, VROTVC, nrange, VR ANGE ( 8 ) 

COMMON / API / IUNITS, ICHOIC, PSI, FLAMDA, NR, R<8), NPLENG « 

T. PLEMGT(P), "/EIGHL, iars, SLASER, DLASER* ELA5ER 
COMMON / A22 / IOIREC, SCANL, SCANR, NS WEEP, SWEEPP(5), IVMAX, 

*■ Twitch, wfilt, pftlev 

common / APO / OFSFTX(P), OFSETY ( 2 ) , OFSFTZ(P) 

OOMmom / API / W T NO ( 3 ) 

OIMFMSION FILTER (50)* VPV(4,P), D I SP ( 3 , ? ) 

OIMC-mETOM f t LTFN ( BO ) 

OAT A FILT^P / 50*0./ 

data GAMMAC* W 8 * UNU, Z* DZERO, RHO* VROTH, VROTHC* VROTV, VROTVC* 
I NPAmOF, VRANC-F, VMAX* v/mIn 

? / 7?9. , 77. P* .0914, 1000.* 50000.* 1.19P* O., 10.* 0.* 10.* 2* 

R BO . , 100.* POo., 1 OOP. , POOO.* 5000.* -1.* -1.* 34., 0. / 

DATA IUNITS* ICHOIC* PSI, FLAMDA* NR* R, NPLENG* PLENGT , WEIGHL* 

1 I Arc, SLASFp, n(_ ASFR , FLASFR 

P. / 1 , 1* 1 .* .0000 106* P, .1* .?* .3, .5, 1.* 2.* 4., 6.* 2* 

P 30., 90., 150., 300., 900., 3*— 1 • , . 5 * 1, O., o., 100. / 

DATA IDIREC, SCANL, SCANR, NSWEEP, SWEEPP, IVMAX, IWEIGH, WFILT, 

1 DETLEV / 1* PO., .5, 3, -15., -7., 0.* 7., 15., 2, 1, .94, .1 / 


HAT 

A 

OP / 9* 

o 

• 


/ 










PAT 

A 

( OESPTv 

( 

I 

) 

<i 

I 

= 1 

* 

p 

) 

/ 

p-if n 

• 

/ 

Z) 

> 

H 

A 

( OFSFTY 

( 

I 

) 


T 

= 1 

* 

P 

) 

/ 


• 

/ 

n at 

A 

( OF c FTZ 

( 

T 

) 

* 

I 

= 1 

• 

P 

) 

/ 

p#n 

• 

/ 

DAT 

A 

( W I ND ( T 

) 

* 

T 

= 

1 

< p 

) 


/ 

P-*0 • 

/ 
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RATA PI. RSO / 3.J41RP26. F7.2GS7R / 

NAMELIST / VORTEX / OAMMAC, Wfl. UNU . Z. D7FRO, RHO. VROTH. VROTHC , 

* VPOTV* VPRTVC. NPANGF, VPANGF, VMAX. VM I N 

NAMELIST / SYSTEM / IUNITS. ICHOIC. PS I . FLAMDA. NR. R. NPLENG. 

* PLFNGT. wfighl. iaps. SLASFP. olaser, flasep 

NAMELIST / SCAM / IDIPEC. SCANL. SCANR. NSWEEP. SWEEPP. I VMAX. 

T> I WEIGH. WEILT. DETLEV. OFSETX. OFSETY, OFSETZ. WIND 
1 A pf A n (F.foo) title 
F°o FORMAT ( 1 ?AA ) 

WRIT r(Aifno) TITL c 
60 A FORMAT ( 1 HI , 7RX . 1 PA6 ) 

PFAD ( F.vORTFX ) 

WP I TF ( 6 . VORTEX ) 

RF An ( F , SYSTEM ) 

'sfR I TP ( 6 . SYSTEM ) 

RFAD ( 5 , c C A N ) 

WRITF(ft.SCAN) 

NFW s 1 

FILT = FLA MPA #WF T LT* • OF ft 
00 1 1 10 = 1,00 
F !LT C N ( IF) = n. 

1 1 FILTER! IF) s O. 

TF(IARS .fo. 1) CO TO 30 
ASS I CM 67 TO JARS 
GO TO 31 

fo ASS I cm 6 C TO JARS 
31 T F ( I t>'F T OH .co. 1 ) GO TO 32 
ASSIGN DO to JWF I GH 
ASSIGN cjj jo KWFIGH 
GO TO "’R 

33 A G s I gn OR TO JWFIGH 
ASSIGN FO TO XWFIGH 
7 7 FONT T Mtjc- 

I F ( I UN I T c .co. 3) GO TO ? OR 
ASSIGN 2 OR TO JUNlTS 
ASSIGN 201 TO KUN ITS 
ASSIGN 111 TO LUMlTS 
ASSIGN 4 2R TO JR I GMT 
ASSIGN 432 TO WIGHT 
GO TG 200 

7PP. ASSIGN 2 a 6 TO JUNITS 
AG g joN 20 3 TO KlIM ITS 
Agc.Tgm 1 1 o TO LUN Its 
30 <o IF ( I P T RFC .FO. 3) GO TO 170 
ASSIGN 1F1 TO JDIRCC 
AGSIGN 1 FO to KD I P r C 
ASSIGN 1 c 7 to LOT pec 
ASSIGN 16 F TO MOT RFC 
ASSIGN 1^7 TO NOIPFC 
go to 1 7] 

1 70 aggIGM 1F2 TO JOlPco 
ASGIGM 1^4 TO KPT RFC 
ASGIGM ICO TO L n I Pec 
ASSIGN ] 6 G TO MO I PRC 
ASSIGN 1 FC TO NPTDFC 

171 IFIIGHOIC .FT. 1) GO TO 12 
ASSIGN 41 TO JCHOIG 
AGgIgn 44 TO KTCHOIC 
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OO T a 1 4 

1? ASSIGN 4? TO JCHOIO 
ASSIGN 4F TO XCHOTC 

14 OO TA ( IS, 1 6« 1 7) ' I WAX 
1 c ASS I ON °1 TO JVM AX 

OO TA Iff 

15 Accij^m op TO JVM AX 
OO TA 1 p 

17 ASS TOM 97 TO JVM AX 
IS NTHFTA = SCAMl./SCAMO + 1. 

IF(NTH C TA .LF. 1 a 0 ) GO TO 20 
WP ITS (A»610) 

610 FORMAT ( /43H *# INPUT ERROR ** SCANL/SCANR IS .GT. 100) 

OO TO 10 
P O OOMT I Ml JF 

t F ( 1 1 INI I TS ,FO. 1) OO TO 421 
niSP(M) = -oLAspo 
OlSP(2« 1 ) = -OLASFO 

O I SO ( p , 1 ) = -c|_ a SEP 

niS°( 1 ,?) = S|_ A CFP 

0 ISP (2,2) = OLASEP 

0 I so (7,2) = FLASFP 

IF( Ap.S( SLASER) + ABS(FLASER) ,LT. 1.) GO TO 420 

C P = SLASFP*sl asfp 

02 = OL ASSR-K-OLASFR 

F? = F!_A.SEP*FLAS c ‘P 

BRl = 2.*S2 + E2 

F>R2 = p.*SOPT(D2 + S? + F P) 

ASSIGN 497 TO JP I GHT 
ASSIGN 471 TO KPT GHT 
GO TO 491 

4 9 a ASSIGN 496 TO JR I GHT 
ASSIGN A 70 TO KR1GHT 
491 CONTINUE 

SC =-scANL/2. - SCAMP 
. no 9 f i=i* nthft a 
sc = SC + SCANP 
2F SCAnd ( I ) = SC 

SMAx = SCAMP (NTHFTA) 

S M I nj - SCANS* ( 1 ) 
no ipn J= l , i F 

GO To ( 2 1 0 « 990 ,9 70 , 30A , p 1 o , .720 » 770 , 340 , 7RO , 760 » 770 * 7SO » 790 » 4 00 « 
S 4 1 A ) , J 
91 a phi = n 
MOM = 1 

JJ = 1 

1 SWF = 9 

IF(MSWFFP .ms. 1) ISWP = 1 

OP( 1 ) = VPANGF ( NPAMGF ) 

PO ( 2 ) = 1. 

PO( 1 ) = TANK VPOTH/R2D) 

PO ( 7 ) = TANK VP0TV/P2D) 

I F ( I GHO I C ,F0. 9) GO TO 211 

PP - D ( NJP ) 

GO TA ~> n 2 

911 DLFMG = PLFNGT ( NPLFNG ) 

GO TA 14Q 
pp G I SWS = N| SWEEP 
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MOM = ? 

on TO )4P 
?? n IFwF = ? 

MOM = 1 

IF(MC.WFFP .LT. 3) TFWF = 3 
OO To 140 

? n O TF(MDfiMGF .OT. 1) OO TO 001 
OP( 1 ) = V/PANGF (NPAMGF + P) 

OO TO on? 

OO 1 OP(l) = VPANGF ( NPAMGF - 1) 

on? jj = ? 

VAL = 0P( 1 ) 

OO TO 73 

OlO vAL = VPANGF (NPAMGF) 

GO TO 1 30 

3?0 IFCNPANOF .FO. 8) GO TO 321 
0P( 1 ) = VPANGF ( NPANGF+1 ) 

IF ( Or>( 1 ) .OT. 0.) GO TO 32? 

021 OP( 1 i = VPANGF (NPANGP-?) 

02? VAL = 0P( 1 ) 

OO TO 73 

030 0P(1) = VPAMOF ( NPAMOF ) 

PO{ 1 ) = TAM((\/R0TH - VPOTHC ) /ROD ) 

JJ = 4 

VAL a VROTH - VPOTHC 
OO TO 73 
040 VAL a VPOTH 
OO TO 1 30 

330 PO( 1 ) = TAN( (VROTH + VPOTHC )/R2H) 

VAL = VPOTH + VPOTHC 
OO TO 80 

?60 PO{ 1 ) = TANK VP0TH/R2D) 

PQ ( 3 ) = TA N ( ( VPOTV - VPOTVO/R2D) 

JJ = F 

VAL a VROTV - VROTVC 
00 TO 30 
070 VAL = VPOTV 
00 to 1 30 

080 PO(0) a TAM ((VPOTV + VPOTVC ) /ROD ) 
VAL = VPOTV + VPOTVC 
GO To no 

700 PO ( 3 ) a TAN ( \/POTV/P?D) 

IF( lOHOIC .FQ. ?) CO TO 393 
TF(Mo . FO. 1) OO TO 391 
PP = o ( no- 1 ) 

OO To 093 
305 PR = R(NP+?) 

390 JJ = 6 

V AL = PP 

79? FACT 1 = (PM*P1*PP**?)**?/FLAMDA**? 

FACT? = FLAMDA/ ( P I *RR-K~*? ) 

I F ( J .FO. 1 ) OO TO 149 
CO TO 80 

ins T F ( NP>|_ FNO .FO. 1) OO TO 396 
PLFMC = PLFMOT ( NPLFNG-1 ) 

CO T o 797 

PLFMC = PL. FMGT ( MPI..FNO+? ) 

707 JJ = 7 
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VAL = PL p NG 

on to ra 

400 !F( Iz-HOIC ,RO. 2) GO TO 401 
VAL = P ( NR ) 

r-o TO 130 

401 v al S PLFNOT(MPLCMG) 

130 MOM - 3 

GO TO 130 

4 1 0 MFW = 3 

TE< IOHOIC ,FO, 2) GO TO 41F 
IFfMP ,FO, 3) GO TO 411 
RR = P < MR+ 1 ) 

1F(PD . GT. O.) GO TO 303 

411 PR = R ( NR— 3 ) 

GO TO 393 

41 « 1 F ( NPLENG .FO. 8) GO TO 416 

PLENA a PLFMGT ( MPLFNG+1) 

T F ( P(_FNG .GT. o.) 00 TO 397 
416 DLFMG = PLFMGT ( NPLFNG—? ) 

OO To 393 

140 VAL = ci,iFFPP( IRvir) 

GO To JOIPF-G , ( 161 , IF?) 

1 F 1 O I ST|_ = VPANG 17 " ( MRANGF ) + SWFFPP 1 TSWE) 

THFTA 1 a -ATAN(S 0 ANL/( 3 .*nTSTL)) 
OTHFTA = - 8 .*THFTA 1 /FLOAT ( NTHETA- 1 ) 

GO TO 86 

1^3 OISTO = VRANGF ( NPANGF ) 

THETA = SWFFPP < I SWF ) / OISTC 

GO TO 86 

3^ GO To MOTRO’O. ( 1 166) 

166 OIGTL = OP ( 1 ) + EK'FFPP ( NFWE c P) 

THETA 1 a -ATAN ( 30 AML/ ( 3. *0 ! STL ) ) 
OTHFTA =-3.*THFTA 1 /FLOAT (NTHFTA-1 ) 



GO TO 

BO 

1 ftft 

D I GTO 

= OP < 1 ) 


theta 

s BWFFPP ( NBWFFP ) /D I BTC 

80 

MOM = 

1 

8F 

I = 0 



GO TO 

KOIP FCi ( 183, 1B4) 

IBB 

THETA 

= THFTA1 - OTHFTA 


GO TO 

1 

1 B4 

OI STL 

= PT^TC - BC ANL* • B - BCANP 

1 B B 

DI STL 

= OIBTL + BCANP 


T = T 

+ 1 

1 B 8 

GO To 

JRTGHTi (426«427»A 28) 


PF’LT A 

= AT AM ( OLABFP/OT BTL ) 


GO TO 

4 OF 

/4P7 

Op - f> I STL*H I ^TL 


CGI = P3 + op + 03 + F3 

003 - -R#RR3 

CC-3 = R3 - 83 - 03 - F? 

CG4 = 3 • * 0 I STL 

4.33 CONTTMUF 

GO TO L.D T RFC » ( 1 R7t 1 RR) 

1R7 THETA = THETA+ OTHFTA 
T = T + 1 

168 GO to FRIGHT. (430,431 .433) 

431 DELTA! = ARGOS ( Rpi /( RB3*SQRT ( S? + ( FLASER-SL ASER*TAN < THETA )) **3 )) ) 
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AP = c.r 1 + r.C2*OOF<DFLTAl ) 

OFLTAR = ARC0FMA2 + Cr 3 ) / ( CC4* SORT ( A2 ) ) ) 

R2 = rr 1 + cc2*cof<pi - oeltam 
DFLTAA = ARCOF ( ( R2 + CC3)/(CC4*FORT<R2) ) ) 
nFLTA = DFLTAR + DFLTA4 
4 30 CRFLTA = FtNtOFLTA) 

CHFLTA = COF(RFLTA) 

TOFL.TA = Fnr-LTA/rnFLTA 
no 1 - frflta + crrlta/trflta 
arp no mo l 2 = 1 « 2 

00 TO JUNITF. (?OR,?Ofi) 

300 niFRL a ( 1 ) - 0. 
n I FPL A ( 2 ) = o. 

r>IFP|.A(R) = o. 

OO To P07 

206 0 I F P|_ A ( M = D I FP ( 1 tL? ) 

DIFP|_A(?) = DI.FP(2«L2) 
niFP| a(R) = D I FP ( 3 • L2 ) 

?07 V\/MAy s -1 . CRO 
VVM T m = 1 .RRO 

VLW = O. 
wft = 0 . 

OSF ( 1 ) = DI STL*COF ( THETA )*COF( PH I ) + DIFPLA(l) + 0FSETXCL2) 

0 F S ( p) - n I STL# C Os < THETA )*STN(PHI ) + DIFPLA(2) + 0FSFTY(L2) 

OFF ( o ) = D I FTL * F I N ( T HF T A ) + DIFPLA(3) + 0FSF.TZ(L2) 

V M OFc = VMAO ( OFF ) 

FO TO JCHOIC, (41 ,42) 

4 1 nniFTW = 2 .*FACT2*VM0FF**2 
OO TO 43 

42 nOIFTH = .F-ttPLFNO 
4 R NWEJOH = DO) T FT'n'/wF I OHL 
NH = 2*NU'F I OH + i 

OIFT'K s \/MOFF - ">F I GHL*FL0AT (NWF IGH + 1) 

TVLMy = o 
IVLMm = o 
no Po T I = 1 , ni><) 

OIFTv.i = DIFTH + "/FTGHL 
GO TO KCHOJCi ( 44 » 4E ) 

44 DL2 = ( F ACT2*\/v0FF-)(-0 I FJW ) **2 

WF = F ACT 1 / ( 1 • + (VM0FS - D T FTW ) ** 2/DL2 ) /D 1 FT W**2 

OO To 40 
4^ "/AT = 1 . 

I F ( M . 0T . NWF I OH ) GOTO 60 

"IF = WAT*FLOAT ( I I -1 ) /FLOAT ( NWF T OH) 

00 To 46 

60 i,/F = "(AT-*f|_oAt ( N"'— I I ) /FLOAT ( N",IF I GH ) 

46 CALL VFL(VL»PHI , THFTA .O I FT"’) 

1 F ( VL . OT. VVM AX) \/VM Ay = VL 
I F ( VL • LT • \/VM IN) VVM IN = VL 
GO TO JARF, (6^,67) 

6R VL = ARF(VL) 

67 I F ( T T • NF . NWFIGH + 1) GO TO R7 

VI (I .3. MOM) = VTOL 
\/2v ( 4 « L 2 ) = VL 

no ( t ' non ) = n 

R7 "IF y - WFT + WF 

OO TO Kl.ic T GH . ( FO» F 1 ) 

R 1 VL"' = VL"’ + VI... *"/*=■ 
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GO T O on 

I F ( VL • LT • O.) GO TO tOF? 

IVL = VL/FU.T + 1. 

FILTFP(IVL) = WF + FILTFP(IVL) 

T F { I \/L • GT . TVLMy) TvLMX = IVL 

GO TO on 

ino.n ivl a -VL/FILT + J. 

FILTP-M(IVL) = WF + FTLTFM(TVL) 

I F ( I \/L . GT. TVL M N> TVL.MN = T VL 

°n gonttmuf 

0-0 TO Jw/cr I ON . ( OF , QO ) 

V?V( 1 *L?) = VLW/WFT 

GO TO GO 
OF 17 = 1 
177 = 1 

IF< IVlJVy .FO. O) GO TO 1 OOF 
no o/, I 6 =I.I\/|_MX 

FILTrO( 1^) = FILTPPI T 6 )/WFT 

04 IF( F TLTFR ( 16) .GT. OFTLFV) 17 = 16 

IF<IWFTGH .no. 1 ) WRITF( 6.670) ( F I LTFR ( I P. ) . I 8 = 1 . I VLMX ) 

67? FORMAT (//noc]?,4) ) 
loop TF(TVI_MN .FO, n ) GO TO of 
po 1 no4 I A - 1 , I VL M N 

F T LTFM (16) = F TLTFM ( I 6 ) /wFT 

1004 TF( FTL.TFNI 16) . GT . OFTL^V) 177 = 16 

IF< IwFIGH .FO. t) WRITn( 6< 16.70) (FlLTFN( 18) < 18=1 » I VLMN) 

167? FORMAT (//(OH -»FH,4,ofi?,4)) 

°F IF( IT7 - 17) 1007,100.0,1001 

100? IF(F 1 LTFP( 17) .OF. F I LTFM ( ! 77 ) ) GO TO 1003 
1.001 \/?V(1.l ..P) = -F ILT* ( FLOAT ( 177) - .n. ) 

GO TO 1 004 

1007 V?V ( 1 iL?) = FILT*(FLOAT( 17) - . F ) 

1004 GONTTMUF 

no 400 IO = 1 , 1 VLMX 

40? F 1 LTFD ( TO) - n, 

no 1400 1 0-1, IVL MM 

14 0? F ILT<“N ( I o ) - n . 

Pi o (70 TO JVMAX , ( o 1 , o? , 07 ) 

Q 1 V?V(7,L3) = VVMAX - V/VMIM 

GO To 07 

O p \/?V(?,LP) = AMAX1 (VVMAX, -VVMIN) 

GO TO 07 
07 GONTTMUF 
o-7 C^NTjN'Mr 

r.n to u iMI T^«(tll*ll°) 

1 1 o rONT T KM )F 

in GO to in INI T ( PC) , pop) 
op i no pm K f = ] # 4 

IF( Kk .FO. O) GO T O pop 
VI ( 

pm roNTiNi f 
GO T o 004 
pop no l i f i<f= i , 4 

TF(l^ .FO. 3) GO TO 11.F 

VI < I = SORT ( VPV ( «5 * 1 ) + ( - V2V ( KF < 1 ) /TOELTA + V 2 V(KF,p>) 

*t. of* po i ) v- P ) 

I 1 F FONT T Ml IF 

T F ( MOM .FO. p ) WP I TP ( 6 , ^PO ) VPV ( 4< 1 ) *\/?V ( 4 « ? ) « VI < I • 4 *N!OM > 
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f~-?r poomat 

CONTT^IUF 

inn i f ( n X h c t a .ft. m no to ndtpfc* ( 157. iff ) 

1 CALL PV(NFW) 

F -0 TO in 
FMD 


C-34 



LMSC-HREC D225936 


i p f t c vfl nrcw' 

FUBPO( IT I NF VFL ( VL. INF « PhTLi thfTA.L * D I FTL ) 

COMMON / A n 1 / OP ( 0 ) , PO ( 3 ) « OISPLA(O) 

COMMON / A OF / OS<3), FP ( 3 ) » VV<3)« VT ( 3 ) i AVT (3), OSS ( 3 ) ♦ VMOSS 


COMMON / AOfi / OAMMAC* ’<'F» UNU » 3 

COMMON / A 07 / Pi, POD 
COMMON / A14 / D 7 FPO, PHO 

COMMON / A 1 F / VL INI* VLIN2« VMvTt VMAVT * D» 
COMMON / A31 / WIND (.3) 

DTMFmfTON TOj ( 3 ) 

\/MOS Fp = niFTL/VMOFF 
OF ( 1 ) = VMOCFd*0 c S( 1 ) 

OF ( p ) = VMOFFO*OFF ( P ) 

OF ( 3 i = vnpc3d*o<= F ( 3 i 

°.p ( i ) = open - oF<n + oifpLa(I) 

FP(3) = 00(3) - OF (3) + n I FOL A ( 3 ) 

FP ( 3 ) = OP ( 3 ) - OF ( 3 ) + O I COLA (3) 

CALL OPOCF(FP,no,VV/) 

VMPO = VMAO(PO) 

VMVV = VMAC, (VV) 
n = \/m\/\//\/mdo 

FX = FXP ( -'■.IF*D**3/ C A • *l JNU*7 ) ) 

VMVT = GAMMAC* ( 1 ■ - FX ) / ( 3 . *P I *n ) 

VMVTO = VMVT/VMVV 
VT ( 1 ) =VMVTP*VV ( 1 ) 

VT ( 3 ) =VMVTP*VV(3) 

VT ( 3 ) =\/MVTP*VV < 3 ) 

VMAvt = 07FP0*FX/( 4 • * P I *RHO*UNU*Z ) 

VMAVTO = V v AVT/\/mPO 

A VT ( 1 ) = VMAVTD*oo(]) 

A VT ( 3 ) = V M AVTO#PO ( P ) 

AVT ( 3 ) = VMAVTD*P0(3) 

VMOF = VMAC (OF) 

PO 1 o 1=1,3 



10 TOT ( l ) = VT ( I ) + AVT ( I ) + W I NO ( I ) 

VTOL = FORT ( TOT ( 1 ) **3 + TOTC?)**3 + TOT (3)**?) 
VLINFI = -OOT ( VT ,OF ) /VMOF 


VLINF3 = -DOT ( AVT » OF ) /VMOF 


VLINfi.i = -DOT Oi/I MD , of ) /VMOF 
VLINf = VL. t NF 1 + VL I NF 3 + VLIMFi,i 

PFTl Ion 
FND 
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TPFrr ny prrx 

C, unpot IT I NF o\/ ( NFW ) 

COMMON / A OP / TTTLE<13) 

A03 / VI ( 100*4*?) , SC ANP ( 1 00 ) « DO ( 100*3) 

AOR / \/m AX , VM IN, SM AX , EMIN, NTHFT A * NOM « JJ, VAL 
A?? / IDIREC, SCANL, SCANR, NSWEEP, SWEEPP<5>» IVMAX, 
WFILTf aftLFV 

RCDX(13*3), <=<CDY(13), ISYMC4), L ARLF (8,7) , SFMRLE ( 6 ) * 


COMMON / 
COMMON / 
COMMON / 

* I wr T OH , 

otmfmsion 


T> RFM(R,R) , CFMP(R,?1 
DATA ( RCDX (1,1), 1 = 1, 12) / 

* 6H PFSPF * 6HCT TO « 6HS I TF 
DATA ( RCDX (1,3), 1=1, 12) / 


6HDATA P , 6H0 I NT P , 6H0S I T I 0 * 6HN WITH* 
L , 6H I NF IN,6H MFTFP , 6HS , 2* 1 H / 

6HDATA P , 6H0 I NT P , 6H0S I T I 0 , 6HN WITH, 


T 6H PFSPF , 6HCT 
S 1 H / 


TO , 6HREFERE,6HCE SWE,6HEP LIN,6HE IN M , 6HETERS , 


/ 


data 

( BCny ( 1 ) 

,1=1,12) / 

6hveloc 1 , 

6HTY IN 

, 6HM/SEC 

, 9*1H 

DATA 

( I ^VM ( I ) 

,1=1,4) / 

IT 

* 

U 

IT 

#■ 

C 

IT 

« 38 / 



DATA 

( L ABLF ( I 

, 1 ) , I = 1 ,B) 

/6 H 

♦ 6H 

9 6H 

, 6H 

T 6H 

9 SH 

SWEEP , 6H POS I T , 6H I ON 

IS / 



HAT A 

( L APL r ( I 

f P ) * 1 = 1 % 8 ) 

/ 6W 

,6H 

«6H 

,6H 

S *H : 

0 ISTA, 6 HNCR 00 , 6HWN 

VOP, 6HTFX 

IS / 



DATA 

(LAHLF( I 

, 3 ) , I = 1 , 8) 

/ 6H 

« 6H 

, 6H 

,6H 

S 6H 

♦ SH 

VOP , 6HTFy PA, 6HNGF 

IS / 




DATA (LARLEI I ,4) , 1=1 ,8) / 6H 


VOR , 6HTEX RO, 6HTAT I ON , 6H ABOUT 


T. 6H 

DFFFP,6H C NCF S« 

6HWFEP 

L,6HINE 

IS / 




DATA 

( L ARLF ( I ,5) , I = 

1 , B ) 

/ 

6H 

,6H 


,6H VORT 

♦6HFX POT 

R 6HAt I ON , 6HAR0I IT , 

6HSITF 

L , 6H I ME 

T S / 




DATA 

( L ABL.F ( I ,6) , 1 = 

1 « 8 ) 

/ 

6H 

, 6H 


, 6H 

,6H 

R 6H 

T , 6HFLESC0 , 

6H D F 

PAD , 6H I US 

IS / 




DATA 

(LARLEI I ,7) , I = 

1 ,8) 

/ 

6H 

« 6H 


, 6H 

• 6H 

T 6H 

, 6'H Pi IL , 

6HSF 

LFN , 6H0TH 

IS / 




DATA 

( SEMPLF ( I ) , I = 1 

« 6 ) 

/ 

6HT = NO 

,6HPM 

TOT 

,6HAL L 

,6H= LINE 

0 6H 

OF SI, 6HOHT 

/ 







DATA 

( SEM ( I , 1 ) , I = 1 , 

3 ) 

/ 

6HX = VM 

, 6HAX 

- V 

, 6HMIN 

/ 

DATA 

( SFM ( I ,3 ) , I = 1 , 

3 ) 

/ 

6HX = MA 

, 6HX( ARS ( 

, 6HVL ) ) 

/ 

data 

( SEMR ( I , 1 ) , I = 1 

* 5 ) 

/ 

6HW = HI , 

6HGH 

ACT , 

6H I VATED , 

6H FILTF« 

DAT A 

( SFMP ( I , 3 ) , I = 1 

* F ) 

/ 

6HW = CE, 

6HNTRCIOi 

3*1 H / 



1 HR / 


1 r 


?! 


oo to ( i o , 20 , pn ) ,nfw 
new = ? 

CALL CAMOA\/(0) 

call p input 

Nl_ = P 
l_ = -0 

IRC ML ,FO. P) L = P 
CALL OU I K3L (L. , S.M IN, SMAX , VM IN, VM AX , ISYM ( \ ) ,RCDX< 1 , I D IREC ) , BCDY , 
T MTHFTA , SCAMP , V 1 ( 1 , 1 , NOM ) ) 
oo pn Hs ?,4 

1 no CALL QUIKRL ( 0 ,DUM,DUM,DUM,DUM, I SYM ( I 1 ) , DUM , DUM , NTHET A , SC ANP , 

* V 1 ( 1 , I 1 , NOM ) ) 

CO TO ( 11 O , 1 ?o , ] SO ) ,NL 

CALL PR INTV (7?, TITLE, 37,1008) 

CALL PRINTV(48,LARLR( 1 , JJ) , PA , 33^) 

CALL LABLVIVAL, 4 74, 32R, -4,1,4) 

CO TO 1 RO 

CALL OP I MTV ( 48, L ABLE ( 1 , JJ) , 34,860) 

CALL LARLV(VAl_, 4 34, 660, -4 ,1,4) 

CO Tp 1 RO 

CALL °P I NT V (48 , L ARL.F ( 1 , J J ) , 34, 993) 

CALI. LARLVIVAL, 434, 903, -4,1,4) 
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1 


1 RR 


PALL or>iMTv/(m,^^MRLF,^?p, Q a?) 

CALL PPIMTV(? c i*PFMB( 1 , I W^IGH) «80R, 90?) 

CALL. PR I MTV ( lA.F-FMI 1 * I\/MAX) *800. lODR) 

Ml. = ML - 1 

!F(N|_ *NF. O) CO TO IFF 
ML = -> 

L = ' -? 

!F(N=-W .FO. 0) CALL. CLP AM 

PFTUON 

FND 
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f! n FTr P T NOI IT nFCK 

FI IRROI IT IMF PJNOUT 
COMMON / AO? / T I TLF ( 1 ? ) 

COMMON / AOP / GAMMAC, Wfl » UNU* Z 

COMMON / AOP / VMAX * VMIN, SMAX, PM IN* NTHFTA* NOM* J J , VAL 
COMMON / A14 / D7FP0, PHO 

COMMON / A20 / VPOTH* VROTHC * VROTV* VROTVC* NRANGE ♦ VRANGEC8) 
COMMON / A?1 / I UN I TP* I CHOI C* PS I* FLAMDA , NR* R(8)« NPLFNG* 
p plengtir). weighl, tabs. PLAPFR* DLAPFR, FLAPFR 
COMMON / A22 / I O I RFC « SC ANL « SCANR * NSWFEP * SWEEPP ( 5 > * IVMAX* 

* I WEIGH* WFILT* n?TLFV 

COMMON / A 30 / OFPFTX ( 2 ) * 0EPETY(2)* 0ESFTZI2) 

COMMON / A 3 1 / W I NO ( 3 ) 

CALL eCOI TV(Q ) 

WR I TP ( 1 P , POO ) 

WP I TP (IP* POO ) T I TLF 
POO FORMAT (1H1.24X.12AP) 

WRITf( 16* pi oigammac* wp* unu* z* dzero. rho* VROTH* VROTHC* 
p VROTV* VROTVC* NRANGE* VRANGF* VMAX* VMIN 
610 FORMAT ( //B 1 X * 1 7HNAMEL I ST ( VORTEX >*// 1 X , 1 1 HGAMMAC. =*E16.8*3X* 

1 llHwR =*F1P.8,3X* 1 1HUNU = * E 1 6. 8 , 3X « 1 1 HZ =* 

2 Fip.R,/1?H D7FR0 = , F 1 6. 8 ♦ 3X * 1 1 HRHO =.F16.8*3X, 

3 11HVPOTH =«F1P.8«3X. 1 1HVP0THC =*E16.8/ 

A 1 2H VROTV =«F16.8,3X, 1 1HVR0TVC = . F 1 6 . 8 , 3X . 1 1 HNR ANGE =. 

8 18, / 

6 12H VRANGF ( 1 ) = , E 1 6 • 8 , 3X * 11 HVRANGE ( 2 ) = « E 1 6. 8 , 3X * 1 1 HVR ANGE < 3 ) =, 

7 

8 3X, i 1 HVRANGE (4) =*E16.8*/12H VRANGE ( 5 ) = * E 1 6 .8 * .3X * 1 1 HVRANGE < 6 ) =* 

9 F16.8.3X. 1 1HVRANGFI 7) = < E 1 6 • 8* 3X ♦ 1 1 HVRANGE ( 8 ) =*F16.8*/ 

A 12H VMAX = ,E1 P.8, 3X* 1 1HVM IN =.E16.8) 

WRITF( 16*620) I UN ITS. ICHOIC. PS I. ELAMDA « MR* R* NPLENG * PLENGT. 

P WEIPHL. I ARP. PLA3FR. DLAPFR . FLAPFR 


FORMAT <//rjx. 

1-7HNAMFLIST 

( SYSTFM) // 


7 

1 2H IUNITF 

= 

.18 / 



1 

12H ICHOIC 


. 18, 1 IX . 

1 1 HPS I 

= ♦ FI 6.84 3X4 

p 

1 1 HFLAMO A 

= 

, E 1 6 • 8 , 3X « 

1 1 HNR 

= ♦18 / 

7 

1 2H R( 1) 

= 

*E1 6.8, 3X* 

1 1 HR C 2 ) 

=.F16.8.3X4 

4 

1 1 WR ( 3 ) 

= 

* E 1 P.8, 3X * 

1 1 HR ( 4 ) 

=♦816.84/ 

rr; 

1 2h R(R) 

= 

* F 1 6.8, 3X. 

1 1HR( 6) 

=§F16.8. 3X « 

6 

1 1 HR ( 7 ) 

= 

, F 1 6,8, 3X. 

1 1 HR ( 8 ) 

=.F16.8./ 

7 

1 2H nplfnc. 

- 

,18 / 



R 

1 2H PLFNGT ( 1 ) 

= 

* F 1 p , 8 * 3X , 

1 1 HPLFNGT ( 2 ) 

=%F16 #8*3Xf 

n 

1 1 HPL.FNGT ( 3 ) 

= 

* F 1 P • 8 » 3X , 

1 1 HPLFNGT ( 4 ) 

=♦816.84/ 

A 

1 2H PLFNGT(R) 

= 

,F 1 6 . 8 , 3X * 

1 1 HPLFNGT ( 6 ) 

=♦816.84 3 Xf 

n 

1 1 HPLFNGT ( 7 ) 

= 

,F1 6.8, 3x. 

1 1 HPLFNGT ( 8 ) 

=♦816.84/ 

r 

12H WFIC-HL 


, F 1 6 , 8 , 3X * 

1 1 HI ARp 

= ♦18411X4 

n 

1 1 HPLAPFR 


, F 1 16 . 8 , 3X *. 

1 1 HOLASFR 

= ♦ F ] 6 . 84 / 

pr 

12H FLAPFR 

= 

,F 1 6.8) 



WR I TF ( 1 P , PRO ) 

IOIRFC. PC ANL . SCANR, NSWFFP, SWFE pp, IVMAX* 


lu/FTGHt u»FTLT 

4 

OFTLFV, OFSFTX* OFSFTY, 

OFRFTZ ♦ WIND 


PRO FORMAT ( //R 1 X » 1 7HNA m FL I PT (PCAN) // 

1 1 2H IOIRFC = , 1 8 « 1 1 X .11H5CANL =*E1P.8*3X. 

? 11HFCANP = »E 1 6. 8, 3X * 1 1 HNSWEEP =.I8,/ 

3 1 2H PWEFPP ( 1 ) =,E16.8, 3X * 1 1 HSWEEPP ( 2 ) =*F16.R*3X. 

/■ 1 1 HSWEEPP ( 3 ) = , F 1 6 • 8 * 3X *11 HSWEEPP ( 4 ) =*F16.8,/ 

B 12H FWFFPP(B) = * F 1 P • 8 / 

P 1 2H IVMAX =« 18*1 IX* 11HI WEIGH =, 18,1 lx , 

7 11HWFTLT = ,E 1 6 • 8, 3X * 1 1 HDFTLFV =.F1P.8 / 
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8 12H OFSFTX(l) = *E 1 6. 8, 3X « 1 1 HOFSFTX ( 2 ) =«F16.8 

9 1 1 HOFFFT Y ( 1 ) =,F16.8,3X« 1 1H0FFFTV(?) =,F16.8 

A 12H OFSETZ(l) = ♦ F 1 6 • 8 , 3X » 1 1 HOFSETZ ( 2 ) =«E16.8 
8 12H WIND( 1 ) = ,E16.8,3X» 1 1 HW I ND < 2) =*E16.8 

C 1 1 HW I NO ( 8 ) = . E 1 6 . 8 ) 

PPTIPN 

FND 


3Xt 

/ 

/ 

3X« 
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SlRFTr CROCS nFCK 


SUBROUTINE CROSS ( VA ♦ VR . VC ) 

C FUBPOl .iT I NF TO PFRFORM CROSS PRODUCT 

0 VFCTOP V A I S CROSSED INTO VFCTOP VR TO FORM VFCTOP VC 

DIMFmcioN VA(R) . \/R<3) «VC<3) 

VC( 1 UVA(?)#VR(F)-VP(?)*\/An) 

VC ( P ) =V A { 3 ) *VP ( 1 ) -VR ( 3 ) * VA ( 1 ) 

\/C < 3 ) = V A ( ) ) *VR ( 2 )-VR( 1 ) *VA ( ?) 

RFTIIOM 

PND 

T I nnr DOT DFCK 

FLJNCT I ON DOT < V A . VR ) 

C FUNCTION TO PERFORM DOT PRODUCT 


CROOP 1 ' 
C R 0 0 o 3 1 
CR0DO3( 
CR0004' 
CROnoF' 

CROOOfti 

CR0.007 1 
CP0008' 
CPOOOP( 


DOTOO 1 < 
D0T0O2( 


DIMFMSION VA ( 3 ) . VP < 3 ) D0T0O4' 

O0T=VA ( 1 ) *VR ( 1 ) +VA ( 2) *VP ( 2 ) + VA ( 3 ) *VP ( 3 ) POTOORt 

C VECTOR VA IS DOTTED INTO VECTOR VB TO VORM THE SCALAR DOT PRODUCT DOT003< 


pftuon 


D0TOO6( 


C M0 DOT007' 

T I PFTO VMAO DFOO 

FUNCTION VMAG(VA) VMA^ op 

C FUNCTION TO DFTFRMIMF VFCTOR MAGNITUDE VMAno^f 

DIMfmrION VA ( 3 ) V M A 0 0 3 ( 

VM.Ac = SORT ( VA ( 1 ) *VA ( 1) +VA ( 2 )*VA ( P ) +VA ( 3) *VA ( 3 ) > yMACOAC 

PFTUPN \/MAOoe-,< 

FND \/m A’OORf 


TPATA 

C ASF No' 24 1 LASFR RANGF 300 LATERAL SCAN FILTERS 

T-VORTFX O7FR0 = O., VROTHC = 30.. VROTVC = 30.. NRANGE = 3. $ 

<r. cystpm « 

RFC AN SWFFPP = -7., -2., R 
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